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Einige Eigenschaften von Lésungen 
der Prandtlschen Grenzschicht-Differentialgleichungen™ 


KARL NICKEL 


Vorgelegt von H. GORTLER 


Ubersicht 


Seit der Begriindung der Grenzschichttheorie (GS-Theorie) durch L. PRANDTL 
im Jahre 1904 [J] sind eine sehr groBe Anzahl von Arbeiten erschienen mit dem 
Ziel, die Prandtlschen Grenzschicht-Differentialgleichungen (GS-DGlen) zu lésen 
(vgl. dazu etwa das Buch von H. ScHLicHTING [2]). Man kennt dadurch heute 
schon eine groBe Zahl spezieller Lésungen der GS-DGlen. Dagegen scheint bisher 
noch kaum der Versuch gemacht worden zu sein, schon ohne Integration allge- 
meine Aussagen tiber die mégliche Gestalt einer Lésung der GS-DGlen zu ge- 
~winnen. Diese Aufgabenstellung soll der folgenden Arbeit zugrunde liegen. Um 
den Umfang nicht zu sehr anwachsen zu lassen, sei dabei nur der Fall der zwei- 
dimensionalen stationaren laminaren Wand-GS in einem inkompressiblen Me- 
dium behandelt. 

Die Arbeit gliedert sich in zwei Hauptteile. Zur Erlauterung der Aufgaben- 
stellung werden in einer vorangestellten Einleitung in heuristischer Form ver- 
schiedene UnregelmaBigkeiten von GS-Profilen geschildert. Die angestrebten 
allgemeinen Gestaltaussagen lassen sich aus den GS-DGlen mit gewissen Ab- 
schatzungs- und Einzigkeitssatzen gewinnen, die aus einem Abschatzungssatz 
von M. Nacumo [3] und H. Westruar [4] folgent. Da die genannten Ab- 
schatzungs- und Einzigkeitssatze auch unabhangig von den speziellen GS-DGlen 
von Interessen sind, werden sie im ersten Hauptteil I fiir eine allgemeine Klasse 
von partiellen expliziten parabolischen DGlen ausgesprochen. Die Anwendung 
der Ergebnisse dieses ersten Teils auf die GS-DGlen fiihrt dann im zweiten Haupt- 
teil II zu den gewiinschten Aussagen iiber die médgliche Form eines GS-Profils. 


Einleitung. Einige spezielle Eigenschaften von laminaren Grenzschichten. 
Heuristische Betrachtungen 


In dieser Einleitung soll versucht werden, anschaulich auf gewisse Eigen- 
schaften von Lésungen der GS-DGlen hinzuweisen, die spater im Teil II streng 


* Von der Fakultat fiir Natur- und Geisteswissenschaften der Technischen Hoch- 
schule Fridericiana Karlsruhe als Habilitationsschrift angenommen. Fiir viele kri- 
tische Bemerkungen und anregende Diskussionen danke ich Herrn Prof. Dr. JOHANNES 
WEISSINGER und Herrn Dr. WoLFGANG WALTER recht herzlich. 

1 Auf die Bedeutung des Nagumo-Westphalschen Satzes fiir die GS-DGlen hat 
zuerst H. GORTLER [5] hingewiesen. 
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bewiesen werden sollen. Gleichzeitig soll an Beispielen demonstriert werden, daB 
diese Eigenschaften nicht selbstverstandlich sind, so daB der mathematische 
Beweis ihres Erfiilltseins also nicht durch den Hinweis auf eine ,,anschauliche 
Evidenz‘‘ ersetzt werden kann. Unter ,,GS‘‘ sei dabei in der ganzen folgenden 
Arbeit immer eine laminare GS verstanden, auf diese Eigenschaft der Laminaritat 
wird spater nicht mehr eigens hingewiesen werden. 


1. Ubergeschwindigheiten 
Ein inkompressibles Medium stréme stationar, zweidimensional und drehungs- 
frei in x-Richtung entlang einer ruhenden und undurchlassigen Wand mit einer 
Geschwindigkeit U(x) und bilde dabei ein GS-Profil w(x, y) aus (vgl. Abb. 1). 


Abb. 1. Grenzschichtprofil Abb. 2. Grenzschichtprofil mit Ubergeschwindigkeiten 


Dann ist aus physikalischen Griinden zu vermuten, daB diese GS stromabwarts, 
d.h. in Richtung wachsender x-Werte, keine Stellen mit , Ubergeschwindig- 
keiten‘“’ u(x, y)> U(x) (vgl. Abb. 2) ausbilden kann. Dieselbe Eigenschaft wird 
man dann auch von der mathematischen Lésung der GS-DGlen erwarten. Dieser 


CM AO ee A 
Ole 
Avsblasen 

Abb. 3 Abb. 4 


Abb. 3 und 4, Erzeugung von Grenzschichtprofilen mit Ubergeschwindigkeiten; durch mitbewegte Wand (Abb. 3), 
durch Ausblasen (Abb. 4) 


Sachverhalt wird sich in ITI.2 streng beweisen lassen. Dabei sind die oben for- 
mulierten Voraussetzungen liber die Art der betrachteten GS durchaus wesentlich, 
denn in einer ganzen Reihe von verwandten Fallen ist diese Vermutung falsch: 


Betrachtet man namlich zuniachst dieselbe Strémung mit dem einzigen Unter- 
schied einer bewegten Wand, dann ist anschaulich klar, da8 hierbei Uber- 
geschwindigkeiten auftreten kénnen (vgl. Abb. 3). Dasselbe gilt offenbar auch 
wenn schrag zur Wandoberfliche abgesaugt oder ausgeblasen wird (vgl. Abb. 4). 

Verzichtet man nun auf die Eigenschaft der Inkompressibilitat des strémenden 
Mediums, dann kénnen schon bei ruhender und undurchlassiger Wand (bzw. 
pordser Wand mit senkrechtem Absaugen oder Ausblasen) Ubergeschwindigkeiten 
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entstehen, wie zuerst W. B. Brown & P.L. Donoucu [6] gezeigt haben. Ein 
solcher Fall tritt z.B. ein, wenn die Wand hoch aufgeheizt wird und strom- 
abwarts ein starker Druckabfall vorhanden ist. 


Aber auch bei inkompressiblen Medien kénnen Ubergeschwindigkeiten von 
selbst entstehen, wenn nur die Strémung instationar ist. Man vgl. dazu etwa 
die Funktionen ¢) und ¢;, in Bild 11.1 bei H. ScHiicutin [2], Seite 170, aus 
denen sich sofort GS-Profile mit Ubergeschwindigkeiten aufbauen lassen. 


SchlieBlich haben A. MacEerR & A. G. HANSEN [7], A. MaGER [8], H. G. Loos 
[9], H. Z. Herzic & A. G. HANSEN [10] und andere gezeigt, daB selbst bei sta- 
tionarer Str6mung in einem inkompressiblen Medium bei ruhender und undurch- 
lassiger Wand Ubergeschwindigkeiten entstehen kénnen, falls die Bewegung 
dreidimensional erfolgt und die AuBenstrémung nicht drehungsfrei ist. 


2. Oszillationen des Grenzschichtprofils 


Von nun an werde die Strémung als inkompressibel und stationar voraus- 
gesetzt, die Wand sei ruhend und entweder undurchlassig, oder, wenn Absaugen 
oder Ausblasen eintritt, dann mége es wenigstens stets senkrecht zur Wand ge- 
schehen. Sowenig es physikalisch zu ver- 
muten ist, daB in einer zweidimensionalen 
GS unter diesen Voraussetzungen von selbst 
Ubergeschwindigkeiten entstehen, ebenso- 
wenig ist mit dem Auftreten von noch 
, noheren“ UnregelmaBigkeiten, z.B. Oszilla- 
tionen (vgl. Abb. 5) zu rechnen. Es ist viel- 
mehr anschaulich evident, daB unter der 
Wirkung der Zahigkeitskrafte sich jede 
Stérung stromabwarts beruhigen wird, so u(x,y) 
da8 also ein GS-Profil der Gestalt von Abb. 5 Abb. §. Grenzschichtprofil mit Oszillationen 
stromabwarts schlieBlich in eines der Gestalt 
nach Abb. 1 oder Abb. 2 tibergehen wird. Der Beweis dieser Vermutung wird 
in II.3 gebracht werden. DaB diese Uberlegung jedoch bei dreidimensionaler 
Strémung nicht mehr richtig zu sein braucht (wenigstens nicht bei drehungs- 
behafteter AuBenstrémung), konnte kiirzlich gezeigt werden [/7]. Dabei handelt 
es sich sogar um ein Beispiel einer exakten Lésung der Navier-Stokesschen 
DGlen, so daB die Ursache fiir die auftretenden Oszillationen also nicht in den- 
jenigen Vernachlassigungen liegen kann, die beim Ableiten der GS-DGlen aus 
den vollen Navier-Stokesschen DGlen getroffen werden. 


uf 


3. Strémungsablosung 


Im folgenden sei die Strémung sogar zweidimensional und die AuBenstrémung 
drehungsfrei. Da8 dann ein hinreichend groBer Druckanstieg in Strémungs- 
richtung eine Strémungsablésung hervorrufen wird, ist schon lange bekannt 
(vel. etwa H. ScHLICHTING [2]). Wie jedoch aus einer Arbeit von H. GORTLER [12] 
hervorgeht, wird eine Strémungsablosung auch schon eintreten, wenn einer kon- 
stanten AuBengeschwindigkeit eine periodische Stérung mit beliebig kleiner 
Amplitude iiberlagert ist, falls nur hinreichend viele Schwingungen tiberstromt 
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wurden. Dagegen wird sich beweisen lassen (II. 4), daB ohne solche Stoérungen 
bei konstanter oder zunehmender AuBengeschwindigkeit das Phanomen der 
Strémungsabliésung niemals eintreten kann. 


4. Wendepunktprofile 
Dagegen ist durchaus zu erwarten, daB ein Grenzschichtprofil sich eine groBere 
Anzahl von Wendepunkten erwerben kann. Denn da nach Abb. 6 der Ubergang 
von Druckabfall zu Druckanstieg (in Abb. 6 durch die Gestalt der Wandkontur 
angedeutet) Grenzschichtprofile mit einem Wendepunkt erzeugt und sich dies 
bei jedem Vorzeichenwechsel des Druckgradienten wiederholt, ist zu vermuten, 
daB ein GS-Profil an einer welligen Wand um so mehr Wendepunkte besitzen 
wird, je mehr Wellungen schon iiberstromt wurden. DaB dies in der Tat zutrifft, 
sieht man aus Rechnungen von A.W. Quick & K. SCHRODER [13], bei denen 
schon nach der ersten Welle 
Profile mit zwei Wendepunkten 
auftreten. (DaB nach den weite- 
ren Wellen das zu vermutende 
weitere Anwachsen der Wende- 
punktanzahl nicht mehr zu be- 
Wendepunk! obachten ist, liegt daran, daB 
WU die GS-Profile durch ein Diffe- 
Abb. 6. Entstehung eines Wendepunktes im Grenzschichtprofil renzenverfahren berechnet und 
bei Druckanstieg aye 6 

also nur an einigen Punkten 
bestimmt wurden. Damit kann natiirlich die zweite Ableitung des GS-Profils 
nicht mehr genau festgelegt werden, wodurch die Beobachtung von nahe bei- 
sammenliegenden Wendepunkten stark erschwert ist.) In II.5 wird die héchst- 

mogliche Anzahl von Wendepunkten abgeschatzt werden. 


5. Einzigkett 
Durch die Anschauung iiberhaupt nicht zu erfassen ist das Problem der Ein- 
zigkeit emer Lésung der GS-DGlen. Ohne den in II.6 zu fiihrenden Einzigkeits- 
beweis ware es allein unter Berufung auf die Anschauung niemals méglich zu 
entscheiden, ob zu einem vorgegebenen Anfangsprofil nicht vielleicht doch zwei 
oder mehr verschiedene GS-Profile gehéren kénnten, von denen physikalisch 
jeweils natiirlich nur eines verwirklicht werden kann. 


6. Differenzierbarkeitseigenschaften 

Der Nachweis der vorstehend in den Nummern 1 bis 5 angedeuteten Eigen- 
schaften von GS-Profilen wird in Teil II mit Hilfe von Abschatzungssatzen fiir 
die GS-DGlen erfolgen. Als Beispiel dafiir, daB in solchen Abschatzungssatzen 
die getroffenen mathematischen Voraussetzungen eine entscheidende Rolle spielen 
kénnen, médgen zwei GS-Profile betrachtet werden, die sich nur ,,wenig“ von- 
einander unterscheiden. Man wird dann i.a. erwarten, daB auch alle daraus 
hervorgehenden GS-Profile stromabwarts nur »gering’ differieren werden. Es 
mu8 jedoch darauf hingewiesen werden, daB diese physikalisch scheinbar wohl- 
begriindete Erwartung wesentlich von den mathematischen Voraussetzungen 
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abhangt, die man iiber die betrachteten GS-Profile trifft. Besitzen namlich die 
GS-Profile die in Abb. 7 skizzierte Gestalt, wo das Profil (4) stetig partiell nach y 
differenzierbar sein soll, wahrend das Polygon-Profil (¢) nur stiickweise differen- 
zierbar ist, so wird die genannte Erwartung nicht zutreffen. Vielmehr wird sich 
bei konstanter AuBengeschwindigkeit das Profil («) stromabwarts etwa in die 
skizzierte Gestalt (#) (Abb. 8) verandern, wahrend man durch Einsetzen in die 
GS-DGlen sofort feststellt, daB das fiir alle x-Werte unveradnderte Profil (@) =(“) 
die GS-DGlen befriedigt, falls die Wand undurchlissig gedacht ist. Da das fiir 
jedes Polygon (*) gilt, kann man also in Abb. 7 das Profil (@) durch das Profil (®) 
beliebig genau approximieren, ohne da8 sich dann in Bild 8 auch (#) und () 
einander nahern werden. Unter der bloBen Voraussetzung stiickweise nach y 


Y(X0,y/ U(X, y) 
Abb. 7. Erklarung im Text Abb. 8. Erklarung im Text 


differenzierbarer GS-Profile u(x, y) kann also solch ein Abschatzungssatz i.a. 
nicht richtig sein, man wird vielmehr mindestens partielle Differenzierbarkeit 
von u(x, y) nach y fiir alle betrachteten y-Werte vorschreiben miissen. 


I. Abschatzungs- und Einzigkeitssatze bei einer Klasse 
von partiellen expliziten parabolischen Differentialgleichungen 


Samtliche im folgenden auftretenden GréBen werden als reell, alle Funktionen 
als eindeutig erklart vorausgesetzt. 


1. Problemstellung bet beschranktem Grundbereich 


Es sei & eine beschrankte, offene Punktmenge des (m- 1)-dimensionalen 
euklidischen (x, )-Raumes der Veranderlichen x und p=(jy,..-,¥»). Das 
kleinste offene Intervall, das die Projektion von ® auf die x-Achse enthalt, sei 
die Strecke x»<x<x,. Fir die Randpunkte von & werden folgende Bezeich- 
nungen eingefiihrt: ©, bezeichne alle Randpunkte von & mit x=% 9. G, sei die 
Menge derjenigen Randpunkte (x, p), zu denen es eine Halbumgebung 0S — x <0, 
|p —y|<62(d>0) gibt, die ganz in © liegt. ©, enthalte die restlichen Randpunkte 
von ® (vgl. Abb.9 und Nacumo & Simopa [14]). Bezeichnet man die abge- 
schlossene Hiille von ® mit 8 und setzt C=C,+€,, so gilt B=G+C+C,. 
sete (iVietes un. Upo. 8-91 Van On on nl ad (4 DD, ww) definiert fiir 
(x, »)€G+G, und beliebige Werte uw, v= (v%, ..., Un), W=(W,---, Mm). Ist die 


m 
2 Es sei wie tiblich | y|?= > yi. 
1 


u= 
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Funktion z(x,y) auf 8 erklart und stetig, existieren in © +(, die partiellen Ab- 


Oz OFZ 


wd’ % 49,= (w=1,...,m) und die hintere partielle Ab- 


leitungen 2y,,= 


; oy Oy2 
leitung te i 
2 PE ie 2(%,)) —4 (~¥—h,) 
[4 04 n—>+0 h ; 


so soll z eine ,,zulassige Funktion“ heiBen. Die Gesamtheit der zulassigen Funk- 
tionen bilde die Menge 3. Zur Abkiirzung werde noch Ba ys ceo ee 


Bey nee nd Cosco ewen. 
Es sollen Aussagen gemacht werden itber Lésungen z(x, y)€3 des Rand- 
a R By = T(x, ), a; 4y) 2y ») in & Fi Gy, (1 Ay) 
“| z(%,) =U(x,y) auf © (422) 


wobei die Randwerte U(x, ) auf © gegeben sein mégen. Dabei werden im fol- 
genden auf der rechten Seite der DGI (1.1) nur solche Funktionen / zugelassen, 


- die in ihrem Definitionsbereich in den m 
aay Y letzten Verinderlichen 1 schwach monoton 


ie) a anwachsen, fiir die also 
e G (— (x, 9, 4, v, ) = f(x, 9, 4, d, 1) 
\ fe) Wee gilt, falls w,2w, (w=1,...,m) ist. Die Ge- 
nt isle samtheit dieser Funktionen / mége mit Z 
% ii —* hezeichnet werden. 
pbion Das Grundechict G und sein Rand Zur Abschatzung der Lésungen des Rand- 


wertproblems R werde noch die Klasse © 
der ,,Oberfunktionen“ (bzw. Ul der ,,Unterfunktionen“) eingefiihrt, die alle Funk- 
tionen z€ 8 umfassen sollen, die den Ungleichungen 


Sey 1\%s 0:2, oa ey) le Or Oye (1.3) 
AS) etl ae 5) a ance ee (1.4) 
genugen. 
2. Nichtbeschrankter Grundbereich 


Ist & beziiglich y nicht beschrankt, dann sind noch ,,unendlich ferne Punkte“ 
einzufiihren, auf denen Randwerte U(x, )) vorzuschreiben sind, weiter ist die 
Stetigkeit von z(x,y) auf dem ,,durch unendlich ferne Punkte abgeschlossenen‘‘ 
Bereich ¥ zu erklaren. Das soll in folgender Weise geschehen: Es sei § ein offenes 


Gebiet des euklidischen (x, ))-Raumes, das die abgeschlossene Hiille G von @ 
enthalt. Durch die Abbildung 


Ws a aaa et Sa) (425) 
werde §) topologisch auf ein beschranktes Gebiet * des euklidischen (x*, y*)- 
Raumes abgebildet. Diese Abbildung fithrt also ® in ein beschranktes Gebict 
(§* und die Randpunkte von © in Randpunkte von @* iiber. 

Ein Beispiel fiir solch eine Abbildung Y ist die folgende: & sei nicht beschrankt 
beztighich der Koordinaten y;,(A=1,..., 1). Man setzt t=; Ve) und habe 
die restlichen y-Koordinaten zu dem Vektor 3 zusammen, 50 daB Sis (= 
(x, x, 3) schreiben la8t. Dabei sei etwa |3|<K fiir (x, »)€B. Das Gebiet'S sel 
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durch die Ungleichungen %)»—1<%<%,+41, |r|<~, 
Wahlt man nun 


<K-+1 abgegrenzt. 


cae’ 


WU: Ger eas AG jae 


fan yest. 
Cees: (1.6 


so ist dadurch eine eineindeutige und in beiden Richtungen stetige Abbildung 
von §) in ein beschranktes Gebiet * gegeben, bei der die unendlich fernen Punkte 
des y-Raumes auf die Oberflache der Einheitskugel im y*-Raum abgebildet werden. 

Man erklart nun die Rander Gy, CF, Gf und €*=CF +C* genau so wie in 
Nummer 1 und setzt 6*=6*+€*+ CF. %$* enthalt neben den Bildpunkten von 
®& noch weitere Punkte, die sich als Bilder der ,,unendlich fernen Punkte“ von & 
deuten lassen. Die durch diese Sprechweise eindeutig festgelegten Bilder von 
_ OF, CF, CF, C* und B* mobgen mit ©, C,,,,€ und B bezeichnet werden. Man 
nennt dann auf diesem ,,durch unendlich ferne Punkte abgeschlossenen‘‘ Be- 
reich 8 eine Funktion z(x,y) stetig, wenn z(x*, y(y*))® stetig auf G* ist. Es 
mége eigens darauf hingewiesen werden, daB durch diese Festsetzungen die Be- 
griffe ,, Konvergenz einer Punktfolge auf $“, ,,Stetigkeit einer Funktion auf $8“, 
,Randwerte einer Funktion auf ¥8“ einerseits gegeniiber dem sonst iiblichen 
Gebrauch abgeandert sind und andererseits von der Abbildung % abhangen. 
Z.B. sind Randwerte U(x, ) auf © aufzufassen als Werte U(x*, y(y*)) auf C*. 
Diese Randwerte U von z sind wegen der Stetigkeit von z(«*, y()*)) auf B* 
insbesondere beschrankt. 


3. Das Nagumo-Westphalsche Lemma 


Die Grundlage fiir die folgenden Uberlegungen bildet ein Abschatzungssatz, 
der zuerst im Jahre 1939 von M. Nacumo [3] bewiesen wurde. Da die Verdéffent- 
lichung in japanischer Sprache erfolgte, scheint sie auBerhalb Japans nicht weiter 
beachtet worden zu sein, bis dieser Satz 1949 von H. WestTpHat [4] fiir m=1 
und einen rechteckigen Bereich 8 unabhangig von NaGuMo noch einmal ent- 
deckt wurde. Inzwischen haben sich eine gréBere Zahl von Arbeiten mit der 
Verallgemeinerung des Nagumo-Westphalschen Satzes befaBt, es sei hingewiesen 
auf M. Nacumo & S. Stmopa [14], G. Propi [15], R. Narasimuan [16], J. SzaRs- 
KI [17], L. Corratz [18], W. Mvax [19]. Fiir die hier mitbetrachteten unbe- 
schrankten Grundbereiche % scheint der Satz bisher noch nicht bewiesen worden 
zu sein, ebenfalls scheinen die anschlieBend gegebenen beiden Zusatze noch nicht 
bekannt zu sein. Der Einfachheit halber, und weil im spateren Teil II nur dieser 
Fall gebraucht wird, werden nur DGlen der Gestalt (1.1) betrachtet. Die Uber- 
tragung auf allgemeinere rechte Seiten /, in denen auch noch die gemischten 
Ableitungen 2,,,y,(u=-”) auftreten, sowie auf Systeme von DGlen macht keine 
besonderen Schwierigkeiten. 


Nagumo-Westphalsches Lemma. Es sei /€X. Die Funktion w€3 sei 
Lésung des Randwertproblems R, weiter gelte v€U, we O. Dann gilt auf ganz B 


v(x, 9) <u(x, 9) < w(x, 9). (1.7) 


Beweis. Es wird zuerst der linke Teil der UnGl (1.7) betrachtet. Man setzt 
d =u —2v, dabei ist nach Voraussetzung d>0 auf €. Ist die Behauptung falsch, 


3 Die Umkehrabbildung von 2% werde mit ¥= x*, )y =)(y*) bezeichnet. 
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dann gibt es Stellen (x’, y’)€G®-+ G mit d(x’, y’) SO. Ist ¥= int4s0 bestimmt 
man eine Folge von Punkten (x, ,) €G + ©, mit d(x, Da) <0 und %,—¥%. Wegen 
der Stetigkeit von d€ 8 und wegen d >0 auf © konnen die Haufungspunkte dieser 
Folge nicht auf dem Rand € liegen, mithin ist diese Folge beschrankt. Daher 
gibt es eine konvergente Teilfolge (%j,, )j,) > (*, yp) € G+, mit d(x, 9) S0. Nach 
der iiber ©, getroffenen Voraussetzung liegt noch eine ganze Halbumgebung 


osxx—x<6, |h-—b|< 3, (@>0) von (%,5) in O7E, 


durch Annaherung von kleineren x-Werten her folgt dann aus Stetigkeitsgriinden 
d(x, ) 20, also d(x, h) =0. Da ad; in 6 +G, existiert, findet man 


fer ap .. d(x, 5) —d(4 —h,9) _. a(¥ —h,¥) 
a, (¥,) =, lim, he ee ee iS ee 


h+>+0 


Ebenso folgen aus d(%, ))20 und d(x, y)=0 wegen der Existenz von 4, und 
dyy in & +, die Beziehungen dy, (%,9)=90 und dy, (%, p) 20 fiir w=1,...,m. 
Also gilt an der Stelle (%, 4) w<0r (1.8) 


und u=1, Uy, = Vy, Uy yy = vy,y, fir w=1,...,m. Daraus folgt wegen /€2 und 
vEU an der Stelle (x, p): 
es ACS 85) oe oS) =f (%,), 4, Uy, i) Sf (%,), 4, Uy, Uyy) = Uz 


in Widerspruch zu (1.8); also ist die linke Seite von (1.7) auf ganz % richtig. 
Die rechte Seite wird genau so bewiesen. 


Bemerkung. Ist fEX, vEU, wEO, dann gilt auf ganz B 
v(x, )) << w(x, y), (1.74) 


wie man sofort sieht, wenn man im Beweis w€ 8 durch w€ 9 ersetzt. Die Gl (1.7a) 
gilt also auch noch, wenn keine Lésung uw € T des Randwertproblems R bekannt ist. 


In dem Sonderfall einer gewéhnlichen expliziten DGI (1.1) (wenn also 
Ch Nd. COTS gtk aleens s ra ie 
Bila Tea O ist fiir w~=41,...,m) ist das Nagumo-Westphalsche Lemma 


sogar noch giiltig, wenn auf €, 


U(%q, 9) S u(x, ) S w(%, v) (1.9) 


zugelassen wird (vgl. etwa E. KAMKE [20] S. 82). Um die Giiltigkeit des Lemmas 
auch unter der schwacheren Bedingung (1.9) sicherzustellen, werden die Funk- 
tionenklassen ZT, 83, O und WU etwas abgeandert zu Klassen mit dem Index 1: 
die Funktionenklasse Z, bestehe aus der Gesamtheit aller Funktionen 


I, Y), u, D, w)E®, 


die sogar noch fiir (x, 5) CG+C,+C, und beliebige Werte w, v, w erklirt sind. 
Die Funktionenklasse 8, bestehe aus denjenigen Funktionen z€8, fiir die z, Ry 
und z,, noch im Durchschnitt des Streifens Xy<X<xX%)+H mit G stetig* sein 


und stetig auf ©, hinauf fortsetzbar sein mégen. Die Klasse ©, (U,) mége alle 


* Bei nichtbeschranktem Grundbereich ist die Stetigkeit im Sinne von Nummer 2 
aufzufassen. 
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Funktionen z€ 8, umfassen, die die Ungl (1.3) sogar in © +€)+(, erfiillen und 
die auf dem Rand € den Relationen 


Ae Oe = Oleeeb) eaul: Co (1.10 
2(x,) <)U(x,y) auf ©, im 


genugen. 


1. Zusatz zum Nagumo-Westphalschen Lemma. Es sei {€Z,. Die Funktion 
u€ 8, lése das Randwertproblem R, wobei die DGI (1.1) sogar in ganz © + €,+€, 
erfiillt sein mége, weiter sei v€ Uund w€ O,. Dann gilt (1.7) noch auf 6+ €,+,. 

Beweis. Genau wie oben betrachtet man die Stellen (x’, y’)© G+, mit 
d(x’, )) =u(x', h’) —v(x’, b’) SO und setzt x =infx’. Dabei ist diesmal x= x, 
zunachst nicht ausgeschlossen; es soll jedoch gezeigt werden, daB tatsadchlich 
X= % ist. Angenommen, es sei ¥ = %), dann bestimmt man wieder eine Punkt- 
folge (*,, 9.) mit d(x,, ),) <0 und x,—¥%, doch soll diesmal (x,, ),) so gewahlt 
sein, daB sogar d, (%,, ),) =0 =(0,..., 0) und dy, (x, 0,) 20 gilt fiir ~=1,..., m. 
Wegen d(x,,))>0 auf ©, ist solch eine Wahl stets méglich. Bei beschranktem 
Grundgebiet & gibt es dann eine konvergente Teilfolge (x;,, );,) >(%o, 9). Wegen 
der Stetigkeit von d, d, und d,, und wegen d(%,, D,,) SO, dy (%%,, De) =9, 
Ayn yy (%y» Dy) 20 ist dann auch d(%,H)S0, dy(%,H)=0, dy, y,(%, 9) 2 O. 
Wegen d(x , )) 20 ist sogar d(x», )) =O und damit gilt 


Vu 


d, (Xo, b) >t (%o, y, u, Uy, Uy») re f(%9,, a Uy» oe) == 0. 


Wegen der Stetigkeit von d; ist dann in dem Durchschnitt D einer Halbumgebung 
XySx<x%y+6', |p-—y|<6’ mit dem Gebiet G+, stets dZ(x,y)>0. Sei 
Xp SasxSp<x%+0”’, 9=p, ein Teilintervall, das ganz in D legt, wobei noch 
(x, 9,) €€ gelten mdge. Dann ist d(«,,) 20 und d;(x,),)>0inaSx<p. Nach 
einem bekannten Satz (vgl. etwa E. KAMKE [20], S. 92/93) ist dann d(x, ),)>0 
fiir «<*<. Da dies fiir jede derartige Strecke gilt, ist d>0 im Innern von D 
in Widerspruch zur Definition von (x,y). Ist das Grundgebiet unbeschrankt, 
dann kénnen dieselben Schliisse nach Nummer 2 in dem beschrankten Gebiet ©* 
durchgefiihrt werden. Mithin fiihrt die Annahme % = x) zu einem Widerspruch, 
also ist ¥ >, und der obige Beweis kann wortlich weitergefiihrt werden. 

DaB auch die Voraussetzung iiber die Stetigkeit von w€ 3 auf dem Rand © 
von % noch gemildert werden kann, zeigt der (im folgenden allerdings nicht 
weiter benutzte) 


2. Zusatz zum Nagumo-Westphalschen Lemma. Die Funktion u(x, ) sei 
nur stetig in © +, erfiille jedoch sonst alle anderen Voraussetzungen, die an 
die Funktionen aus der Klasse 8 gestellt werden; weiter sei ve U und wE€®. 
Fiir jede Punktfolge (x,, y,)€ G+, mit (x,, y,) (x, y) CW gelte stets v(x, 9) < 
inf u(x,, »,), sup “(%,, 9,)<w(x, 9). Dann gilt (1.7) auf ganz B. 

Der Beweis folgt fast wértlich dem oben angegebenen Beweis des Nagumo- 
Westphalschen Lemmas. 

4. Abschatzungssdaitze 

Aus dem Nagumo-Westphalschen Lemma lassen sich sofort Schranken fiir 
die Differenz zweier ungefahrer Integrale von zwei DGlen der Gestalt (1.1) 
angeben. Es gilt namlich der folgende 
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Abschatzungssatz®. Es sei fE%, gc, ue ZB und vB. Die Funktionen 
&, (x), €(#), E, (x), Eo(x), 6(*) und A(z) seien in %S%*S% stetig, die Funk- 
tionen w(x, ¢) und Q(x, t) seien fir »S*S%, |t|< co beschrankt und stetig. 


Auf 6 +©, gelte 


& (x) Suz —F (4, 0,4, My, My) SE (A), (1 
&4(x) Sv; — g(%, Y, UV, Uy» Uy) SEs (%), (1 
6(%) Sg(%, 0, 0, Uy, My) —1(%, 0, Vy, Upp) S A(X), (1 
w(%,v — 4) Sf(%, 9, 2, %, Uy) — } nD FOr 55 1) <O(x,v—u). (1 
Auf © sei 
pSv(x,9) — u(x, 0) <M. (! 


Die Funktion s(x) sei in %»< +X, definiert als Minimalintegral der DGl 
y’ == Min (0, w(x, y) — Ex (*) + €2(x) + 6(x)) (1 


unter der Anfangsbedingung 


Y (Xo) = 4; (1 
S (x) sei dort das Maximalintegral der DGl 
y’ = Max (0, Q(x, y) — e(*) + E,(x) + A(a)) (1 
unter der Anfangsbedingung 
y (%o) = M. (1 


Dann gilt auf ganz 8 die Abschatzung 


$(%) Sv (%, 9) —4(x, oy) = S(a). Ct 


Beweis. Man betrachtet die Differenzfunktion d=v—vwu. Dafiir gilt 


d, =U, — U, = 8%, 9, v, %, Dvn) a aoa Vy» Oe -- 

+ f(%, 0, v, Ops Vy y) a) (ema a or) + 

+ 1(%, 0,4, Uy, Uyy) — 1(%, 9, 6, Uy, Uy) — & (4) = Ea(x) 
<f(%, 9, 4, u%, + dy, Urniat Oyu fy Doll, Uy , Uy») —. 

+ 22(%, d) — e,(x) + E, (x) + A(x) + «, 


| 


1) 
12) 
13) 
14) 


15) 


16) 


A7) 


18) 


AQ) 


21) 


wo e>0 sein soll. Auf dem Rande © gilt nach Vor. (1.15) d<M+e. Nach 


dem Existenzsatz von PEANO hat die DGl 


y’ = Max (0, Q(x, y) — e,(%) + Ey (x) + A(x)) + .¢ 
unter der Anfangsbedingung 


¥(%) =M+ « 


mindestens eine Lésung in x»<~x<.x,. Eine solche werde fiir jedes e>0 aus- 
gewahlt und mit S(x, «) bezeichnet. Nun ist offenbar S (x, e) eine Lésung der 


® Vgl. dazu L. GruLtano [21], L. Corratz [18] und J. Szarsxi [17]. 
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DUnGl - 

e By = Pe, Y, U, U = 4y > Uyy ae aan) oe ACD Y, U, Uy» Uy») + 
+ 22(%, 2) — 4 (x) + E,(x) + A(x) + «, 

die sich aus (1.21) durch Ersetzen des <-Zeichens durch das =-Zeichen ergibt. 

Auf dem Rande € ist weiter d<M+eXS(x, ¢). Da schlieBlich noch die rechte 


Seite von (1.22) eine Funktion der Klasse & ist, gilt nach der Bemerkung zum 
Nagumo-Westphalschen Lemma auf ganz 8 die UnGl 


d(x, 9) <S(x, 8). 


Fiir e—0 strebt S(x, e€)>S(x) (man vgl. dazu etwa E. KAmKE [20], S. 83ff., 
wo zwar der Anfangswert festgehalten wird — das wiirde hier S(x,2)=M 
bedeuten —, doch laBt sich der dort angegebene Beweis sofort auf den Fall 
S(%9, 8) =M-e iibertragen). Also gilt auch noch 


VAC Pas) eo Ca (1.23) 


(1.22) 


In genau derselben Weise laBt sich auch zeigen, daB auf ganz % gilt 
Sys ala, by, (1.24) 

und die beiden UnGln (1.23) und (1.24) zusammengenommen ergeben gerade 
die Behauptung. 

Beispiele. 

a) Wahit man —o(x, t) =Q(x, t)=K(x) mit einer in %Sx<x, stetigen 
munktion K(x)" und. setzt noch —pw==M)--— é,(*%) =E, (x), — 23%) =Es{%), 
— 6(x) =A(x), so findet man aus (1.16), (1.17) und (4.20) die Abschatzung 


Jw, ) a(x, I SM+STKO +E +E) + AQ] a. (1.25) 


Diese Abschatzung wird fiir den Sonderfall konstanter Funktionen K(x), E, (x) 
und £,(x) und fiir A(x) =0 schon von H. WEsTPHAL [4] S. 692 angegeben. 
b) Wahlt man 
=O, ty = Ne, t= KK (&) b-ER (%) 


mit in (%»), x,> stetigen Funktionen A(x) und k(x), und setzt wieder —u=M, 
— e,(x) =E, (x), —&,(x) =E,(x), — d(x) =A(x), so erhalt man aus (1.16), (1.17) 
und (4.20) die Abschatzung 


|v (x, 9) — u(x, »)| 


x” x t (1.26) 
<exp| (Kai) {M+ [AO +60 +E.) +40] (exp—SK(6) ds) at] 


und in dem Sonderfall 
TS 9 res 


k(x) + E, (x) + E,(*) + A(x) = e = const. 
die vereinfachte Formel 


| (x, 0) —m (x, oy SM RO 4 5 [eB — 1], (1.27) 
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Die Abschatzung (1.27) mit M=0 wird bei H. Wesrpuat ([4], S. 693) zur Her- 
leitung eines Einzigkeitssatzes benutzt. Bei gewohnlichen expliziten DGlen 
4. Ordnung ist (1.27) bekannt, man vel. z.B. E. KAMKE [20], S.94 fiir d=0. 

Ersetzt man in (1.14) die Funktionen w(v, ¢) und Q(x, ¢) durch allgemeinere 
Funktionen w(x, y,f) und Q(x, 9, #), so lassen sich leicht weitergehende Ab- 
schatzungen angeben, bei denen jetzt statt der gewOhnlichen DGlen (14.16) und 
(1.18) partielle explizite parabolische DGlen auftreten. Zwei einfache Aussagen, 
die sich auf diese Weise herleiten lieBen, sollen in der folgenden Nummer 5 for- 
muliert, jedoch direkt bewiesen werden. Ebenfalls soll auf eine weitere Verall- 
gemeinerung verzichtet werden, die man dadurch erhalt, daB man statt des 
_Abstandes‘’ d=v—w zweier zuldssiger Funktionen we 8 und v€ 8 eine verall- 
gemeinerte Abstandsdefinition d (wu, v) mit d=0 nur fiir w(x, h) =v(x, y) einfiihrt. 
Ein allgemeiner Einzigkeitssatz, der sich auf diese Weise gewinnen lieBe, wird in 
Nummer 6 direkt hergeleitet werden. 


5. Randwertsdtze 


Unter gewissen Voraussetzungen gilt auch bei DGlen der Gestalt (1.1) das 
bekannte Maximum-Minimum-Prinzip, d.h. die Aussage, daB eine Lésung des 
Randwertproblems R ihre Extremwerte nicht allein im Innern von 8 annehmen 
kann. Dieses Prinzip ist als Sonderfall enthalten im folgenden Satz: 


Erster -Randwertsatz. Es sei f€2. g(x) und G(x) seien im ajoxe 
erklart, stetig und in x»<*< x, differenzierbar. Fiir (x, y) € G+, und beliebige 
Werte v gelte 


g(x) Sf (x, , v, 0, 0) S G(x), (4.28) 


wobei unter p= (0, ..., 0) der m-dimensionale Nullvektor verstanden sei. Seien 
ferner « und M so bestimmt, daB auf © die UnGl 


+ 8(%) S U(x, ») SM + G(x) (1.29) 


erfullt wird®. 


Dann gilt fiir jede Lésung w€ 8 des Randwertproblems R auf ganz 8: 
+> &(%) Su(x,) SM+ G(x). (1.30) 


Beweis. Man betrachtet fiir e>0 die Funktion v(x; 2) = M+ G(x) + 
é(1tx—x). Esist vC 8, weiter gilt Up = Vp y= 0, Vz = G'(x) +e>H(x, 0, U, Up, Vpn): 
Wegen (1.29) ist noch v>U auf © und damit v€D. Nach dem Nesumoawese 
phalschen Lemma gilt dann auch u(x, y)<v(x; e) auf ganz 8 und fiir e>0 
ergibt dies die rechte Seite von (1.30). Ganz analog folgt auch die linke Seite 
womit (1.30) vollstandig bewiesen ist. 


Fiir g(x) =G(x) =0 erhalt man daraus das 


Maximum-Minimum-Prinzip (vgl. auch L. Corzatz [18]). Es sei ii 
; : » IWS Gel pa £ 
(x, ») €G +, und beliebige Werte v gelte (78) fe be 


1 (% 9, v, 0,0) =0. (1.28 a) 


° Wegen der Beschranktheit von U auf © ist das stets mdglich. 
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Ist ferner auf € noch die UnGl 
MS U(x,5) SM (1.29a) 


erfiillt, dann gilt fiir jede Lésung w€ 8 von R auf ganz &: 
MSu(x,y) <M. (1.30a) 


Zweiter Randwertsatz. Es sei f€Z. Die Funktionen g(x), G(x), k(x) und 
K,,(%) (w=1,..., m) seien in x» x<x, erklart, stetig und in x<x*< x, diffe- 
renzierbar. Fiir (x, ))€@®-+, und beliebige Werte v und » gelte 


eoie <f(x,,v,0,0)<G ic Naas (1.31) 


= f=1 


Seien ferner « und M so bestimmt, da8 auf © die UnGl 


p+elx ) +O hal) U(x, 9) SM+ G(x ) DK x) Vu (1.32) 
Lt _ 
erfiillt wird’. ve 


Dann gilt fiir jede Lésung w€ 3 des Randwertproblems R auf ganz 8: 


m m 


e+ g(x) +d ky (*) y,S u(x, 9) SM+ G(x Nee als Lae (1.33) 


w=1 = 


Beweis. Der Beweis verlauft mit der Funktion 


v(x, );e) =M+G(« ) DK, lx *) Vu + E(1 + x — Xp) 


=1 


ganz entsprechend zum Beweis des Ersten Randwertsatzes. 


6. Einzigkettssatze 
Die Grundlage zu den Einzigkeitsaussagen dieser Nummer bildet folgendes 


Einzigkeitsprinzip®. Es sei f€&. Zu jeder Lésung w€ 8 des Randwert- 
problems R gebe es eine Folge von ,,approximierenden™ Ober- oder Unter-Funk- 
tionen v,€EC O(U) (k=1, 2,...) mit y%(x, 9) u(x, 9) in G+, fiir kR> ow, 

Dann hat das Randwertproblem & héchstens evme Lésung we 3B. 


Beweis. Es werde nur der Fall der Oberfunktionen betrachtet, der Beweis 
fiir Unterfunktionen geht analog. Seien w€ 8 und w€ 3 Lésungen von R und 
an der Stelle (x,,9,)€ G+, sei w(x, ,) >u(%,,9,). Nach Vor. gibt es dann 
eine Folge v,€ © mit v,—u, also gibt es eine Zahl K derart, daB v, (x1, 9;) <u (x1, 94) 
ist fiir alle R>K. Nach dem Nagumo-Westphalschen Lemma gilt aber fiir alle 
v,(k=1, 2,...) auf ganz B stets v,>%, also insbesondere v,(%,, ))) > 4 (%, 04) 
fiir k>K, womit der gewiinschte Widerspruch erreicht ist. 


? Das ist fiir beschrankte Bereiche 8 stets méglich. 

8 Man vel. dazu bei gewohnlichen expliziten DGlen erster Ordnung die notwendige 
und hinreichende Bedingung fiir die Einzigkeit einer Lésung von K. Kunua1, die 
von T. YosIE [22] mitgeteilt wurde. 
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Zusatz. Ist f€3,, wE€ 8B, Lésung von R und v,€ ©, (U,), so gilt dieselbe Be- 
hauptung. Beweis wie oben mit dem 1. Zusatz zum Nagumo-Westphalschen 


Lemma. 
Als Beispiel zu diesem Zusatz folge die Ubertragung des Satzes von NAGUMO 


(vgl. z.B. G. Sansone [23], S. 102ff.), von gewohnlichen DGlen auf das Rand- 
wertproblem R: 
Einzigkeitssatz nach NAGUMO. Es sei /¢ Z,, weiter gelte 


ES Y, v, D, 1p) — f(x, Y, U, v, ) << fiir (x, yn E€G+C6,4+ G, (1.34) 
deca) 
—co<u<v<oo und beliebige Vektoren » und w. 


Dann hat das Randwertproblem R héchstens eime Lésung #€ 31. 


Beweis. Man setzt v4, =u+ ee fiir kl) 2h Eesist Saar 
Ov G2v ‘ - 
®+6€,+6, ist stets v,>u, weiter gilt dort fn a ae =, und schlieBlich 
noch 
o-v = 1 Vp —U 
a = Ug Ses L(%, 9, U, Uy, Uyy) + re Xo > 1X0, Ups My» Mn) 


= dvp OP Up 

=i(x, Y, UR» ye 2 
Also ist v,€ 0, und wegen v,—w fiir k—>oo wird uw von den v, approximiert. 
Damit sind alle Voraussetzungen des Einzigkeitsprinzips erfiillt, also gilt die 
Behauptung. 


Aus dem Einzigkeitsprinzip ergibt sich folgender 

Allgemeiner Einzigkeitssatz®. Die Funktionen «(x,), «(x,»), B(%, »), 
(x, h) seien auf B erklart und in © + G, gelte «(x, y) Sa (x, y)<z(x, 9) <B(x, »)S 
(x, 9) fiir jede Lésung z€ 8 von R. Die Funktion 0(x, y, u,v, w) sei stetig fiir 
(GPS, ox, )) SB (4, Y), ae, )) <v<B (x, yp), —co<w<-+ oo und sei fiir 
( 


Z 
p 
x, ny EG1C,, axu<p, i<v<fp, —0o<w<-+ co zweimal stetig partiell diffe- 
renzierbar. Weiter gelte dort noch o(x,, u, u,0)=0, e,>0 und e,<0. Die 

Gleichung 
0(%, 9, u,v, w) =0 (1.35) 


sei fiir (%, 9) €B, a<u<f, —co<w<-+ oo (und wegen 9,>0 dann auch ein- 
deutig) nach v auflésbar, ihre Lésung werde mit v =o(%, ), u, w) bezeichnet. 
Die Funktion w(x, y, u,v, tv) sei fiir (x, y) €B und beliebige Werte w, v, w 
erklart. Die Folge «2,0 (k=4,2,...) und die auf © erklarte Funktionenfolge 
Bx (%, 9) S,0 (k=1, 2, ...) seien beide Nullfolgen. Die explizite partielle paraboli- 
sche DGl 
Red (FY); & 253 Spy) or OG (1.36) 


habe unter den Randbedingungen 
2(x,)) = B,(*,)) auf € (1.37) 
° Bei gewohnlichen expliziten DGlen erster Ordnung entspricht dieser Satz un- 


gefahr einem Einzigkeitskriterium von R. Sard [24] (S.4, Théoréme 3) 
G. ZwiRNER [25], der die Voraussetzungen noch etwas allgemeiner faBt. 


, vgl. auch 
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fir k=1,2,... auf B stets eine Losung e(x, );k)€8 mit e(x, y; 2) (<0) and 
es gelte dim, é(x, 0; k) =x(x, y) =0. Es sei f€T und fiir jede Funktion u€ 3 


und alle hinreichend groBen Werte & gelte auf 6+: 


0, + 0, f(x, o, w, Uy, Uy) + Oy f(%, 9, 6, Oy, Spy) + Oy @(%, Y, ura bee One (loys) 
Dabei sind in @,,@,,@, und g, als Argumente (x, b, u,o(x, y, u, €), €) einzu- 
tragen, die Argumente von é¢ sind (x, y; A). 
Dann hat das Randwertproblem R hdéchstens eine Lésung u€ 8. 
Beweis. Es sei w€ 3 eine Lésung von R. Man setzt 
v(x,9; k) =o(x, 9, u, e(%,9;%)) mit veg. 
Wegen 
>u fir w>0 
ee Sages ist o(%,),u,w) =x fir w=—0 
<u fir w<0 


und damit strebt v->w fiir k— oo; dabei gilt auf ganz %, also auch insbesondere 
auf ©, noch v (2,4. SchlieBlich ist noch in © +, 


< 


Gs 
0= FZ, 0(%, 0, u(x, y), v(x, 9; 2), e(x, 9; 2)) 


bes eae Taal 
— 07, Oy Ux r Qy Uy FT Ow Ex 


= 02 T Ou t(%, , 4, Uy » Uy») ap Qv Vy or Ow [o (x, ), &, Ey» Ean) ai Ot. | , 


oder nach Voraussetzung fiir hinreichend groBe k 


1 
vy Oe zi Lox = Ont (4% Y, u, Uy, Uy y) a Ow (w - 0t,) | a f(x, ), v, Yy> Vy») F 


Also ist v(x, »;) fiir hinreichend groBe Werte von k eine Folge von Ober- 
(Unter-) Funktionen mit vw fiir k-~oo. Nach dem Einzigkeitsprinzip gilt 
dann die Behauptung. 

Beispiel. Man setzt o(%,),u,v,w)=v—u—w. Damit ist o,=1>0, 
Ow = —1<0, o(%, 9, u, wv) =u+w. Fir w(x, 9, 4, 0, w) =a (x, wu) wird (1.36) zu 
der gewohnlichen DGIl y’=q@(x, vy) +a, bei der jetzt natiirlich statt Rand- 
werten auf © nur noch Randwerte auf ©, vorzuschreiben sind. Es sei w(x, f) 
stetig und beschrankt fiir 7»<*< 4, und beliebige ¢ sowie w(x, 0)=0; die DGl 
y’=a(x, y) habe unter den Anfangsbedingungen y(%») =0 in ¢%9, *,> nur die 
triviale Lésung y (x) =0. Ferner seien die beiden Nullfolgen «, und f, (k =1, 2, ...) 
mit «,>0 und f,>0 vorgegeben. Nach dem Existenzsatz von PEANO hat die 
DGl y’=o(x, y) +a, unter der Anfangsbedingung y(x))=£, fir alle & in 
%p << xXx, mindestens eine Lésung. Fiir jedes k werde eine dieser Lésungen aus- 
gewahlt und mit e(x; k) bezeichnet. Nach einem Abschatzungssatz fiir gewohn- 
liche DGlen, der analog zum Nagumo-Westphalschen Lemma lautet (vgl. 


10 Diese Voraussetzungen sind z.B. erfillt, wenn w¢Z ist, wenn @(¥%, b, u, 0, =O 


gilt und «= - p= 7 gesetzt wird, weil dann offenbar e(%, y; k) = 
pra —*o >0 eine Lésung von (1.36) unter den Randbedingungen (1.37) ist. 
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KamgE [20], S. 82), ist dann auch e(x; 2) >0 auf %»<*<x,. Nach einem weiteren 
bekannten Abschatzungssatz gilt auch e(x; k) >y(*) =0 fiir kR—> oo (vgl. z.B. 
Kame [20], S. 87; dort ist — abweichend von den hier getroffenen Annahmen — 
w(x, t) zwar in einem offenen Gebiet als stetig vorausgesetzt, doch laBt sich w 
offenbar stetig auf die ganze Ebene fortsetzen, so daB auch diese Bedingung 
befriedigt werden kann). Mit ¢,=é,,=9 wird dann die Forderung (1.38) zu 


1 (x,0, 6+ €(%5R), My» My) — F(%, 0s Ms My My y) SO (se(ash)ye 


Diese Bedingung ist unter den Voraussetzungen des folgenden Satzes von 
L. Grut1Ano [21] (siehe auch J. Szarsxi [17]) sicher erfiillt, so daB sich der 
Einzigkeitssatz von GIULIANO als Sonderfall des obigen allgemeinen Einzigkeits- 
satzes erweist: 

Einzigkeitssatz von GIULIANO. Es sei fC2. Die Funktion w(x, ¢) mit 
w(x, 0)=0 sei stetig und beschrankt fiir x» *<%, und beliebige ¢. Die DGl 


y= w(x, 9) (1.39) 
habe unter der Anfangsbedingung 
y(%) = 0 (1.40) 


in %SxX%, nur die triviale Losung y(x) =0. Weiter gelte noch 


(x, 9, v, v, ) — f(x, y, uw, v, W) S w(x, v — 4) (1.41) 


fiir (x, 9) €G+C,, — o<u<v<-+ oo und beliebige Vektoren » und 1. 
Dann hat das Randwertproblem R héchstens eine Lésung u€ 3. 


Bemerkung. Der Einzigkeitssatz von GIULIANO 1aBt sich auch sofort aus dem 
Abschatzungssatz von Nummer 4 ablesen, wenn dort /=g, ¢&,=e,=£,=E, 
=A 1 — M=0 gesetzt wird. 


II. Einige allgemeine Eigenschaften von zweidimensionalen laminaren und 
stationdaren Wand-Grenzschichten in inkompressiblen Medien 


In diesem Teil II werden einige Eigenschaften von Lésungen der GS-DGlen 
hergeleitet, die sich hauptsachlich auf ihr geometrisches Bild, das ,,Grenzschicht- 
profil’, beziehen. Diese Aussagen werden dabei nicht aus einer expliziten Dar- 
stellung der allgemeinsten Lésung der GS-DGlen gewonnen, auf die Frage der 
Existenz einer solchen Lésung wird sogar iiberhaupt nicht eingegangen. Es 
gelingt vielmehr, durch Anwendung der allgemeinen Ergebnisse von Teil I auf 
die speziellen GS-DGlen alle Aussagen direkt aus den definierenden Gleichungen 
abzuleiten. ‘ 

Die dabei getroffenen Voraussetzungen beziiglich Stetigkeit und Differenzier- 
barkeit der auftretenden Funktionen sind gering und diirften wohl in allen 
praktisch vorkommenden Fallen erfiillt sein. Dagegen muB auf drei physikalische 
Einschrénkungen hingewiesen werden, die fiir das Folgende wesentlich sind: 
Finmal wird die AuBenstrémung als drehungsfrei vorausgesetzt, zum andeved 
wird nur der Fall der ruhenden Wand betrachtet und schlieBlich sind zwa 
beliebige Absauge- und Ausblasegesetze zugelassen!, jedoch soll dabei oe 


Scheu 3 . 
1 Bis auf die Nummern 4 und 5, in denen nur Absau 


rae en zugel 
undurchlassige Wand vorausgesetzt wird. S gelassen bzw. sogar 
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senkrecht zur Wand abgesaugt bzw. ausgeblasen werden. Wird eine von diesen 
drei Annahmen verletzt, dann sind die formulierten Ergebnisse i.a. nicht mehr 
richtig. 

Es mége schlieBlich ausdriicklich betont werden, daB die folgenden Aussagen 
liber eine Grenzschicht in einer stroémenden Flissigkeit fiir die physikalische An- 
schauung alle mehr oder weniger ,,selbstverstandlich sind. Um diesen Zusam- 
menhang mit der Anschauung noch starker hervortreten zu lassen, werden die 
nachfolgenden Ergebnisse (Nummer 2 bis 6) in Nummer 7 in einer physikalisch- 
anschaulichen Sprechweise noch einmal zusammengefaBt und durch Skizzen 
erlautert. 


1. Die Grenzschicht-Differentialgleichungen 
a) Die Prandtlschen GS-DGlen. Es sei U(x)>0 in x, SxS, erklart, 
stetig und stiickweise stetig differenzierbar; vg (x) sei dort erklart und stiickweise 
= als 
jl—S 


Funktion von s in 0XsS1 stetig differenzierbar, weiter gelte lim a (y) = U(x»). 
y— oo 


stetig. Die Funktion %(y) 20 sei in 0Sy< oo erklart und es sei ii( 


Man betrachtet die DGlen der zweidimensionalen laminaren und stationdren 
GS in einem inkompressiblen Medium der kinematischen Zahigkeit y (vgl. z.B. 
H. ScHLICHTING [2], S. 93): 


UU, + UU =UU'+rn,, (2.1) 


Ut =0. (222) 


Dabei sei »>0, w=u(x, y),-v=v(s, vy); 2 und y mdgen auf dem Bereich 
Bs {%)SxS%,, OS y< co} variieren und u,v médgen noch den folgenden An- 
fangs- und Randbedingungen geniigen: 


u(%,¥) =%(y), (2.3) 
u(%, 0) =0, (2.4) 
u(x, 00) = U(x) ®, (2.5) 
v(x, 0) = v9(x). (2.6) 


Das Randwertproblem (2.1) bis (2.6) sei wieder mit R bezeichnet. Anstatt ein 
Funktionenpaar u(x, y), v(x, y) als Lésung von R zu betrachten, soll im folgenden 
v(x, y) vermége der Glen (2.2) und (2.5) als Abkiirzung fir den Ausdruck 


CN eee (2.7) 


stehen und damit allein u(x, y) als ,,Lésung“ von (2.1) unter den Bedingungen 
(2.3) bis (2.5) aufgefaBt werden. 


12 Diese Schreibweise soll lim u(x, y) = U(x) bedeuten und wird durch (2.5a) noch 
=> CO ; 
naher prazisiert werden. Es laBt sich zeigen (W. RHEINBOLDT [26], K. NicKEL [27]), 
daB diese ,,auBere Randbedingung“‘ unter gewissen Voraussetzungen stets von selbst 
erfiillt ist und daher weggelassen werden kann. 
Ach.. Rational Mech. Anal., Vol. 2 2 
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Der Grundbereich % fiir die unabhangigen Verdnderlichen x und y ist nicht 
beschrankt. Man geht daher gemaB I.2 etwa durch die eineindeutige Trans- 


formation 
ih 


7 BM, 
ates mit der Umkehrung ath FG (2.8) 
reese Soa och 

iiber auf den beschrankten Bereich 8*: {x»S7S%, OSs}. 


Damit werde eine Funktion u(x, y) dann als Lésung des Randwertproblems R 


bezeichnet, wenn gilt: w(x, y) sei auf % erklart und =O,“ (r, <5) sei als Funk- 
tion von rv und s sogar auf §* erklart und stetig. Die partiellen Ableitungen 
Uy, Uyy und w,y, seien als Funktionen von x und y in &©+ 6,” erklart, u, [, <] 


sei auf $* noch erklart und stetig!4. Weiter gebe es im Intervall <9, x,> eine 
Zerlegung %)=Yo<Vi<0+'<Yn = % derart, daB uw, und w,, auf jedem der Streifen 
Vj;S4XSYj41, OS y< 00 existieren4 und wu, dort stetig ist. w(x, y) gentige zu- 
sammen mit der durch (2.7) definierten Funktion v(x, y) auf 8 der DGl (2.1) 
und erfiille die Anfangs- und Randbedingungen (2.3) bis (2.5), wobei (2.5) als 


se eer (7, aa) = U(x) (2.5 a) 
aufgefaBt werden soll. Die Gesamtheit dieser Lésungsfunktionen des Randwert- 
problems R sei mit & bezeichnet. Zwar ist ein groBer Teil der nachfolgend 
formulierten Ergebnisse schon unter schwacheren Voraussetzungen giiltig, doch 
sollen der Einheitlichkeit halber im folgenden nur noch Funktionen u€& be- 
trachtet werden. 


b) Die v. Misessche Transformation. Nach v. MisEs [28] mégen in (2.1) 
bis (2.6) neue Koordinaten eingefiihrt werden durch: 


SY tana 
(x, ¥) mo u(x, t) dt + n(x) (2.9) 
mit 
no(x) = — frvg(t) at. (2.10) 


Dadurch geht % tiber in den (unbeschrankten!) Bereich 
Be {x SES Mm, 9(E) S7< co} 


der €, 7-Ebene. Da die Funktionaldeterminante Se ti Ast. erweletiercite 
v. Misessche Transf “ Oy) ; coe 
atari: ranstormation (2.9) fiir w>0 als im kleinen eindeutig umkehrbar; 


si pene hen epee ©, ©, ... seien wie in I.1 und I.2 eingefiihrt 
abei kommt auf den Randern von 8 nur einseitige Differenzierbarkeit in Frage 


‘i sd : San) 
Wegen der Stetigkeit von u (, Bee auf 8* und (2.5a) ist offenbar 


lim (4%, \) Ook 
y> oO 
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wegen der speziellen Gestalt von (2.9) gilt dies Ergebnis sogar noch im groBen. 
Die danach fiir w=--0 auf % eindeutigen Umkehrfunktionen zu &, ™ mogen mit 
x (€), y (§, 7) bezeichnet werden. 
Setzt man nun 
g (En) = U? (x (8) — #(% (6, 7), (2.41) 
so wird 
UE 


eee (2.12) 


f= 2|UW uu, vn, |, g, == 24,,  8,,= — 2 7 


und (2.1) geht damit iiber in die partielle explizite parabolische DGl. 


=v |) OF e,.- (2.13) 
Weiter sind noch die folgenden Anfangs- und Randbedingungen erfiillt: 


8 (%9,) = U? (xq) — a? (¥ (xo, n)), 
g (&, m0 (€)) = U7 (6), (2.14) 
eg ele i (", <=) pur 


Wegen wE& ist (EW). n(r, <=3)} auf §* stetig und auf 6+: {ig by, 
No (€) << oo} existieren g,, die hintere partielle Ableitung gz und — soweit 
u>0O ist — auch g,,. Also ist. g(&, 7) fir wc 8 und w>0 in & + ©, im Sinne von 
Abschnitt I eine Lésung fiir das Randwertproblem (2.13), (2.14). Da die rechte 
Seite von (2.13) wegen v>0 zur Klasse Z gehort?®, lassen sich fiir w-- 0 in 6+ G 
die Satze des Teil I auf das Randwertproblem (2.13), (2.14) anwenden. 


c) Die CRocco-Transformation. Nach L. Crocco [29] kann man in (2.1) 
bis (2.6) auch auf andere Weise neue Koordinaten einfiihren, namlich durch 
Q(x) =* 
a(x, y) = u(%, 9). 


% geht durch (2.15) iiber in einen beschrankten(!) Bereich der @, o-Ebene, z.B. 
fir uy=0 in 8: {4<e<%, 0So<U(o)}. Die Funktionaldeterminante ist 


(2715) 


ae Z =u,, und fiir u,+-0 ist (2.15) sogar im GroBen eindeutig umkehrbar. 
gay ° “A A 
Die danach fiir u,-+-0 eindeutigen Umkehrfunktionen mégen mit %(@), ¥(@, 0) 


bezeichnet werden. 


Setzt man nun 


7 (9,0) = uy (%# (9), 9(@,9)), (2.16) 
so wird 
Ux U. U uU,,U —uU 

T) = Uzy — ae mie ties v7, a Vy (2:17) 


16 Die rechte Seite von (2.13) ist zunadchst nur fiir g =U? als reelle Funktion 
definiert, doch 1aBt sie sich durch 
) iybbe oe (Ue 
(ens tred={,yormee, fir 4S0" 


fiir (€, 7)€G +, und beliebige Werte g, g,, &,, erklaren. 
2% 
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und (2.1) geht damit iiber in die partielle explizite parabolische DGl 


7 (GRU 
gt es Ve erg stn (2.18) 


An Anfangs- und Randbedingungen gelten 


1 (0, 0) = u, (%(@), 0) ”” (2.19) 
t (0, U(e)) =0. 
Wegen wE& ist (0 (7), (7, a) auf $* stetig. Da fiir u,>0 in & durch 


Ss 


o(7), ofr, ‘ die beschrankten Bereiche 8* und %$ topologisch aufeinander 
1S) 


abgebildet werden, ist damit t(o, 0) auch stetig auf ®. Weiter existieren fiir 
u,>0 auf 64 6: {%<eSm, 0<o<U/(0)} die partiellen Ableitungen 1,, Too 
und die hintere partielle Ableitung 7,, d.h. fiir w€2 und w,>0 ist t(Q, 0) im 
Sinne von Abschnitt I eine Lésung fiir das Randwertproblem (2.18), (2.19). Da 
noch die rechte Seite von (2.18) wegen v>0 eine Funktion €@ ist, lassen sich 
die Satze des Teil I auf das Randwertproblem (2.18), (2.19) anwenden, sofern 
u,== 0 bleibt. 


2. Ubergeschwindigkeiten 


Definition. Man setzt e= sup [u(y)]—U(%)). Es ist stets eZ0. 
0sy<0~ 


Satz. Sei uC, dann gilt aut % 


u(x, Vy) SVU? (x) + 2e U(x) + e2 38. (2.20) 


Beweis. Es sei zunachst u(x, y)>0 in © angenommen. Man betrachtet die 
Funktion g(é,7) nach (2.11) auf dem durch die vy. Mises-Transformation (2.9) 
erzeugten Bereich 8. Da die rechte Seite von (2.13) fiir g,,, =0 identisch ver- 
schwindet, da ferner g(é, 7) auf % stetig ist und gz, g, und £n, dort existieren, 
gilt auf 8 das Maximum-Minimum-Prinzip von I.5. Damit gilt auf ganz % 
die UnGi 


g(é,m) 2 inf g(6,m) = inf g(E,m) = — e(e + 2U(m)). (2.24) 


Gabe es nun eine Stelle (x’, y')EB8 mit u(x’, y’)> VU? (x’) +2 U(%) + e?, dann 
gilte an der Stelle f= 8 (x'), 1 = (a ay) -aaen g (E', n')<—el(e+ 2U(%)) in 
Widerspruch zu (2.21). Also ist (2.20) auf ganz & erfiillt. 

. Gibt es Stellen aus & mit ~=0, dann werde etwas anders geschlossen. Es 
sei G = ink eie.7) =O. Angenommen, alle Punkte (&’, 7’) mit g(&’, 7’) =G mogen 


in G+ ©, liegen. Man bestimmt & = inf £’ und sucht eine Stelle 7n mit g(E,7) =G. 


1? Dieser Randwert ist i.a. nicht vorgeschrieben, zur Aufl 1 
; | _W ; osung von (2.18) wird 
es AG die Randbedingung 1(@, 0) [vt (@, 0) — vo (@)] ae ee 
ie sich aus (2.1) fiir y+0 ergibt. Fiir die foleenden A ii j ie 
eight Sarees g genden Aussagen gentigen jedoch die 
#8 In einer Arbeit von H. WirttinG [30] wird fiir kleine Stéru 1 i 
; | torungen eines GS-Profil 
ein lokales Anwachsen gezeigt. In Erganzung dazu wird dealt (2.20) fiir oe 
liebiges Anfangsprofil %(y) eine globale Schranke aufgestellt. 
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Es sei y;S7S7n, das gréBte, 7 enthaltende Intervall mit g(&,7)=G. Man be- 
trachtet ein Rechteckgebiet G': {f —e<&<€&, y,—e<yj<np+ €} mit dem Rand 
©’+ © nach 1.1, wobei e>0 so klein gewahlt ist, daB einerseits GW’ <G gilt und 
andererseits g(&, 7) — U?(&)<0 in W’ ist. Wegen GS0 ist das sicher mdglich. 
Dann ist in dem Urbild &’ von ©’ in der x, y-Ebene stets u(x, y)>0, also 
existiert g,,, nach (2.12) auf B’= G’+ ’+ €. Auf dem Rand W’ gilt nach Kon- 
struktion g(&,7)>G. Damit laBt sich das Maximum-Minimum-Prinzip auf ¥’ 
anwenden, mithin ist auch g(&,7)>G auf $8’. Dies steht in Widerspruch zu 
g(&,7) =G, also wird G auf dem Randstiick € von B angenommen. Also gilt 
auch hier (2.21) und man zeigt wie oben das Erfiilltsein von (2.20) auf ganz %. 


Bemerkung. Aus (2.20) folgen die drei Sonderfalle 
a)  u(%, 9) S UGB fire=20;, 
De) = Uae tar -U (a= Us), 
Cewayy eee. tur Ua) == (45). 


3. Die Schubspannung u, (x, ¥) 


a) Satz (Maximum-Minimum-Prinzip fiir w,). Es sei w€ 2 und ’ eine offene 
Teilmenge von & mit dem in I.4 und I.2 definierten Rand €’+€;. Dann nimmt 
u, seine Extremwerte auf 6’=©’+’+(, auf dem Rande ©’ an. 

Beweis. Es sei M= sup u,. Angenommen, die Behauptung sei falsch, d.h. 

auf B’ 
alle Punkte (x’, y’) mit u,(x’, y’)=M mégen in &’+, liegen. Man betrachtet 
x% = infx’ und eine Stelle y mit u,(x,y) = M. Es sei noch y,Sy< yp das gréBte, 
y enthaltende Intervall mit u,(x, y)=M. Zunachst werde M<O angenommen. 
Man betrachtet ein Rechteckgebiet ©”: {x —exix<x, yy—Eexy<y.te} mit 
dem Rand &’’+ Cj nach I.1, wobei e>0 so klein gewahlt wird, daB Gc G’ ist. 
Da auf 8” =6"+C"+C7 stets u,<0 und damit w=£0 ist, laBt sich dort die 
Crocco-Transformation (2.15) anwenden. Dadurch geht %”’ iiber in einen Be- 


reich 8” der o,o-Ebene und aus u,(x, y) wird durch (2.16) die Funktion t(@, o) 
erhalten. Auf dem Rand ©” von %” gilt u,<M, also ist auch t<M auf Ce 
Man betrachtet die Crocco-DGl (2.18) auf 8”. Da die rechte Seite von (2.18) 
wegen o=0 auf G+ CY definiert ist und fiir t, = 1,, =0 identisch verschwindet, 
da ferner t (0, 0) auf ®” stetig ist und T), T, und t,, dort existieren, kann auf B” 
das Maximum-Minimum-Prinzip von I.5 angewandt werden mit dem Ergebnis 
t<WM auf 8”. Dies gilt auch insbesondere in dem Bildpunkt (%,@) der Stelle 
(x,y), wahrend andererseits nach Konstruktion von (x,y) dort t=M gelten 
sollte. Aus diesem Widerspruch folgt die Behauptung fiir M<0O. 

Ist M>0, so wahlt man « so klein, daB u,>0 auf 8” ist und schlieBt genau 
so. Ist schlieBlich M—0O, dann wahlt man das Rechteck 8” so, daB auf ©” stets 
uy<0, also etwa uw, —6<0 bleibt. Verkleinert man %’’ zu einem Rechteck- 
bereich 8": {%¥—eSxSx<%, y,—eSySy.+ 6}, so lassen sich auf 8’ die 
obigen Uberlegungen wieder anwenden; sie ergeben w,<—o auf ganz 8’. 
x—>Xx gilt dann auch u,< —6 auf 8” in Widerspruch zur Voraussetzung. Genau 
so schlieBt man auch beziiglich inf u,. Damit ist die Behauptung vollstandig 
bewiesen. 
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b) Wahlt man 6’=G, so geht die vorstehende Behauptung iiber in den 
Satz. Esseiu€ &, M,Su'(y)SM,in0sSy< cound 4, (%, 0) SN in %S bee 


Dann gilt auf 8: 
M, <u, (x, y) S Max(M,,N), 


d.h. , nimmt seine Extremwerte auf dem Rand © von B an. 


c) Definition. Bei festgehaltenem FEK<K,, Hy > Sel EAS) die Anzahl der 
Extrema von u(x, y) in 0<y< oo. Ist u konstant in einem groBten y-Intervall 
<<, so wird dies als en Extremum gezahlt, und zwar nur dann, wenn 
u(%, y) Su (x, ;) fir yp-e<y<ot € ist bei geeignetem e>0. 


Satz. Fir w€2 und 4% %, gilt 
EG) = E(%,). 


Beweis. Fiir E(x») =o ist nichts zu beweisen, es sei also E(x) < co ange- 
nommen; auBerdem geniigt es offenbar, die Behauptung fiir «= x, zu beweisen. 
Zwischen zwei benachbarten Extremstellen von u(%9, y) hat u,(%9, y) dasselbe 
Vorzeichen; also ist E(x») gleich der Anzahl der Vorzeichenwechsel von , (%9, ¥) 
und das entsprechende gilt fiir E(x,). Man betrachtet auf den Halbgeraden x = %9 
bzw. «=x, die gréBten Intervalle iy bzw. i,, in denen u,=-0 ist (ihre Anzahl 
ist héchstens abzadhlbar unendlich). Es sei etwa i, ein derartiges Intervall mit 
u,<0 und @’ die gréBte i; enthaltende zusammenhangende Punktmenge mit 
Uy<0. Wiirde &’ kein Intervall von Gy enthalten, dann ware auf dem Rande W’ 
von &’ iiberall w,=0 in Widerspruch zu a) (auf dem ,,oberen“ bzw. ,,unteren“ 
Rand von © ist u,=0 bzw. u,=20). Also enthalt ©’ mindestens ein Intervall 
ip CC mit w,<0. Sind i, und ij’ zwei solche Intervalle mit u,<0 und W’, G” 
die zugehérigen zusammenhangenden Punktmengen, so kénnen ©’ und &” nicht 
dasselbe Intervall ig <€, enthalten. Denn ware das der Fall, so gabe es von tg 
zwei Polygonziige p’ und p” nach ij und i/. Im Gebiet zwischen diesen beiden 
Wegen ware nach a) dann w,<0, wahrend in Widerspruch dazu die Intervalle i, 
schon als maximal vorausgesetzt waren. Es seien also den beiden Intervallen 
i, und ij zwei punktfremde Intervalle ig und ij zugeordnet und es sei auBerdem 
My <0 zwischen ig und ig. Nach a) erkennt man, da8 dann auch w,<0 zwischen 
t, und ty gilt. Ausgehend etwa von ij betrachtet man nun das gréBte Intervall 3, 
mit ™, (x,y) <0. Nach obigen Uberlegungen ist jedem der Teilintervalle i), i%, ... 
mit ij’ Cp und w, (xo, y) <0 auf if? héchstens ein Intervall if, i/, ... zugeordnet. 
Ist 9, das kleinste, alle i?) enthaltende, abgeschlossene Intervall, so ist nach den 
vorstehenden Schltissen auch v,<0 auf §,. Denkt man sich zu jedem dieser 
(nach Voraussetzung nur in endlicher Anzahl vorhandenen) Intervalle §, das 
zugehorige Intervall 3, aus der Halbgeraden x=x,, yO weggenommen, so 
bleiben dort nur Punkte mit w,=0 tibrig. Hieraus folgt die Behauptung. 


4. Abloseprofile 


Satz. Es sei w€%, %(y)>0 in 0<y< co und u'(0)>0, sowie U'(x)=0 und 
ONe) =O aul xo 4S x7. 


Dann gilt u(x, y)>0 fiir »<x%<%,, 0<y< co sowie Uy (x, 0)>0 auf (x9, x1). 
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Beweis. Man betrachtet die v. Misessche Form (2.13) der GS-DGlen auf dem 
Bereich 8. Wegen der gleichmaBigen Stetigkeit von s(é (vr), (r, sal auf $* 
und wegen g(, cc) =0 gibt es zwei Zahlen « und G mit 0<e<U?(%), 0<G 


derart, daB g(&,7)<e bleibt fiir 729 (&)+G und alle &€ (x, x). Dabei sei G 
so gro8 gewahlt, daB noch 


U? (4%) mee 


&(%,) SU? (x) — aC Ge Silt fun Ose 7 SG: 
 Bedeutet ©, den nach I.1 gebildeten Rand des Gebietes 
Be: (MSS, mE) Sy Smo(é) + Gh, 
so gilt wegen 7 (% 9) =0 nach Konstruktion 


U2(&) — 
gn) S02) — TO * (n —ng(@)] auf Gg. (2.22) 
Leitet man die rechte Seite von (2.22) partiell nach ¢ ab, so gilt dafiir wegen 
U'=0, nSm+G, n= —%y=0 auf ganz Be die Abschatzung 


; U(é) U’ U2(é)—e » 
os2u@) ue —28 2 [yng] +7 =O. (2.23) 
Da die rechte Seite von (2.13) fiir g,,,—0 identisch verschwindet, kann auf dem 
Gebiet Be auf die DGl (2.13) der zweite Randwertsatz von I.5 angewandt wer- 
den. Dabei sind die UnGlen (4.31) und (1.32) durch (2.23) und (2.22) zu ersetzen. 


Danach gilt dann (2.22) sogar auf ganz B_. Damit ist g(&, 7) < U?(€) fiir n>) (€) 
U2— 
und g,(&, (&)) S— ie 220) 


Ist nun “(x’, y’) =O an einer Stelle (x’, y’)€%, dann ist nach (2.11) an der 
zugehorigen Stelle (&’, 7’) der &, 7-Ebene auch g= U?; wie soeben gezeigt wurde, 
gilt dies jedoch nur fiir die Stellen (£’, 7). Wegen g, (€’, 7o(é’)) <0 und wegen 
(2.12) ist dort weiterhin wu, (x’, y’)>0. Wegen der Stetigkeit von u,, wegen u=0 
auf 8 und u(x, 0) =0 ist das nur méglich, wenn y’=0 gilt. Damit ist die Be- 
hauptung vollstandig bewiesen. 


Zusatz. Ist U’(x)>0 auf <%, %,>, so bleibt der obige Satz richtig, auch 
wenn v,(x) beliebiges (wechselndes) Vorzeichen hat. 
Beweis. Wahlt man G in dem Beweis des vorstehenden Satzes so groB, daB 
U2—e 
= ie 
OSUU'+ 6 
oben gefiihrt werden. 


no(&) gilt, so ist (2.23) wieder richtig und der Beweis kann wie 


5. Profilwendepunkte 
Definition. Die Anzahl der Extrema von U(x) in %<*<*%<¥%, werde mit 
G (x) bezeichnet!®, so daB also O<G(%) SG (%) S co gilt. Die Anzahl der Extrema 
von u,(%, y) in 0<y<oo mége W(x) genannt werden’; es ist OS W(x) So. 
Haben #'(y) und U(x) an der Stelle x =x), y=0 beide gleichzeitig ein echtes 


19 Zur Definition dieser Anzahl vergleiche man 3.c). Die Randextrema von U(*) 
und u,(%, y) werden nicht gezahlt. 
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Randmaximum oder ein echtes Randminimum™, dann soll. d=1, in jedemg! 
anderen Fall d=0 gesetzt werden. r 

Satz. Es sei v)(x)=0 und wE€&. In einer Umgebung %S*S%,, OS y<6 © 
der x-Achse existiere noch u,,,, und — eventuell bis auf die Ausnahmestellen 
x =y,; — auch w,,,; die Funktion u,, sei in jedem der Streifen y;S%* Sy;+1, 
O0<y<oo stetig. Die Funktion yy (1, =) sei auf 6* stetig. Fiir festes 


% CXxX, %> und O< Y< oo habe die Funktion uy (x, y) in O<ySY nur endlich 
viele Extrema. 


Dann gilt auf x<*7<% stets 
W(x) S W(%) + G(x) +4. (2.24) 


Beweis. Offenbar geniigt es, den Satz allein fiir x=, und W(x%)<~, 
G(x,)<.co zu beweisen. Zum Beweis wird jedem Maximum (Minimum) von 
uy (x4, y) entweder ein Maximum (Minimum) von , (%, y) oder/und ein Maximum 
(Minimum) von U(x) zugeordnet. 

Das erste Extremum von u,(x,, y) liege bei y,>0, es handele sich z.B. um 
ein Maximum. Wie in Nr. 3 soll das bedeuten: w,(x,, y) nimmt seinen Maximal- 
wert J; auf einem gréBten Intervall y;<y,<,y;' an, und es gibt ein 6>0 mit 
uy (%,, y) ST, fir yj —d<y<y+6. Man betrachtet fiir e,>0 das gréBte zu- 
sammenhangende Teilgebiet ©, CG, in dem u,(x, y)>T—e, gilt und das die 
Stelle (x, y,) enthalt. Die Zahl e, sei dabei so klein gewahlt, da 7,— e, gréBer 
ist als das nachste Minimum von u,(x,, y) und auch noch gréBer als u, (x, 0) 
bleibt. Aus 3.a) folgt, daB &, bis zum Rand € von & reichen muB, und daB das 
Maximum von #, auf ©, nur auf dem Rand © angenommen werden kann. Dies 
sei etwa an der Stelle (%,,.0) des ,,unteren Randes“ der Fall, es ist x9 54,<0%- 
Dann gilt offenbar noch in einer (eventuellen Halb-) Umgebung von *, die UnGl 
Uy (%,0)>1,—&. Es sei Xp Sxy<xX<X,<-x, die gréBte derartige Umgebung. 
Dann ist ,(%,,0)=4,(%,,0)=T—e. Wegen u,(x, y)Su,(%,,0) fiir alle 
(x, y)€@, ist ferner w,, (%,,0)<0 und auch w,,(%,,0)=0. Setzt man y=0O in 
(2.1) ein, so findet man die Gl. 


y ty, (x, 0) = — U(x) U'(x), (2.25) 


also ist U'(x,)20. AuBerdem ist wegen (2.25) und der Annahme Uyy (%,, Vy) 20 
aut 0S y= y, auch WU" (4) = 0: ue 
Es soll zunachst gezeigt werden, daB es in *,< «< x, mindestens eine Stelle x’ 
gibt mit U'(x’)<0. Ist U’(x,)<0, dann ist x’= x, eine derartige Stelle. Es sei 
nun U’(x,)=0 und es werde angenommen, es gebe keine solche Stelle x’ mit 
U"(x')<0. Dann ware also U’(x)=0 in <%,, 45. Nun ist 0<u (4, 0)<7Ts 
u,(%,, 0), also gibt es in ¢%,, ¥,> sicher Stellen mit P< Manvelenlee “ce 
derart gewahlt (wobei ¢, héchstens verkleinert werden mu), dab een 


» U(x) habe genau dann ein ,,echtes Randmaximum bei i 
‘ 2 a um" bei ¥= 4%), wenn es ein 
groBtes PUBS) %S* S*’<-%, und eine Zahl «> 0 gibt derart, dab U(x) =U (4%) 
gilt auf (%», ¥ > und dals U(*) SU (4) ist in 44% <x%’+ 68’. Der Fall Ce — Sone 
auf <%, %,> wird also nicht ais Extremum gewertet, d.h. dann ist d=0. Die Bés 


zeichnungsweise fiir ein ,,echtes Randmini st tae Gh 
; ” imum” von U(x) und fiir die F ion 
sei entsprechend festgelegt. (*) unktion %(y) 
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Uy (%,,0)<O0 ist und daB dabei x, nicht mit einer der Ausnahmestellen y, zu- 
sammenfallt. Es gibt dann eine von y,-Stellen freie rechtsseitige Umgebung 
%SxSx,+ 0 <x, derart, daB u,, (x, 0) wegen seiner Stetigkeit auf dieser ganzen 
Umgebung stets <0 bleibt; wegen wu, (x, 0) 20 gilt dann dort sogar wu, (x, 0)>0; 
mithin ist 0<,(x,+ 0", 0) <u, (x,,0). Man betrachtet einen in ©, verlaufenden 
Polygonzug p, der die Punkte (,, 0) und (x,, y,) verbindet. Das Gebiet ©j sei 
begrenzt nach unten durch y=0, nach rechts durch x =¥%,+ 0’, nach oben und 
links durch p. Da auf p stets u,=7,—e, gilt und da auf dem Rand y=0 von @) 
tiberall u,,<0 bleibt, wird das Minimum von wu, nach 3.a) auf & in (x,+0”, 0) 
erreicht, d.h. es ist ,,(%,+0, 0) 20. Nach (2.25) und wegen U’(x,+ 0’) =0 ist 
aber w,,(%,+ 9’, 0)S0, also gilt u,,(x,+0,0)=0. Leitet man (2.4) partiell 
nach y ab und setzt y=0, dann findet man w,,,, (x, 0) =0 fiir alle xE <x, x). 
Zweimaliges partielles Ableiten von (2.1) nach y gibt fiir y=0 die Gl 


% (x, 0) Wey (x, 0) = Uy yyy (x, 0). (2.26) 


Da u, auf ®&%, sein Minimum in (%,+ 0’, 0) annimmt, muB dort Uy yyy (X1+ O,-0) 20 
sein. Wegen %,,y(%+ 0", 0) <0 und u, (%,+ 0’, 0) >0 gilt jedoch w,,,,,, (4+, 0) <0 
nach (2.26). Aus diesem Widerspruch folgt die Behauptung, da es in <4, x1) 
Stellen ~ mit) U/(%’)<0 gibt. 

Daraus folgt einerseits: Liegt das Maximum von uw, auf &, auf dem ,,unteren 
mang = V— (0, dann ist Ur(7)3=0 aul <%,%,>. Andererseits: gilt; Ist, 7,4, 
dann hat U’(x) an einer Stelle xf ein Maximum. Ist xj>x,, dann handelt es 
sich um ein inneres Maximum, also gilt G(x,)=1. Ist jedoch xs = xy, dann tritt 
ein echtes Randmaximum von U(x) auf; da in diesem Falle jedoch auch %(y) 
bei y=0 ein echtes Randmaximum besitzt, gilt dann d=1. 

Ist jedoch £,>% 9, gibt es also zwischen x) und *, Stellen mit u,,(x, 0)>0, 
so zeigt man 4hnlich, daB es dann links von %, Stellen x*’” mit U’(x’’)>0 gibt. 
Dazu wahlt man ¢;>0 so klein, daB u, (x%, 0) <u, (%,, 0) — e; ist und bezeichnet 
mit @’ dasjenige zusammenhangende Teilgebiet von ©,, das (*,, 0) enthalt und 
in dem stets w, (x, y) >, (#1, 0) —e1 bleibt. Ist (x1, 0) der linke Randpunkt von 
@’ fiir y = 0, so denkt man sich auBerdem ¢; noch derart gewahlt, daB w,,(«, 0) >0 
bleibt. Ist dann #’>0 so klein gewahlt, daB noch u,,(x, 0) >0 bleibt in x SxS 
x, +8” (die y,-Stellen sollen wie oben ausgeschlossen sein) und bezeichnet man 
mit @ den Durchschnitt von ©’ mit dem Streifen x»<*<+«;+0", so kann 
man in ®, wie oben weiterschlieBen. Mithin nimmt auch in diesem Falle die 
Funktion U(x) an einer inneren Stelle xy € XX, %,> ein Maximum an, d.h. es ist 
auch hier G(x) =1. 

Ist der erste Extremwert von u,(x,, y) ein Minimum, so 1a8t sich ihm auf 
genau dieselbe Weise ein Minimum von U(x) zuordnen, auch in diesem Falle 
ist dann entweder G(x,)=1 oder d=1. Falls also das Extremum von w, in ©, 
auf dem unteren Rand von © angenommen wird, ist stets G(x) +d=1. 

Man fiihrt nun zu den nachsten Stellen (%,, y,) (v=2,3,...), im denen 
Uy (%,, 'y) einen Extremwert 7; annimmt, dieselben Uberlegungen durch, kon- 
struiert also zusammenhangende Gebiete ,, fiir die u,>7Z}—e, (Maximum) 
bzw. u,<T,+¢, (Minimum) gilt. Dabei seien die ¢, so klein gewahlt, daB die &, 
punktfremd bleiben. Im Falle eines Minimums bei y,>y, wird das Minimum 
von u, auf &, wieder auf © angenommen; es liege abermals auf dem unteren 
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Rand von G, etwa an der Stelle (%,, 0). Es gilt #,<%,. Genau wie beim ersten 
Schritt sieht man auch hier, daB U(x) an einer Stelle x7 <x, x¢ €<%_,%> ein 
Minimum besitzt, und daB dann G(x) +d=2 ist. 4 

Bei den Gebieten @, ®,,... fiihrt man dieselben Uberlegungen durch und 
findet Stellen x*<x*,, in denen U(x) einen Extremwert annimmt. Nach 
héchstens G(x,) Schritten kann daher die bisherige Annahme nicht mehr erfiillt 
sein, daB das Maximum (Minimum) von wu, auf &, an einer Stelle (%,, 0) ange- 
nommen wird. Tritt dieser Fall nach K Schritten ein, dann gilt G(x,)+d2K. 
(Fiir U’(x) =0 in <x», %> ist offenbar K =0.) 

Fir »=K+1, K+2,... werden dann die Maxima (Minima) von uw, auf den 
Gebieten @, im Innern der Halbgeraden x = %), y>0O angenommen, das bedeutet 
gerade, daB a’(y) dort ein Extremum hat. Dieser Fall kann also héchstens 
W(x) mal eintreten. Insgesamt gilt dann also die Abschatzung W(x,) [ W(x) + 
K<W(x,) + G(x) +d, womit die Behauptung (2.24) vollstandig bewiesen ist. 


6. Einzigkett 

Satz. Es sei w#(y)>0 in O<y<oo. Seien wEX und wE€& zwei Lésungen 
von R mit den Eigenschaften u(x, y)>0, u(x, y)>0 fiir x» Sx<x, und y>0; 
mit u(x, 0)>0 und w, (x, 0)>0 in »S*<%,; sowie u,,<C und u,,<C auf 8. 

Dann ist uw=%a auf %, d.h. das Randwertproblem F hat héchstens eine ein- 
zige Lésung u€ 271, 

Beweis. Wegen u(x, y)>0" auf G': {xy9< x< %_< 4%, O<y< oo} ist dort die 
v. Mises-Transformation umkehrbar eindeutig, als Umkehrfunktion von € und 7 
nach (2.9) findet man 


FCS erat ON (Ey) 


nol& 


(2.27) 


dt 
é VO2(® —ge(& 2) 


Zum Beweis der Behauptung werde der Allgemeine Einzigkeitssatz von I.6 fiir 
m= 1 und die v. Misessche Form (2.13) der GS-DGl unter den Randbedingungen 
(2.14) auf dem transformierten Bereich 8’ benutzt. Die in I.6 verwandte Funk- 
tion f ist also 


0) fie ge = 
eran one | 


vVU2(é)—ge,, fir g<U?, 
Fiir die in I.6 benutzte Funktion @ setzt man 

p(é,n, g,h, w) = VU?— g —- VU? -h—w. (2.29) 
Sie ist stetig fiir — 0o<g, h<U?(&), —co<w<+oo. Fir g, h< U?() gilt weiter 


(2.28) 


on Ui Lai eat ie 
VU?—g-JU2—h 

beth hy 1 
2/U2— ¢’ On 


 ] 


(2.30) 
aaa SS ——- > Q, has . 
aVOPaA Ow —ai()) 


Ka : " ‘ : he: 
; Man vergleiche dazu H. GOrTLER [6], der als erster einen derartigen Einzig- 
keitssatz angab. 


a ce Die folgenden Aussagen werden stets fiir formuliert, sie sind genau so giiltig 
ir %. 
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Als Auflésung von o(&, 7, g, h, w) =0 nach fh findet man 


o(é,n, g,w) =g + 2w)U?— ¢— wv 
mit 


VW—o =|U?—g—w. 
Weiter werde gesetzt 


wo(f,n,g,h,w) =yeg (2.31) 


mit einer Konstanten y> 0, iiber die spater noch verfiigt werden soll. Als Folgen 
«, und f, werden 


1 é 


gewahlt. Damit ergibt sich als Lésung der DGI (1.36) unter den Randbedin- 
gungen (1.37) die Funktion 


c(E.ni = B=— Se er 4) b= 1,2,..) 
mit «<0 auf B, ¢,=¢,,=0 und jim e(7, £5) =7=0. Setzt man ¢ in o ein, 
so findet man die Apiotaneen fie 
d € 
dn o(, 1, §) 2 Engg 6? 
ae 2(U? ao Te gee) 
o(é, Up §,€ é) = San a Bn 


2/U?—¢ 
Damit sind alle tiber 9 und w getroffenen une sen des genannten Satzes 
von I.6 erfiillt, es bleibt nur noch die Gl (1.38) nachzupriifen. 

Zur Abkiirzung werde die linke Seite von (1.38) mit L bezeichnet: 


B= 06 (S518 70,'8) +0¢7 (E57 Bs Gis Sua) + Ont (E: 9s Os Ons Onn) =P Ow (E57, 29 Ss Ena) 


Dabei sind in e,, @, und e, dieselben Argumente wie in g; einzutragen, weiter 
tego (¢, 7); 8, «) und «== 6(€;%);.k). Mit (2.30), (2:31) smd; (2.32),ergibt sich 
dann 


, U = == /U2 —g 
ee OE ip Gee =H 5 Onn Cis VE 
ee D2 U? 8). nae | 
| 
Negeri Sega tia: =e cara 


Macht man nun die v. Mises-Transformation gema8 (2.11) und (2.12) wieder 
riickgangig, so findet man 


yuz—uU|U'|—vuC , | 
; ag ey 4 


=|e| [A(z 9) +7] 
mit 
A(x, y) = ee —uU|U'| —vuC]. 


Wegen u(x, y)>0 fir y>O und der Stetigkeit von u und w, auf B ist A(x, y) 
gleichmaBig beschrankt auf 8’—G’, wobei ©’ einen Sess Koss Maas 
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O0<y<6 bedeuten soll. Wegen » us, (x, 0) >0 und u(x,0)=0 und der Stetigkeit 
von w und w, kann dabei 6 so klein gewahlt werden, daB auf GS’ sogar A(x, y) 20 
gilt. Also ist A(x, y) auf B’ gleichmaBig nach unten beschrankt, es sei dort etwa 
A(x, y)=C’. Setzt man nun y= —C’, dann gilt auf ganz %’ und damit auch 
aut B’ stets LO und dies war gerade die letzte noch fehlende Voraussetzung 
(1.38) zum Allgemeinen Einzigkeitssatz von I.6. Damit gilt dessen Aussage, 
d.h. die v. Misessche DGl (2.13) hat unter den Randbedingungen (2.14) hochstens 
eine Lésung auf %’. Macht man nun die v. Mises-Transformation wieder riick- 
gaingig, dann findet man w=w auf BW’, dh. die Einzigkeit einer Lésung der 
GS-DGlen (2.1) und (2.2) unter den Randbedingungen (2.3) bis (2.6). Da diese 
Aussage fiir jeden Bereich Q’ gilt, ist sie auch noch richtig fiir 7.2%, d.h. auf 
ganz %, und dies war gerade zu beweisen. 


Zusammenfassung 


Es sollen die Ergebnisse der Nummern 2 bis 6 in Worten und Skizzen veran- 
schaulicht werden unter Benutzung der Sprechweise der Grenzschichttheorie. 
In den beigefiigten Bildern stellt jeweils der linke Teil das Anfangsprofil, der 
rechte ein stromabwarts liegendes Geschwindigkeitsprofil dar. Der durch- 
strichene Pfeil <> zwischen den beiden Teilen soll bedeuten: aus dem darge- 
stellten Anfangsprofil kann nicht das dargestellte rechte Geschwindigkeitsprofil 
entstehen. Aufgetragen sind jeweils als Abszisse u(x, y), als Ordinate y von der 
Wand y=0 an, wie das bei dem Anfangsprofil von Bild 10 vermerkt ist. 


Zu 2. Ubergeschwindigkeiten kénnen a) nicht ,,von selbst“ entstehen, sind sie 
jedoch entstanden, so kénnen sie b), c) nicht weiter anwachsen. Genauer gilt: 


a) Besitzt das Anfangsprofil keine Ubergeschwindigkeiten, so auch kein strom- 
abwarts liegendes Geschwindigkeitsprofil (vgl. Abb. 10). 


b) Besitzt das Anfangsprofil eine Ubergeschwindigkeit der GréBe e und 
nimmt die AuBengeschwindigkeit nicht ab (in Abb. 11 durch die Gestalt der 
Wand angedeutet), dann kann stromabwarts keine Ubergeschwindigkeit gréBer 
als e auftreten (Abb. 11). 


c) Besitzt das Anfangsprofil eine Ubergeschwindigkeit der GréBe e und 
nimmt die AuBengeschwindigkeit U(x) nicht zu, dann kann stromabwiarts nir- 
gends die Geschwindigkeit U(x ))-+e tiberschritten werden, wenn U(x») die 
AuBengeschwindigkeit im Anfangsprofil ist (Abb. 12) 


Zu 3. a) Die Schubspannung nimmt ihre maximalen und minimalen Werte 
im Anfangsprofil und an der Wand an. Das bedeutet also 


b) Ein ,,Aufsteilen“ eines Geschwindigkeitsprofils im Innern der Strémung 
ist nicht médglich (Abb. 13), und 


Monotone Geschwindigkeitsprofile bleiben auch stromabwarts noch monoton 
(Abb. 14). Weiter gilt sogar 


c) Die Anzahl der Extremwerte eines Geschwindigkeitsprofils kann sich strom- 
abwarts héchstens vermindern (Abb. 15) 


Zu 4. Kiickstromung. Ist U'(x)=0 (nimmt die AuBengeschwindigkeit strom- 
abwarts nicht ab) und wird an der Wand héchstens abgesaugt (wobei undurch- 
lassige Wand noch zugelassen ist), so kann 
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Abb. 10 Abb. 11 


7, Vp 
Abb. 12 Abb. 13 

Ge as iy Uh 
Abb. 14 Abb. 15 


Abb. 16 Abb. 17 
Yr UV Y 


Abb. 18 Abb. 19 


Abb, 10—19. Veranschaulichung der Ergebnisse 
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a) im Innern der Strémung die Geschwindigkeit « nirgends auf Null ab- 
sinken, d.h. dort kann keine Riickstrémung eintreten (Abb. 16), und 

b) auch die Wandschubspannung kann nicht zu Null werden, d.h. auch 
Stromungsablosung ist ausgeschlossen (Abb. 17). 

Fiir U’(x)>0, wenn also die AuBengeschwindigkeit stromabwarts dauernd 
zunimmt, darf an der Wand sogar noch beliebiges Ausblasen zugelassen werden. 


Zu 5. Profilwendepunkte. Undurchlassige Wand! Wechselt der Gradient U’ 
der AuBengeschwindigkeit zwischen einem Anfangsprofil und einem Endprofil 
G-mal das Vorzeichen, dann hat das Endprofil héchstens um G+ 1 mehr Wende- 
punkte als das Anfangsprofil. G+1 kann durch G ersetzt werden, wenn das 
Anfangsprofil an der Wand konvex ist und dort U(x) ein relatives Minimum 
besitzt, sowie wenn das Anfangsprofil an der Wand konkav ist und dort U(x) 
ein relatives Maximum hat (Abb. 18). 


Zu 6. Einzigkeit. Bei vorgegebener AuBengeschwindigkeit U(x) ist die ge- 
samte GS durch das Anfangsprofil eindeutig bestimmt (Abb. 19). 


Literatur 


[1] Pranptt, L.: Uber Fliissigkeitsbewegungen bei sehr kleiner Reibung. Verh. d. 
III. Int. Math.-Kongr. Heidelberg 1904, Teubner 1905, S. 484—494. Wieder 
abgedruckt in: L. PranpTtL-A. Betz, Vier Abhandlungen zur Hydrodynamik 
und Aerodynamik, Gottingen 1927, S. 1—8. 

[2] ScuLicutine, H.: Grenzschicht-Theorie. Karlsruhe 1951. 

[3] Nacumo, M.: Note in Kansu-HOteisiki No. 15 (1939) (in japanisch). Zitiert 

nach [14]. 

[4] WesrpuHaL, H.: Zur Abschatzung der Lésungen nichtlinearer parabolischer Dif- 

ferentialgleichungen. Math. Z. 51, 690—695 (1949). 

[5] GOrTLER, H.: Uber die Lésung nichtlinearer partieller Differentialgleichungen 

vom Reibungsschichttypus. ZAMM 30, 265—267 (1950). 

[6] Brown, W. B., & P. L. Donoucu: Tables of exact laminar-boundary layer 

solutions when the wall is porous and fluid properties are variable. NACA 

TN 2479 (1951). 

[7] Macer, A., & A. G. Hansen: Laminar boundary layer over a flat plate in a flow 

having circular streamlines. NACA TN 2658 (1952). 

[8] Macer, A.: Three-dimensional laminar boundary layer with small cross-flow. 

J. Aer. Se. 21, 835—845 (1954). 

[9] Loos, H. G.: A simple laminar boundary layer with secondary flow. J. Aer. 
Sc. 22, 35—40 (1955). 

[10] Herzic, H. Z., & A. G. Hansen: Experimental and analytical investigation of 

; secondary flows in ducts. J. Aer. Sc. 24, 217—231 (1957). 

pea ee a Eine einfache dreidimensionale Grenzschicht. ZAMM 36, 303—304 

1956). 

[12] GortLER, H.: Einflu8 einer schwachen Wandwelligkeit auf den Verlauf der 

laminaren Grenzschichten. Teil I, II]. ZAMM 25/27, 233—244 (1947); 28 

13—22 (1948). a 

[13] QUICK, A.W., & K.ScHRODER: Verhalten der laminaren Grenzschicht bei 

d periodisch schwankendem Druckverlauf. Math. Nachr. $,°217-—238 (4052)3 

[14] Sete ea Se ke as linégalité différentielle concernant les 

u type parabolique. Proceedings of the Ja 
| od von, q g Japan Academy XXVII, 
[15] Propt, G.: Questioni di stabilita per equazioni non lineari alle derivate parziali 
di tipo parabolico. Accademia Nazionale dei Lincei: Rendiconti 
scienze fisiche, matematiche e naturali) (8) 10, 365—370 (1951). 


(classe di 


[26] 


(27 


[28] 
[29] 


[30] 


Strenge Grenzschichttheorie 31 


NARASIMHAN, R.: On the Asymptotic Stability of Solutions of Parabolic Dif- 
ferential Equations. Journal of Rational Mechanics and Analysis 3, 303—313 
(1954). 

SZARSKI, J.: Sur la limitation et l’unicité des solutions d’un systéme non-linéaire 
d’équations paraboliques aux dérivées partielles du second ordre. Annales 
polonici mathematici II.2, 237—249 (1955). 


] Cotratz, L.: Fehlerabschatzungen fiir Naherungslésungen parabolischer Dif- 


ferentialgleichungen. Anais da Academia Brasileira de Ciencias, vol. 28, 
NOMENA 95 6: 


] Mrax, W.: Differential inequalities of parabolic type. Annales polonici mathe- 


matici III.2, 349—354 (1957). 
Kamke, E.: Differentialgleichungen reeller Funktionen. 2. Aufl. Leipzig 1945. 


] Grurano, L.: Sull’unicita della soluzione per una classe di equazioni differenziali 


alle derivate parziali, paraboliche, non lineari. Atti Accad. naz. Lincei, 
Rend., Cl. Sci. fis. mat. natur., (8) 12, 260—265 (1952). 

Man beachte die dazu im Referat Zbl. f. Math. 47, 92/93 (1953) von J. SzArskI 
gemachte Bemerkung. 


] Yosrs, T.: Uber die Unitat der Losung der gewohnlichen Differentialgleichung 


erster Ordnung. Japanese Journ. of Math. 2, 161—173 (1925). 


] SANSONE, G.: Equazioni differenziali nel campo reale. Parte seconda. Seconda 


edizione. Bologna 1949. 


] Saté, T.: Contribution a l’unicité de la solution d’une équation différentielle 


ordinaire. Japanese J. Math. 13, 1—6 (1936). 


] ZwiRNneER, G.: Criteri di unicita per gli integrali delle equazioni differenzial del 


primo ordine. Rend. Sem. mat. Univ. Padova 19, 273—293 (1950). 

RHEINBOLDT, W.: Zur duBeren Randbedingung bei den Grenzschicht-Gleichun- 
gen. 50 Jahre Grenzschicht-Forschung, S. 328—333. Braunschweig 1955. 
Siehe auch RHEINBOLDT, W.: Zur auBeren Randbedingung bei den Grenz- 
schichtgleichungen. ZAMM 36, 153—154 (1956). 

NIcKEL, K.: Die 4uBere Randbedingung der Grenzschicht-Differentialgleichung. 
Erscheint demnachst in der ZAMM. 

Misss, R. v.: Bemerkungen zur Hydrodynamik. ZAMM 7, 425—431 (1927). 

Crocco, L.: Lo strato limite laminare nei gas. Monografie scientifiche di aero- 
nautica, Numero 3, Dicembre 1946, Roma. 

Wirtinc, H.: Uber die Instabilitaten der Prandtlschen Grenzschichtgleichungen. 
50 Jahre Grenzschicht-Forschung, S. 334—342. Braunschweig 1955. 


Institut fiir Angewandte Mathematik, 
Technische Hochschule Karlsruhe 


(Eingegangen am 8. April 1958) 


Die Zustandsgleichung 


der klassischen statistischen Mechanik 


RUDOLF KURTH 


Vorgelegt von C. TRUESDELL 


§ 1. Das Problem 


KuyincuIn hat in seinem Buche ,,Mathematical Foundations of Statistical 
Mechanics‘‘ [2] eine strenge Herleitung fiir die Zustandsgleichung 


(1.1) P=SPTP 


eines einatomigen idealen Gases gegeben. P bedeutet hierbei den Druck, k die 
Boltzmannsche Konstante, I die Temperatur und D die Anzahl-Dichte des 
Gases. Die wesentlichsten Voraussetzungen jener Herleitung sind: Die Summe 
aller Wechselwirkungsenergien zwischen je zwei Partikeln ist im Vergleich zur 
Gesamtenergie des Systems sehr klein, und das System befindet sich im Zustande 
statistischen Gleichgewichts. Letzteres besagt: Seine statistischen Eigenschaften 
kénnen mittels einer zeitunabhangigen Wahrscheinlichkeitsverteilung beschrieben 
werden. 

Was wird nun aus der Zustandsgleichung (1.1), wenn man jene Voraus- 
setzungen (die in KHINCHINs Theorie methodisch durchaus wesentlich sind) fallen 
laBt? Diese Frage zu beantworten, ist die Absicht der vorliegenden Note. 

Wir gehen hierzu von den folgenden Voraussetzungen aus: Das mechanische 
System besteht aus m Massenpunkten (,,Partikeln‘) mit den Massen m1, m?,...,m”. 
Ihre Ortsvektoren, bezogen auf ein System rechtwinkliger kartesischer Koordi- 
naten 9,, 92,93, selen 


(1.2) 7 = (01, 92, 93)» 
ihre Impulsvektoren 
(1.3) r= (fi, p>; Ps) « 
DaB das System mechanisch abgeschlossen sei, wird nicht vorausgesetzt. Die 
auf die einzelnen Partikel wirkenden Krafte kénnen ganz beliebig sein. — Wir 
fiihren nun die folgenden Abkiirzungen ein: 
— 1 2 wee ie 
(1.4) PP Pee Py 


g=(7,9,.--.9") 
sowie 


(1.5) x = (p,q). 
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Die 6n-tugel x bezeichnen wir als die Phasenpunkte des Systems und die Menge 
all dieser 62-tugel als seinen Phasenraum J’. Ferner nehmen wir an: Die Wahr- 
scheinlichkeit dafiir, daB der Phasenpunkt des Systems zur Zeit ¢ irgendeiner 
bestimmten Borelschen Teilmenge B von J" angehért, kann durch ein Integral 


(1.6) f w(x, t)dx 


iiber eine Wahrscheinlichkeitsdichte w(x,¢) dargestellt werden, die iiber I” 
Lebesgue-integrierbar ist, wobei 

fw(x,)dx=1. 

whe 
Ohne Einschrankung der physikalischen Allgemeinheit diirfen wir annehmen, 
daB w(x, ¢) nur in einer beschrankten Teilmenge von Null verschieden ist. 

Um unsere Frage zu beantworten, haben wir die Begriffe von Dichte, Druck 

und Temperatur zu definieren, und zwar in Analogie und Verallgemeinerung zu 
den entsprechenden Begriffen fiir Systeme im statistischen Gleichgewicht. 
Sodann haben wir zu sehen, ob und wie sich die Zustandsgleichung (1.1) mit 
diesen Begriffen verallgemeinern 1aBt. 


§ 2. Definition der Dichte 


Man wird die (Anzahl-) Dichte D(q°,#) des Systems im realen Raume I}? = 
{9°} = {(@2, @8, 98)} zunachst so erklaren wollen, daB das Integral 


(2.1) ape t) dq® 


iiber die beliebige Borelsche Teilmenge B? von J7° stets die Anzahl aller zur Zeit ¢ 
in Be befindlichen Partikel liefert. Indessen sieht man sofort, daB es eine solche 
Dichte D nicht geben kann: Die Mengenfunktion (2.1) hangt von B} stetig, die 
Anzahl dagegen hangt von Bf unstetig ab. Man wird daher die Anzahl durch 
eine geeignete stetige Approximation zu ersetzen suchen, z.B. durch ihren 
Erwartungswert. 


Hierfiir fiihren wir die folgenden Abkiirzungen ein: 


Tie his: 
(2.2) T= {pr 
M={x}={(0".9)}, 
(2.3) i age) ee (aerate yet cas foi") > ASW 


Ferner sei Bi, die zu BY kongruente und kongruent liegende Menge in I; und B’ 
die iiber BY in I’ errichtete Zylindermenge, so daB also B’= By x (I/I7). (Das 
liegende Kreuz bedeutet Mengenmultiplikation.) SchlieBlich bezeichnen wir 
charakteristische Funktionen stets mit dem (,,Kronecker“-) Symbol 6, z.B. die 
charakteristische Funktion einer Teilmenge M von I’ mit 6, (x); sie ist definiert 
durch die Gleichungen 

4, wenn «cM, 


2. Ona) = 
ee a (%) 0, wenn x«C€l'—M. 


Arch. Rational Mech. Anal., Vol. 2 


Ww 
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Wir erklaren nun in J" eine Funktion N(x) durch 


(2.5) 


Offenbar ist N(x) bei beliebig vorgegebenem Phasenpunkt x die Anzahl der 
Partikel in BY. 
Ferner sei ‘ 
(2.6) wi, (q”, t) = J w(x, t) d(x/¢’) 
HIT. 
die Wahrscheinlichkeitsdichte des y-ten Partikels in seinem Ortsraum J’. Aus 
(2.5) und (2.6) folgt unmittelbar: 
(2.7) Npo(t) & f N(x) w(x, t)dx =f Dd) wy (9, t) dq. 
e Bo y=1 
Das heiBt: Die Anzahldichte D(q°, 4) im Punkte g® zur Zeit ¢ wird durch 


n 


(2.8) D(q?, t) = es (7°, ¢) 
dargestellt. Sie ist ,,fast iiberall‘‘ eindeutig bestimmt. — Ahnlich erhalt man fiir 
die drtliche Massendichte: 


(2.9) 0 (9°, t) = Dm’ w4 (9°, t). 
y=1 


v 


Die Ausdriicke (2.8) und (2.9) stimmen mit denen von NOLL iiberein [6]. 


Im vorangehenden haben wir angenommen, daB die momentane Anzahl in 
jedem Augenblicke durch ihr Phasenmittel hinreichend approximiert wird. Unter 
welchen Voraussetzungen das nun tatsachlich auch zutrifft, ist nur in einem 
speziellen Falle bekannt — nur bei Systemen, deren Wahrscheinlichkeitsdichte 
zeitunabhangig ist und allein vom Energieintegral abhangt [4]. Dariiber hinaus 
kann man folgendes sagen: N(x) ist nach (2.5) eine ,,Summenfunktion‘ [2], 
d.h. eine Summe von m Gliedern, deren y-tes Glied allein von x” abhangt und 
einen positiven Erwartungswert besitzt. Ist nun zur Zeit t=0 der Korrelations- 
koeffizient je zweier Funktionen Ope (g“) und Ony(q’), =-v, sehr klein (in der 
Ordnung von 1/m), so ist er es auch in einer hinreichend kleinen Umgebung des 
Zeitnullpunktes. In ihr ist dann die Wahrscheinlichkeit, da8 die Funktion N(x) 
von ihrem Erwartungswert erheblich abweicht, ebenfalls klein. 


§ 3. Definition des Druckes 
Wir betrachten in einem Punkte g® des realen Raumes T;° ein kleines orien- 
tiertes Ebenenstiick S mit dem Flacheninhalt s und dem Normalenvektor 
% =(%,,%,%3). Wahrend des Zeitintervalls ¢...¢+At werden moglicherweise 
einige Partikel dieses Flachenstiicks nach derjenigen Seite hin durchfliegen, nach 
welcher der Normalenvektor weist. Jedes dieser durchfliegenden Partikel tragt 
einen Impuls mit sich, der eine Komponente in der (positiven) Normalenrichtung 
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besitzt. Der Druck P, des Partikelsystems auf das Flachenelement in der Rich- 
tung « wird nun erklart als der Erwartungswert fiir die Summe all dieser Impuls- 
komponenten, reduziert auf die Flachen- und die Zeiteinheit. Dabei sind nach 
der Division durch s und At die Grenziiberginge s+0, At->0 gleichzeitig zu 
vollziehen. P_ wird analog definiert. Als Druck P schlechthin bezeichnen wir 
den Ausdruck 


(3.4) P ee PEA Pit Pee Ps 


Die Definition ist von KHINCHIN iibernommen [2], aber in abgeadnderter Weise: 
KHINCHIN nimmt fiir S ein reales materielles Ebenenstiick (bei uns ist es nur ein 
vorgestelltes ideales), an dem die Partikel elastisch reflektiert werden. Daher 
istepeL thin ===. Unda? 32.2 & 

Um P, und P. zu berechnen, kénnen wir unter zwei Abanderungen KHINCHINs 
Gedankengang einfach iibernehmen. Diese Abanderungen sind: Die Wahrschein- 
lichkeitsdichten w*(x’,/) in J” sind bei Kuincurn vom Boltzmannschen Typ, 
wahrend bei uns keine spezielle Form vorausgesetzt wird; und insbesondere 
brauchen fiir uns die Impulsverteilungen nicht mehr kugelsymmetrisch zu sein, 
so daB P,, P. und P durchaus von der Richtung « des Ebenenstiicks S abhangen 
k6nnen. 


Die Ergebnisse sind: 


23> we | BY (6.9.1) ap", 


sass _{ rie *(p", g°, t) dp’. 


no 
i 
Me 
iM- I 


I 
fan 


(3.2) 


a 

I 

| 
- 


ll 
ran 


Hierbei bezeichnen H,' und Hy die beiden Halften des Impulsraumes /?, in die 
I? durch eine Ebene durch den Nullpunkt mit dem Normalenvektor « zerlegt wird ; 
und zwar ist H,* diejenige Halfte, in die der Normalenvektor « hineinzeigt. 

Die Abhangigkeit der Ausdriicke P, und P_ von « kann also recht kompliziert 
sein, weil ja auch die Integrationsgebiete H,” und H, von « abhangen, so daB 
sich im allgemeinen auch die Integrale mit « andern. Der Gesamtdruck P= 
P.+|P.| dagegen ist einfach eine quadratische Form 


3 3 
(3.3) Pp =>) > PB; (9°, t) Oy 


in g, deren Koeffizientenmatrix 


(3.4) Be) 


yf Pe (p", 9, t) dp’ 


wir als den Drucktensor bezeichnen. Wegen 
n 3 As 
1 v y vi bY 0 v 
P=> 5 [(d Pai) doa) ee" #0 de 
p 


ist die Form (3.3) positiv definit. 
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Den Ausdruck (3.4) fiir den Drucktensor kénnen wir noch etwas anders dar- 
stellen, indem wir den Begriff des értlich bedingten Erwartungswertes einfiihren. 
Diesen erklaren wir fiir irgendeine integrable Funktion /(p’, @’, ¢) durch den 
Ausdruck 

LF (B, 9°, t) w (Bg, t) ap” 
ay rT? 
oo MO = Tarr 0 ae 


ae 
rT, 


Den Nenner hatten wir in §2 mit w%(g®,?) bezeichnet. Wir kénnen daher 
schreiben: 
Meiers ck 
(3.6) BAO) = Doge PPh DD 
Bisher haben wir das Ebenenstiick S als unbewegt angesehen. Jetzt nehmen 
wir an, es bewege sich mit einer (momentanen) Geschwindigkeit v = (v1, v2, vs). 
Den Druck definieren wir auch in diesem Falle relativ zu S; daher haben wir 


den Ausdruck 


ms p° fiir die Geschwindigkeit eines Partikels relativ zu S jetzt 
m 
durch (—; pv) zu ersetzen. Somit ergibt sich aus (3.6) fiir den Drucktensor 


“ 1 v v Vv v vy 
(3.7) Pi = > 5 (BR mv) (Bf — mv) wy (9, 0). 
v=1 


Noir gibt in [6] einen zu (3.5) aquivalenten Ausdruck fiir den ,,kinetischen 
Anteil des Spannungstensors’’, zu dem er aber noch einen Wechselwirkungsanteil 
hinzufiigt. Wir ziehen es hier vor, diesen Wechselwirkungsanteil den Kraften 
zuzurechnen (vgl. auch [5]) und den Drucktensor allein durch den Erwartungs- 
wert des Impulses zu definieren, weil gerade bei dieser Definition sich die Zu- 
standsgleichung leicht verallgemeinern laBt. 


In Wirklichkeit wird ja nun der Druck nicht so gemessen, wie wir, sondern 
wie KHINCHIN ihn erklarte: Die Partikel durchfliegen nicht ein ideales Ebenen- 
stiick, sondern werden an einem realen Ebenenstiick (sagen wir: elastisch) reflek- 
tiert und teilen ihm dabei einen Impuls mit. Nicht P,+|P_| wird also gemessen, 
sondern 2P,. Der Unterschied P,—|P_| ist prinzipiell immer berechenbar, in 
praxi allerdings nur in Ausnahmefallen, z.B. bei Systemen im statistischen Gleich- 
gewicht (ftir die er verschwindet). Nort definiert den Druck mittels der Kraft 
auf eine beliebige, offenbar ideelle geschlossene Flache im Innern des Systems 
(der Flissigkeit). Wie KuiNncHIN ziehen auch wir es vor, phanomenologische 
und statistische Gesichtspunkte streng zu trennen. 


§ 4. Definition der Temperatur 


Fiir Systeme im statistischen Gleichgewicht ist der Erwartungswert fiir die 
kinetische Energie eines Partikels vom Orte unabhangig und fiir alle Partikel 
gleich. Bis auf einen konventionellen Faktor erklart man die Temperatur in 
diesem Falle durch jenen Erwartungswert. Wie ist nun der Temperaturbegriff 
bei allgemeineren Voraussetzungen zu definieren? Offenbar hat er unter anderem 
folgenden Forderungen zu geniigen: 
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(i) Fir Systeme im Gleichgewicht soll sich wieder die urspriingliche Definition 
ergeben. 


(ii) Die Temperatur soll im allgemeinen Falle eine Funktion des Ortes g® und 
der Zeit ¢ sein. (Bei Gleichgewicht ist sie von beiden unabhangig.) 


(iii) Die Temperatur soll dann und nur dann verschwinden, wenn die Ge- 
schwindigkeitsstreuung verschwindet. 


(iv) Die Temperatur soll mittels einer Summenfunktion dargestellt werden. 
(Grund: Messungen liefern nicht Erwartungs-, sondern Momentan- oder Zeit- 
mittelwerte. Fiir Summenfunktionen kénnen wir am ehesten die annahernde 
Gleichheit von Erwartungs- und Momentan- bzw. Zeitmittelwerten erwarten [4].) 

(v) Die Form der Zustandsgleichung (1.1) soll nach Méglichkeit gewahrt 
bleiben. 


Die Forderungen (i) bis (iv) sind zum Beispiel erfiillt, wenn wir 
eo oy th Pee ey 
(4.1) feed ath 2 we (p? — m’d) 
setzen, wobel 
er 
2 = — 
(4.2) 3 Ym 


die ,,mittlere Str6mungsgeschwindigkeit“’ der Materie des Systems im Punkte q® 
zur Zeit ¢ bedeutet. Die rechte Seite von (4.1) ist durch die (etwas vage formu- 
lierten) Forderungen (i) bis (iv) nicht eindeutig bestimmt. Fiigt man namlich 
noch passende Glieder hinzu, die zeitliche Ableitungen irgendeiner Wahrschein- 
lichkeitsdichte als Faktoren enthalten, so sind (i) bis (iv) immer noch erfiillt. 

Die in (4.1) definierte Temperatur ist richtungsunabhangig, wahrend der 
Druck (3.7) im allgemeinen von der Richtung abhangt. Soll eine Zustands- 
gleichung von der Form (1.1) gelten, so miissen aber beide GréBen entweder 
richtungsabhangig oder richtungsunabhangig sein. Man hat also noch entweder 
die Druck- oder die Temperaturdefinition abzuandern. 

Tatsachlich fiihren auch beide Wege zum Ziel. Erklart man das Richtungs- 
mittel der quadratischen Form }) P,, «;«; als den skalaren Druck J im Punkte q° 
zur Zeit t, so wird 47 


. it Saori Ls eee 
my —_ (p”— m’ 0)? wi (9°, t) , 


(4.3) y= 


und durch Vergleich mit (4.1) sieht man ohne weiteres: Sind die Ausdriicke 


_ (pr— mo)? , y=1,2,...,”, oder sind die Funktionen w;(9°, #), y=1, 2, ...,%, 
Mm 
simtlich von vy unabhangig, so gilt eine Zustandsgleichung von der Form 


(4.4) PR=kTD. 


Die beiden hinreichenden Bedingungen besagen: Die kinetische Energie der 
Individual-Geschwindigkeiten (Streugeschwindigkeiten) ist auf alle Partikel 
gleichverteilt (was fiir Systeme im statistischen Gleichgewicht tatsachlich zu- 
trifft), bzw. alle Partikel sind aquivalent. Letzteres soll heiSen: Ihre Massen 
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sind gleich, und die Anfangs-Wahrscheinlichkeitsdichte w (%, OV = 0 (04, 55.45% ,0) 
ist in den Orts-Impulsvektoren 2}, x?,..., x” symmetrisch. (Es ist dann nam- 
lich auch die Wahrscheinlichkeitsdichte w(x, ¢) in x1, x?,..., x” symmetrisch, 
und es sind alle Funktionen w%,(q°, t) gleich.) 

Die Symmetriebedingung fiir w(x, 0) im Falle gleicher Partikelmassen be- 
deutet praktisch nur eine unwesentliche Einschrankung, da im allgemeinen 
die Partikel sich nicht unterscheiden und daher auch in der gewahlten Funktion 
w(x, 0) nicht unterschieden werden diirfen. Das heiBt aber: w(x, 0) ist sym- 
metrisch anzunehmen. 

Anders allerdings liegt die Sache bei kristallartigen Systemen. Bei ihnen ist 
jedes Partikel in seiner Bewegungsfreiheit auf einen kleinen (abgeschlossenen 
zusammenhangenden) Raumteil beschrankt, und all diese Raumteile sind punkt- 
fremd. Zwar gilt dann die Zustandsgleichung immer noch, aber in der Summe 
fiir die Dichte D ist héchstens ein Glied von Null verschieden — die Partikel 
sind nicht dquivalent. In diesem Falle deckt sich der experimentelle Dichte- 
begriff nicht mehr mit dem der statistischen Theorie, wie er durch (2.8) definiert 
wurde, so daB nun die Erfahrung der Gleichung (scheinbar) widersprechen muB. 

Wir betrachten jetzt die zweite Moglichkeit, bei der wir Druck und Tem- 
peratur beide als richtungsabhangig ansehen. Die Zustandsgleichung muB8B dann 
die Form haben 


(4.5) P=kT;;-D, 
wobei die Temperatur 7 durch 

3 
(4.6) Ef = XT, C5 Ob; 


gegeben ist. Nehmen wir nun an, es seien alle Partikel aquivalent, so folgt aus 


(4.5) und (3.7) 
(4.7) AD eo 


und diese Definition (4.6), (4.7) wird man ganz allgemein beizubehalten suchen. 
Die Temperatur ist dann proportional zu der Komponente der mittleren kineti- 
schen Energie eines Partikels in der Richtung «. 

Sind zwar die Partikel nicht alle aquivalent, sind aber wenigstens die Aus- 
driicke 


1 
m?” 


vy 
a 


(pb; — m’0;) (bp; — m’0,) feed Neer 5 

samtlich gleich (d.h. ist fiir jede Richtung die entsprechende Komponente der 
kinetischen Energie auf alle Partikel gleichverteilt), so gilt wiederum die Zu- 
standsgleichung (4.5) mit (4.7). 

Die Voraussetzungen fiir (4.5) sind enger als die fiir (4.4). Daher gilt (4.4) 
stets, wenn (4.5) gilt, wahrend die Umkehrung nicht richtig zu sein braucht. 
In dem einzigen explizit bekannten Beispiel, dem des statistischen Gleichgewichts 
wird jedoch die Unterscheidung gegenstandslos, weil die Impulsverteilung iets 
symmetrisch ist. 
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§ 5. Die allgemeine Zustandsgleichung 
In praxi sind gewohnlich nicht alle Partikel Aquivalent; jedoch gibt es dann 
meistens nur einige wenige Klassen von sehr vielen unter sich aquivalenten Par- 
tikeln. In einem solchen Falle haben wir 


1 
Le =e s janet 
(5.1) va i 
A 2 Alia A Ayes A 


wobei die oberen Indizes 4 die jeweilige Klasse bezeichnen und n”) die Anzahl 
der Partikel in der Klasse 4; ferner 


U ( 

T,;= ES JE 

(5.2) ete 
PEO | (A = ye) (pe! =a 

OS = m4) (P; m v;) (p; m v;) 

sowie schlieBlich 
pa$pe 
= 


(5.3) 


(5.4) BO = kT, D® 
und 
(5.5) PB, =i Ee dD”). 


Wir haben also nicht nur Partialdichten und -drucke eingefiihrt, sondern auch 
Partialtemperaturen. Gesamtdichte und Gesamtdruck sind die Summen der 
Partialwerte, wahrend die Gesamttemperatur als das (nach dem Umfang der 
Klassen gewichtete) Mittel der Partialtemperaturen erklart wird. Mit diesen 
Definitionen gilt die Zustandsgleichung (1.1) fiir jede Partikelklasse und fiir jede 
Richtung. 

Eine Illustration dieser Verhaltnisse bietet das MilchstraBensystem: Fiir jede 
Massenklasse von Sternen gibt es in ihm eine richtungsabhangige Partialtem- 
peratur. (Vgl. z.B. [Z].) 

Die Zustandsgleichung gilt weitgehend unabhangig von der Art der Krafte, 
die auf die Partikel einwirken. Auf Grund dieses Ergebnisses méchte der Ver- 
fasser seinen friitheren Zweifel widerrufen, daB die Zustandsgleichung auf das 
System der Atomkerne und Elektronen im Innern der Sterne anwendbar sei [3]. 

Allerdings ist bei den Anwendungen stets zu bedenken, daB der Druck in 
der Zustandsgleichung durch P,+|P_| erklart ist und nicht durch 2P,, durch 
den Impuls frei sich bewegender Partikel und nicht durch den auf ein festes 
Flachenstiick iibertragenen Impuls. Dieser Unterschied kann Diskrepanzen zwi- 
schen Theorie und Beobachtung zur Folge haben. 
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I. Introduction 


In 1866, when MAxweELt discovered that quadratures of the Maxwell-Boltz- 
mann collision integral could be effected only for a molecular model with an inverse 
fifth power repulsive force, he thought that experimental data implied such a 
force law [14, p. 41]. He disavowed this belief in the second of his famous pa- 
pers [15, p. 692] on the kinetic theory of gases, but retained the model because 
of its mathematical advantages. Upon the theoretical discovery of thermal dif- 
fusion by. Ensxoe [7] from the study of more general models, it was recognized 
that MAXWELL’s molecules fail to exhibit this important phenomenon. At the 
present time, when it is still not known whether the Maxwell-Boltzmann equation 
possesses any rigorous solution describing non-uniform states, it is recognized 
that no mathematical model for which approximate calculations have been made 


42 D. E. Jonnson & E. IKENBERRY: 


gives a fully satisfactory qualitative description of all phenomena. With the 
exception, however, of thermal diffusion, the Maxwellian model is on a par with 
other models for the qualitative description of such states. 

Extensive exact calculations for the Maxwellian model, without linearization 
of the collision integral, have recently been made by IKENBERRY & TRUESDELL 
[11; 20], who based their work on a principal theorem showing that the “spherical 
moments” P,,,,, of order g=27-+s, which are defined by 


(1.1) Pr png ite fre? oY, (C) Hoa e, 


where Y,(e) is a solid homogeneous spherical harmonic of degree s in the com- 
ponents of the molecular velocity e, satisfy equations of transfer of the form 


(1.2) lee See re, Cons Fone Gane 


Here the dots stand for certain linear and bilinear forms in spherical moments 
of orders 1,2,¢—1, 9, and g-+1, and their gradients, Cy,), is a constant, and 
Qers is a bilinear form in spherical moments of orders less than g. 

Earlier it had been shown by WANG CHANG & UHLENBECK [2/] that products 
of solid spherical harmonics and certain Sonine polynomials, special cases of 
which had been introduced into the kinetic theory by BuRNETT [4], are eigen- 
functions of the linearized collision operator for Maxwellian molecules. It is 
of course not entirely accidental that the eigenvalues of the linearized operator 
are identical, except in sign, with the C,,), appearing in the theorem (1.2) on 
which IKENBERRY & TRUESDELL based their work. The C,,,, determine the rate 
of approach to equilibrium, in both the linear and in the non-linear theories. 
TRUESDELL [20, Ch. VI] has shown that, for the rigorous theory of a Maxwellian 
gas, no spherical moment relaxes more slowly than the flux of energy*. 

Justification of the linearization of the collision integral is to be sought in 
the insight it affords toward treatment of the non-linearized integral and in the 
extent to which conclusions drawn admit experimental verification. Meaning of 
the results derived from the linearized equation is admittedly problematic. 

The earliest use of orthogonal functions for a series representation of the 
distribution function seems to be Jarr&’s [12] expansion in spherical harmonics 
with coefficients which are functions of the molecular speed**. An expansion in 
Hermite polynomials, which are orthogonal over all velocity space with respect 
to the Maxwellian distribution as a weighting factor, appeared with GrRap’s 
famous ‘‘thirteen moment approximation” [8]. Hermite polynomials were also 
used independently by IKENBERRY [/0] in an expansion of the distribution 


* For the problem of relaxation, the forms indicated by the dots on the left-hand 
side of (1.2) vanish identically ; when (1.2) is linearized in the moments, Q»2 ;\5 1s dropped, 
so that (1.2) may be integrated to give 


(*) Bs (2) == Iaeg((O)) e Carist 


By (5.25) follows the result asserted in the text. Thus for the linearized theory 
(Caria is the relaxation time of P,,,,, but for the exact theory this is not true in 
general, since (*) has to be replaced by a more complicated solution. 


*x The advantages of an expansion in an orthogonal set over a power series such 
as initiated by BortzMann are remarked by IkENBERRY [10]. 
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function in phase space. These polynomials seem to offer no particular advantage 
except when it is desired to emphasize spatial symmetry of the Maxwell-Boltz- 
mann equation. 

In an approximate treatment of the propagation of sound in a monatomic 
gas, in which they neglected spatial nonhomogeneity, WANG CHANG & UHLEN- 
BECK [21] used the eigenfunctions of the linearized collision operator for Max- 
wellian molecules in a series representation of the distribution function. Mort- 
SMITH [/7] found the same eigenfunctions useful in a discussion of special cases 
which include boundary conditions in the conduction of heat and in the flow of 
gases at low density. 

In obtaining relations of physical interest from the Maxwell-Boltzmann 
equation, one has at the outset a choice between BoLttzmann’s [3, pp. 184—194] 
method of introducing a series representation of the distribution function and 
MAXWELL’s method of considering the transport equations for certain moments 
[15, §§ 7—14; 3, pp. 115—153; 5, pp. 47—52]. The two methods, each of which 
requires the selection of a basic set of functions with which to work, are not 
unrelated. Following BoLTzMANN, one may substitute the assumed series re- 
presentation into the Maxwell-Boltzmann equation and obtain an infinite set of 
partial differential equations for the coefficients in the expansion. When an 
orthogonal set of functions with appropriate weighting factor is used for the 
expansion, the set of differential equations for the coefficients in the assumed 
expansion constitute a set of transport equations. 

In this paper, we consider only the linearized collision integral for Maxwellian 
molecules, but for it we construct an exact and general theory, using as our 
basic set of functions the eigenfunctions of the linearized collision operator for 
Maxwellian molecules. An approach similar to ours has been made by Mott- 
SmituH [17], but we employ the complete set of eigenfunctions, in order to describe 
azimuthal as well as radial and latitudinal variations in the velocity distribution. 

The external forces which we consider are those due to an electromagnetic 
field. Such fields are becoming increasingly important in the kinetic theory of 
gases because of their application in the development of “‘magnetic bottles” and 
in the study of hydromagnetic waves, plasma oscillations, and cyclotron resonance. 
In addition, the general theory of an ionized gas subject to a magnetic field is 
important in the theory of radio propagation [4, p. 322] owing to the part played 
by the ionized layers of the upper atmosphere. The theory is useful, too, in solar 
physics [5, p. 322] where the sun’s general magnetic field plays a vital role. 

The standard form of the linearized collision integral which we obtain has 
been given by CHAPMAN [5, p. 129] and Ensxoc [6]; however, their method is 
obscured by complicated notation and unnecessary intermediate substitutions. 
Our first objective is to obtain this form by a method which, although equivalent 
in principle to CHAPMAN’S, has a rationale which may be useful in further re- 
searches. Our next objective is to show that the eigenfunctions of the linearized 
collision operator for Maxwellian molecules are of the form obtained by WANG 
CHANG & UHLENBECK [27]. In their proof, doubly infinite power series are 
summed by interchanging the order of summation. Here we introduce a method 
which depends only on properties of the transcendental functions involved and 
avoids the use of series expansions. It is expected that this method may shed 
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light on the eigenfunctions of the collision operator for molecular models other 
than Maxwellian. These first two objectives are accomplished in Chapters II, 
be VV and V- 

Our final objective, which is twofold, is accomplished in Chapters VI and VII. 
In Chapter VI we represent the distribution function as a series of the eigen- 
functions of the linearized collision operator, and we give a physical interpretation 
of some of the coefficients in the expansion. To achieve this physical interpretation, 
we operate with mean values of the eigenfunctions. In Chapter VII we substitute 
this series representation of the distribution function into the Maxwell-Boltzmann 
equation. By equating coefficients of like terms we establish partial differen- 
tial equations for the unknown expansion coefficients. 

In Chapter VIII we derive various recursion relations needed for the develop- 
ments in Chapter VII. 


II. The Maxwell-Boltzmann Equation 
A. Basic Notation 


The kinetic theory of gases is based on a velocity distribution function 
F =F (a, €, t), representing the number density of molecules at the point x at 
time ¢ having velocity c. The function F(a, ¢, ¢) is defined so that the number 
of molecules in the range x to x-+da, having velocities between e and c+ de, 
at time #, is given by F(a, ¢€, ¢) dx de [5, Ch. 2]. 

Consistent with the definition of F(a, c, 4), the mass density o and the number 
density ~ of molecules in position space are given by 


(2.1) o=mn=m){Fde, 


where m is the mass of a molecule. The mean value Q of a quantity Q is defined by 
(2.2) nQ=f{OQOFde. 


The components #,;; of the pressure tensor and the components g; of the 
energy flux are defined by 


(2.3) — 
G@=0C,C=imfC,CFde, 


where the peculiar velocity C is the difference between the molecular velocity € 
and the mean velocity wu=e: C=ce—u. The hydrostatic pressure # and the 
temperature T are proportional to the trace of the pressure tensor: 
(2.4) P=nkl = sp; 
where k is the Boltzmann constant. 
The Maxwell-Boltzmann equation for F is [5, p. 46] 
(2.5) DF = — 
where 


(2.6) eee a A Ae 
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a,F 
at 


= f vl (v, y) (F' F*' — F F*) dQ. 


In these equations, a quantity associated with the second of a pair of molecules 
is denoted by an asterisk and quantities associated with the outcome of an en- 
counter are denoted by primes. Accents when attached to function symbols 
refer to their arguments: F’=F(e’), etc. The external force acting on a molecule 
is given by mf. I(v,q) is the differential cross-section corresponding to a 
turning of the relative velocity » =e*—ce over the angle @ into the solid angle 
dn’ in the direction of v’=e*'—e’. I(v, gy) is determined by the intermolecular 
force law. We may represent the geometry of an encounter as in Fig. 1, where 
g and ¢’ are polar angles for v’ relative to wv as axis. 


Line of recession — 


Effective cross-section of 
second molecule (considered fixed) 
Line of approach VA 
of first molecule 


& 


Polar axis 


Center of 
second molecule 


Fig. 1. The Encounter Diagram 


Let 6 and ¢’ be polar coordinates in a plane through the center of one molecule, 
perpendicular to the line of approach of the other molecule. If @ is the angle of 
deflection, then for spherically symmetric molecules, g is a function of 0 only. 
If the intermolecular force is a monotonic function of the distance of separation, 
this function can be inverted [8, p. 339] to give b as a function of gy, b=b(q). 


The differential element of area in the cross-section is b db deb -dp de’. 
Since 


(2.8) dn'=singdgde', 
we have 
(2.9) bdbdel —— 2 © any — Tv, ») dn’, 


sing dp 
where I(v, y) is the differential cross-section. The integration over dQ is quintuple: 


(2.40) dQ —de*dn’. 


B. The Boltzmann Collision Diagram 


For particles of equal mass m, the Boltzmann collision diagram [3, p. 19] 
takes the form given in Fig. 2, in which the relative velocity vectors v and v’ 
are diagonals of a rectangle. We take k, k,, n, and n’ as unit vectors in the direc- 


tions of e’—e, c*’—c, »=c*—e, and v'’=c*'—c’, respectively. 
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We define non-negative parameters « and «, by writing 


ce =c-+eak, 


(2°14) c*’ — Cc = ley 


Consistent with the equations of conservation of momentum and of energy, we 


obtain from Fig. 2 
ct = 8 ak, 


(2.12) ef Zea 
and 
(2.13) v=asecd = asec 0, =v’. 


a Ly 


Evidently the collision diagram 
is completely and uniquely de- 
termined by specifying either 
é, ef, and nn‘, or €, ¢ ,.and i 
or c, c*’, and n. 


C. The Linearized Collision 
Integral 


The collision integral (2.7), 
which appears as the right-hand 
member of the Maxwell-Boltz- 
mann integro-differential equa- 
tion (2.5), is frequently linearized 
by replacing F by A, (1+ Q) and 
neglecting quadratic terms in Q. 
Here Fj is a Maxwellian equilibrium distribution function: 


Fig. 2. The Boltzmann Collision Diagram for Particles of Equal Mass 


(2.14) hee Ae 


‘ 1 Dx + Ik ; i ar 
Since Ig [¢"’=) Fk," as a consequence of the conservation of energy in collisions, 
the result of the linearization is 

0,F 


(2.15) are tod (Oe 
where /(Q) is the “linearized collision operator’, 
(2.16) J(Q) = J oT(v, ¢) R*(Q' + Q*'— Q — Q*) dQ. 


Both Cuapman [5, pp. 129—133] and Ensxoe [6, p. 140] have given J in a 
standard form: 


(2.17) J(Q) = — Ky(e) Q — J K(e, e*) Q* Fy de*, 
where re 
(2.18) Ky (ec) = f v1 (v, vy) R* dQ 


and where the kernel K(e, c*) is symmetric in e and e* and is invariant with 
respect to orthogonal transformations in which c* rotates around e at constant 
angle X (Fig. 2). 
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It is readily seen from (2.16) and from the conservation of momentum and 
energy in collisions that the eigenvalue equation 


(2.19) Ip) = 4p 


has the quintuply degenerate eigenvalue 4=0, with eigenfunctions 1, e, c2. It 
has been anticipated by several writers for some time, without proof, that for 
Maxwellian molecules all eigenfunctions are of the form* 


(2.20) ple) = St 3 (ac?) Y.(e), 


where Y, (€) is a solid spherical harmonic of degree s, and S%, ; («) is a Sonine poly- 
nomial of degree m defined by [4, p. 123] 


‘ 


erties tame a2) Sx) a: 
m=0 


This has been recently proved by WanG 
CHanc & UHLENBECK [21, pp. 8—17]. 

As was stated in the introduction, 
our first two objectives are to obtain 
(2.17) and to show that the eigenfunc- 
tions of J are of the form given by (2.20). 
In obtaining (2.17), we make full use of 
the Boltzmann collision diagram (Figs. 2 
and 3). By using known integrals, we are 
able to obtain (2.20) without employing 
any series representations of transcen- 
dental functions. 


Fig. 3. Details of the Collision Diagram 


III. Derivation of the Enskog-Chapman Integral Equation 


A. Transformation of Velocity Differentials 


The Enskog-Chapman transformation of /(Q) is accomplished by specifying 
c, c’, and n for the part of J/(Q) that contains Q’ in the integrand and by specifying 
c, c*’, and n for the part that contains Q*’ in the integrand [4, pp. 129—133]. 
We therefore proceed to express d2 =dc*dn’ first in terms of de’ and dn, then 
in terms of de*’ and dn. 


From the Boltzmann diagram we have 
(3.1) c*=c+uvn. 
Hence, for fixed e, 
(3.2) dce*=vdvdn. 
From (2.11) and (2.13) we obtain, with ¢ fixed, 
de’=c2da dk =v? cos? 6 dvdk, 
63) dc*’ = of da,dk, = v? cos? 6,dv dk,. 


* The eigenfunctions ¢ and c? are linear combinations of functions of the form (2.20). 
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Also from the Boltzmann diagram, using m as polar axis, we see that 
dn'-=sing dg dé’, 

(3.4) dk =sind do dé’, 
dk, = sin 6,d6,de’. 

Since p= —26=26,, we find from (3.4) that 

(3.5) dn' = 4cos 6dk = 4cos 6, dk,. 

Combining (3.2), (3.3), and (3.5), we obtain 

(3.6) dQ = 4Asec2d de’ dn = 4sec? 6, dc*'dn. 


When we specify e, ce’, and n, we take k as polar axis. When we specify 
c, c*’, and n, we take k, as polar axis. From the details of the Boltzmann dia- 
gram (Figs. 2 and 3), we see that 


dn =sin6 do de, 


G.7) dn'=sin 6,d6,d¢,, 


where ¢ is the angle between the plane of k and e’ and the plane of k and n, 
and ¢, is the angle between the plane of k, and c*’ and the plane of k, and n. 
By means of the relations p=2—26=206,, we find from (3.6) and (3.7) that 


(3.8) afd = sec? 6 singde dgde’, 
dQ = sec? 6,sin pde,dq dc*’. 
B. The Enskog-Chapman Integral Equation 


The desired transformation of (2.16) is best accomplished by treating separ- 
ately the four integrals 


(3.9) H(Q) = — QJ vl(v, y) Fr dQ, 

(3.10) K(Q*) = — ful(v, y) Fo Q* dQ, 

(3.11) L(Q') = J vl(v, 9) FY Q'aQ, 

and 

(3.12) M(Q*') = Jv (v, g) BY O*' dQ. 

It is readily seen that 

(3.13) FQN EONS (O08) ar 1O,) cv Os. 


H(Q) becomes the first term in the right-hand side of (2.17), where 
(3.14) Ky (ec) =f v1 (v, y) Fo* de* dn’. 


K(Q*) requires no transformation since it is already in the desir 
y ed form. 


We proceed to consider transformation of L(Q’) and M (Q*’) separately. For 
the desired transformation of L(Q’), it is necessary to express Jo* as a function 
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of ce’ and n, since the integrand contains Q’ as a factor. We begin by observing 
from Fig. 2 that 


(3.15) c*=c’—ak-+on. 
Since n -k=cos 6 and v=«sec 6, this implies 
(3.16) c¥® — ¢/2+ tan? 6 + 2asecd(n-c') — 2a(k-c’). 


By a well-known theorem in spherical trigonometry, we obtain from the dihedral 
angle formed by ec’, k, and n (see Fig. 3) 


(3.17) n-c' = (e'-k)cosé + c'sindsin(e’, k) cose. 
From c’=c+ «k (see Fig. 2) we also obtain, by taking the cross-product e’ x e’, 
(3.18) a|sin(c’, k)| = c|sin(e, e’)| =c|sin X’|, 


where X’ is the angle between e€ and ce’. Finally, by combining (3.16), (3.17), 
and (3.18), we have 


8 12 2 2 , * , 
(3.19) i = fle h(c* +o? tan? 6+ 2cc’ tandsin X cos €) 


We substitute (3.19) and the first of (3.8) into the right member of (3.11) and 
perform the integrations in the order ¢, m, and ec’. Using a known formula [J6, 
p- 1323] which relates the Bessel functions /J,,(z) to the Fourier cosine coefficients 
of exp (7z cos e), we can perform the integration with respect to ¢ in the evalu- 
ation of L(Q’), for any spherically symmetrical molecular model. The result is 


L(Q') = 224A f dpdc' asin sect 6 I(«sec 6d, g) X 

.20 24 tan? 
B20) igre cos. osm ye ‘PRS Oy, 
where J,(z) = Jy(¢z) and 
(3.21) 7 = 2hetanoseco. 
Since « =|c’—e|, v=|c*—el, (e, c’) =X’, and (ce, c*) =X, we may replace ce’ 
by c* in (3.20) if at the same time we replace «a by v and X’ by X. The factor 
A exp(—hc’?) in the right member of (3.20) becomes Fo* when we replace c’ 
by c*. We then have 


L(Q’) = 2x f dy de* vsect dsin py I(vsec 6, y) X 


x 1g ies cososin XK). f* 8 oR OF. 


The transformation of M(Q*’) is started by observing from Fig. 2 that 


(3.23) c*¥ = c*'— ak, +n. 
Since k, -n=cos 6, and v=v’=«, sec 6,, we have, using (3.23), 


(3.24) c¥2 = c¥/2 4 g2 tan2d, + 20, sec d, (n+ e*’) — 204 (Key: €*’). 
Arch, Rational Mech. Anal., Vol. 2 4 
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By a relationship analogous to (3.17) (see Fig. 3), we have 

(3.25) n-c*' = (e*’- k,) cos 6, + c*’sin 6, - sin (c*", ky) cos &. 
From c*’=c-+a,k,, taking e*’ xe*’, we have 

(3.26) a, |sin (c*!, ky) |= cl sim Xe], ‘ 

where X*’=(c*’,c). Finally, from (3.24), (3.25), and (3.26), we have 
3.27) k*=Aexp —hA(c*’?+ af tan? d, + 2cc*’ tan 6, sin X*’ cos &). 


Upon substituting (3.27) and the second of (3.8) into the right member of (3.12), 
and integrating with respect to ¢,, we obtain 


M(Q*') = 227A f dp dc*' a sin psect 6,1 (a, sec 6,, y) X 


.28 . #12 2 2 , 
Be Kloet C0s0, siete. ee 
where 
(3.29) iia DO LAN SEe Ogs 


Since a,=|e*’— e|, v=|e*—el, (e, ec) = X*, and (e, e*) = X> we may replace 
c*’ by e* in (3.27) if simultaneously we replace a, by v and X*’ by X. We also 
replace p by x — q,. The factor A exp (—/hc*’?) becomes Jo* when we replace c*’ 
by c*. Noting from Fig. 2 that 0, is the same function of g, as 6 is of g, we drop 
the subscripts from 6, and g,. The result is 


M(Q*') = 22 f dy de* vsing sect 6 I (usec 6, —q) X 
x [on (kOhcososin X\cn ee ere 
Upon substituting into (3.13) the quantities H(Q), K(Q*), L(Q’), and M(Q*’), 


given by (3.9), (3.10), (3.22), and (3.30) respectively, and upon comparing the 
result with (2.17), we obtain 


(3.30) 


27 
3.31) K(e, e*) = 2x f v {1(v, p) — sect O[I(usec 6, y) + I(v sec 6,%—q)]| x 
Din 0 


< I(r c* cos 6 sin X) e~" oN singdg. 


From (3.21) it is seen that the argument of the Bessel function in the right 
member of (3.31) is equal to 2hcc* sin X tan 6. Now, noting that 


(3.32) uv? = |e* — c|? = ¢2 + c¥2 = 2¢c* cos X, 6=i(x= 9), 


we see that A(€, €*) is a function of c, c*, and X only, is symmetric in e and e* 
. . . “ 
and is invariant as e* rotates around € at constant angle X. 


IV. The Integral Equation for the Radial Eigenfunctions of J 


From the general theory of linear integral equations with symmetric, ortho- 
gonally invariant kernels, it is known that the eigenfunctions wy(e) of J are of 
the form [19, pp. 360—363 | 


(4.1) ple) = S(2)¥,(e), 
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where Y,(e) is a solid spherical harmonic of degree s. To obtain the integral 
equation which the radial part S(c?) must satisfy, we substitute S (c?) Y,(e) for Q 


in (2.17). Using c*, X, and * as spherical coordinates of c* with e€ as polar 
axis, we have 


(4.2) de* = dg* dX dc* c*¥*sin X = —dy* du dc* c*?, 

where 4=cos X. We perform the integrations over y* and X in the order 
indicated, viz. dy* then dX. The result which we shall obtain was proved by 
WANG CHANG & UHLENBECK [2/, pp. 7—12], but we present a more lucid deriva- 


tion avoiding the use of infinite series. Using MAXwELL’s integral relation [15, 
p. 688], we have 


27 
(4.3) of f ¥,(e*) dp* = 2xc** Y,(c) F.(u), 
0 


where P.() is the Legendre polynomial of order s [23, p. 303]. We obtain from 
(2.17) and (4.2), since K(e, c*) is independent of 9g, 


(ArAN cc? (Sw) == — 1K, (€) SAY 2a Y, if io Sea Fo) EMC INOGUN NIE ee 
Defining a linear operator fi, by vet 

@5) PR) — — PK le) f —20 fadc* ***G,(0, oF) BY, 

where ; 

(4.6) Glee) = J du P.(u) Keer), 


we see from (4.4) that S(c?) Y,(e) satisfies the eigenvalue equation J(y) =Ay if 


Ss 


S(c?) satisfies the eigenvalue equation 
(4.7) (SC) as 


In general, K(e, e*) is such a complicated function of c, c*, and X through 
its dependence on v that no further progress has been made in finding the radial 
eigenfunctions S (c?), except for the Maxwellian model and for a Lorentzian gas [1]. 
As MAXweELt discovered [14, p. 41], the function vJ(v, y) is independent of v 
for molecules with an inverse fifth power repulsive force. For this model it is 
shown in the next section that the radial eigenfunctions S(c?) are the Sonine 
polynomials S¥'., (hc?) defined by (2.21). 


For Maxwellian molecules, vJ(v, g) is defined as follows [//, p. 22; 14, p. 41]: 


i [G da 
(4.8) v1(v, 9) = Vege 
where «(y) is given in terms of a parameter # by 


at = 2cot?2%, 


and G is the force constant; K(k) is the complete elliptic integral of the first kind. 
4* 
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V. The Radial Eigenfunctions for Maxwellian Molecules 
For Maxwellian molecules, vI(v, y) is independent of v. Let us define func- 
tions 7,(g) and 7,(p) by 
T.(p) =vI(v, 9) sing, 


ee Tin = Th bee. 


From (2.8), (2.10), (2.18), and (5.1), we obtain 


(5.2) Ky (ce) =JS de dp T(g) J de* Fo* = 22009 J T,(y) 4, 
where 
(5.3) ny = f de* Rf = Alani). 


For K(e, e*) we have, from (3.31) and (5.1): 
(5.4) K(e,c*) = 22 f dp[T,(~) — Ty(¢) sec? 6 Ip (r.c* cos db sin X) e*" 49], 


after observing that v sec 6 I(uv sec 6, y) sin p= 7, (q). 


We may perform the integration over m in (4.6) by substituting (5.4) into 
(4.6) and making use of (3.21) and the known integral [/6, p. 1325]: 


fit eerx 7 (zsin isin. X) P(cosX) sin. 1 a 
0 


(5.5) : 
=i" (2.1/2) B™ (cos) Jni4(2)- 
On letting z= —irc*, 0=$n—6, cosX =p, m=O, and »=s, we obtain from 
(5.5); 
1 
(5.6) f P,(u) Ip(re* cos dsin X) ef *9X dy, — (= J) B (sin d) Tnug(rc*), 


—1 
where I,,(z)'==(— 2)" J, (72) 
Upon substituting A(c, e*) from (5.4) into the integral (4.6) for G, (GC e=)e 


substituting for v? from (3.32), and using (5.6) to perform the integration with 
respect to uw, we obtain 


G, (e, c*) = 2a fag * 
(5.7) 


0 
20 \t 
x 2 0: T,(v) — ( a] P,(sin 6) - Ty (p) sec? 6 Teae(vth) Geese Nae e 


We proceed now to show that for 


(5.8) Pa] K he) = ? ee et hol(1—x) = y's 


n=0 
we have 


(5.9) HR, (/) 7 Ps hee Sopa (A c*) a 
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where 4,,, is independent of c and x. Since §, is a linear operator we shall then 
have 


(5.10) DH, [ Sy He] a = Y Aye Shy (he?) x" 
n=0 n=0 

Equating coefficients of x”, we shall have 

(5.11) R, [$4 (oe)] = Ans Shy (he). 


To obtain (5.9), we substitute Fy from (2.14), G,(e, e*) from (5.7), and f from 
(5.8) into (4.5) and perform the integration with respect to c* by using the known 
integral [22, p. 394, eq. 4] 

fore) y —72/4 p2 
—p?t? wt Oe YT € 
(5.42) L,(rt)e t He perrat 


Here we replace ¢ by c* and »y by s +4. The result is 


s = Say ye f. 
ae RS (/) ¢ le) 


3 apa ~ he? tan? 6 
- ana f An > R(sin 6) T, (g) sec? 0 e4? 


i! pba pe oa 


sec? (1 — x sin?6). 


& h 
(5.14) p=h(t+>2— + tants) 7" 


By use of (3.21) and (5.14), the coefficient of —P (sin 6) - J(y) sec? d in the 
integrand on the right side of (5.13) may be written in the form 


- (Fees) 
(5.15) Fale eee Oe es o esiny eos?'d en ges? 
(Le att ip ae hh (4 — x sin? 6)S+4 


From this we see that the integrand in (5.13) is a function of and of x sin?6 = 
x cos?3. By use of (5.15), (5.13) becomes, after division by c’‘, 


hc) | = — Ke IE) — 
(5.16) I, [f (x, hc?) | te) Ila ) 


yaa 


By use of (5.3) and (5.8), (5.16) becomes 


— (Sessa) 
smn". (sino) 7, (eye" 7 re 
(14 — x sin? )s*? 


0: L(y) — dy. 


PO eae Bales ee) 
— 27 fiat T,(¢) —sin‘ 6 P(sin 6) f (x sin?6, hc) T,(y)] dg. 
Noting that 


(5.18) 65 L(y) =X 8 On Th (p) Sty (hc?) , 


m=0 
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and using (5.8), we see that (5.17) may be written in the form 


(5.19) we] Sty (he?) a] =D Ins Sea lh) 2, 


where 


(5.20) As = — Kole) — 2m f dp[02 On T(p) — (sin d)*”"** F (sin 6) 7, (9)]. 
0 
Substituting for K,(e) from (5.2), we have 
(5.21) Ams = — 2709 f [(1 + 8 On) L(g) —.(sin 6)?" P, (sin 6) 72 (y)] de. 
i) 


Replacing 7;(p) by its value from (5.1) and recalling that 6 = 3 (a — @), we obtain 
by elementary reductions, 


(5 2) Ans = — 27M) J [1 7 oY Om — P, (cos 39) (cos$@)?""* = 
OD, 0 : : se 
— P(sin3¢) (sing y)”""*] L(g) @¢@. 
We have therefore shown that 
(5.23) T[SS4 (0?) Y, (€)] = Ams Soa (he?) ¥,(€), 
where 4,,, is given by (5.22). It is readily found that 
(5.24) hoo = 401 = 0 = O- 


IKENBERRY & TRUESDELL [J/, pp. 21—23] express /,,, as a linear combina- 
tion of certain /,,9, enabling one to carry out the calculations for Maxwellian 
molecules by using only tables of circular functions and of complete elliptic 
integrals of the first kind. IkENBERRY & TRUESDELL also show that Amo is 
ultimately proportional to mt when m is large and that, except for the special 
cases in which 4,,,=0, 


(5.25) hg Dia de eeiay byes 


where « is the coefficient of viscosity [11, pp. 36—38]. 


VI. Expansion of F in Terms of Sonine Polynomials and Spherical Harmonics 
A. Formal Expansion of F 


Let us assume that there exists a uniformly convergent series expansion of 


the distribution function F in terms of Sonine polynomials and solid spherical 
harmonics of the form 


(6.1) te VEG W Set SH DINE Aa ager (c), 


ms Wins 
m=0 s=0 k=—s 


where 


6.2) Ym s(€) = Is — & +1) S% 4 (he?) VF (e) 


Y,' (c) = c P* (cos 8) e'*®, m=O, —sSkssinrseo. 
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Spherical coordinates of € with respect to a fixed rectangular coordinate system 
are taken as c,%,g. The Nj, are arbitrarily introduced normalization factors 
which will be determined so as to make 


(6.3) A®,= SF yi,.de=ny,.. 


The PR+*(u) are associated Legendre functions defined in terms of the Legendre 
polynomials P.(u) by [23, p. 323; 9, p. 99]: 


(6.4) 


The second of equations (6.4) becomes nugatory when |k|>s. We define 
Wms(€) =O when m<0, or s<0, or |R|>s. 

Since Y,*(e) is a solid spherical harmonic of degree s in €, we see from (5.23) 
that pk,, is an eigenfunction of J with eigenvalue /,,,, 7.¢., 


(6.5) T (wh s) Fe Apis Yess 


where 4,,, is given by (5.22). Equating P=j(1+Q) used in linearizing the 
collision integral (Chapter II) to the value of F given by (6.1), replacing /y by 
its value in (2.14), and dividing out the factor exp(—hc?), we have 


foe) Ss 


(6.6) ANON SO. DNAS hohe: 


m=0 s=0 k=—s 


Operating on (6.6) with the linear operator J and recalling that J(1)=0, we 
have from (6.5), 


(6.7) AU OM epi Mal Aig! he, sh 


From (2.14), (2.15) and (6.7), we obtain 


s 


(6.8) i Bey i De NE Age Aah: 
m=0 s=0 k=—s 


We proceed to determine the arbitrary constant multipliers Ni, to make 
(6.3) valid. Making use of known integrals for products of Legendre functions 
[23, p. 324] and for products of Sonine polynomials [8, p. 124], we find that 


i Yn Ynse  ° de 
(2s +1) (m+1) At? 2 


(6.9) 


Multiplying both sides of (6.1) by yj, and integrating term-wise over €, we 
obtain by means of (6.9), 


(6.10) [Fug do=Aze 
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provided we let 
(—1)* (2s +41) Pim +1) Wtt 
anT(stk+1)0(s—k+1)P(stm+ 3) * 


(6.41) Ni = 


Thus we obtain (6.3). From (6.11) follows 


(6.12) N,¢ = N;E,. 
Noting that 
(6.13) Ym = (—1)" Wns)* 


where the asterisk denotes ‘complex conjugate of’, we see from (6.3) that 
(6.14) Ans =(—1) (Ags)* 


since F is a real function of ¢. 


B. Interpretation of the Coefficients in the Expansion of F 


We proceed now to obtain Ws for certain values of m,s, and k, and by 
applying equation (6.3) to these special cases to give a physical interpretation 
to some of the coefficients in the expansion of F. From the generating function 
(2.21) or from a known explicit expression [5, p. 124] for the Sonine polynomials 
Si (x), we find that in particular 


Sn(x) = 1, 
(6.15) one) he 
Sh(*) =F (mn +1) (n+ 2) — (n+ 2) x42 x2. 


By repeated use of (6.2), (6.13), and certain Legendre functions as needed, and 
by considering that c, 8, y are spherical coordinates of €= (c,, Cy, C3), we obtain 
the following special cases: 


0) Deere 

Yoo 1; 

Lie tees 

Wo1 = °3,> 

ey es " 
Ya =A tte), 


OR Omen Ao 
Yio 2 he, 


(6.16) via = (3 — he’) cs, 


pit = + (§ — he) (cq. bic,), 
Wyse Bs tod ct 


Vox =+3(q+ Cy) Ce, 


Yor =3(CU—cg+ DUG; Cp). 


Linearized kinetic theory 57 


We now apply (6.3) to the quantities given by (6.16), making use of (2.1), 
(2.3), (2.4), and the equations €=u, e=C-+u. We obtain 


0 
Ayy= ee 
1 
oe 
Ue (tine 043, 


: { . 
Agi = +5, (@um+tiom), 
1 (3 
Alo = mdlex: — 3hp — gh), 


ENS 
(6.17) ae m 3 @ Uz — h(293 + 263; 4; + 3p Ug + 0, 42)| ; 


5 ; F 
Att = £ 5 |5 ol tim) — 2h +g) — 
— 2h(pi1 Et Pin) U; — 3h P(t tu) — hg u®(uy + in), 


at ate 30us — ou), 


salt 
m 
AR= +t = (bis tP es + 0% My + 10 Mn Ms) 
3 i 
m (P11 £ 21 Pia — Poe OU, — QUy + 279m Ue). 


VII. The Differential Equations for the Expansion Coefficients 
A. DF in the Absence of External Forces 
In this section we express the left-hand side of the Maxwell-Boltzmann 
equation DF, as given by (2.16), as a linear series of the p*,,, for the case of no 
external force. In this case (2.16) becomes 


F 


For the series representation (6.1) of F&, the two terms on the right in (7.1) are, 
respectively, 


OF na St PAms 
(7.2) Ot =e pa pa 2 Nn at Wms 
and 
ars —he as “ Ne dAne k 
(7.3) (Eo Pog ee a pa be ms Ox; Cems: 
Making the change of variables 
Vy = % +1 Xe, 
vst X3» 
and 
Ey =O + 1p, 
(7.5) &y = Cy — 1p, 


et C3, 
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we find that 
aA ears 
(7.6) bag = Om OF 


It is shown in Chapter VIII that 


1 
UO Qst1)h (mn 1)(s—k)(s—k 1) Vesa — 
(7.7) 
4| 
= (ict S das (s — k) (s—k —1) pitt + Aye — Ae seal, 
| Bs. 
(7.8) #2¥ms = (2s+1)h [(m 1) (6 42 Be) (Sap = 1) ea os = 
We 
1 at = 
: (m + Ss x) (4 k) (s t 1)yht + hyd — bye sta] 
and 
1 1 
(7.9) Ym = Tasqaye (lm +8 +5) (9 — 2) Yims a 


— (m -L 1) (s? << k?) Corres Sia Wenge as, - J ae : 


Making use of (7.6), (7.7), (7.8), and (7.9), we may write (7.3) in the form 


A k 
« {Se [(m +4) (8 — 8) (8 — #1) yi, a 
a Wm+1,s—1 


(m +s +>) (s—)(s—k — 1) yhtt + byt — Pina, s¢1| — 


7.10 Ams 
ely) = SSE [EA E+ 8) (ok Ay 


1 
= (m + Ss ate =) (s i k) (s + k =3 1) ety rs h te =a hYm—, 51] ary 
pA 
OV3 


! (mm +s + 5) (8? = 2) yh, a — 


= (4) (S? eB?) We Aral aig Te ae aa a |}. 


Making use of (6.11) we may write (7.10) in the form 


(7.11) CUR ny Wns Dol Am aes 


ms 


where Qp is a linear operator which will be given explicitly in (7.34) 


In the absence of external forces we have, from (7.1), (7.2) and (7.41) 


7.12 DFaeYy SS ye [24ms - 
oe eM 2 a Mee (Gr $2 ol4nd)| vo. 


Linearized kinetic theory 59 


B. Terms Due to an Electric Field 


We now consider an external force mf, due to an electric field E. We write 
(7.13) , mf: =E, 


where E is independent of e. In this case, the force term in DF becomes 
(7.14) fot 


If we make the change of variable given in (7.5) and define new force components 
by 
mé,=EH,+7E,, 


(7.15) m & = E,—iEs, 


me, = F., 
then (7.14) becomes 


OF 
(7.16) fe: VF = 6, 5. 
It is shown in Chapter VIII that 


4) —h oc? é rc co = 
(7.17) Gee Yims = oeqay Am ssa + (m +1) (8 +h) (8+ 1) ving, sa] > 


—hc® é a 
(7.18) 08, a Ys a (2s-+1) Lae h Wels (m is 1) (s k) (s k— Tip sani | , 


and 


(7.19) Gee Wms = Taspay Ln, s+ + (m +1) (s*— B?) Wnt, sa] - 


Using (6.1), (7.16), (7.17), and (7.19), we may write (7.14) as 


co wo 


Pun ABest Sie Sy Sik st al icin. 


m=0), s=0 k=—=s (2s+1) 


x {E, [hye dia + (m +1) (s + 2) (s + —1) veh, al — 
— Eh Ym ser + (m k—1) ys al + 
ae 2é;[— WPraty (m 1) Os kh?) Ponsa lt 


(7.20) 


—s 
SS, 
— 

wn 

=a 
SZ 
— 

“vn 


Making use of (6.11), we may write (7.20) in the form 
(7.24) fe ViF=e"" 3X YY Nye €lAns) Pns. 
m=0 s=0 k=—s 


where © is a linear operator which will be given explicitly in (7.35). 


C. Terms Due to a Magnetic Field 


We now consider an external force mf, resulting from a magnetic field H. 
The Lorentz force mf, is given in electromagnetic units by 


(7.22) mfg=eocxH, 
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where ¢é, is the charge on the molecule. We have for the force term in DF, 
(7.23) fy VF =—2H- (ex VF. 


Noting that the operator —i(exV,) is obtainable by transcribing into € 
space the quantum mechanical angular momentum operator * 


(7.24) M=rxp=—i(rxV,), 
whose properties are well known [2, Ch. 14], we write (7.23) in the form 
(7.25) fu VF = —i—2 (H-M)F 
after replacing r by ¢ in (7.24). Introducing the shift operators 
(7.26) M,=M,+71M, 
and hypercomponents 
(7.27) : 

—_ 0 

Hs ae Hz, 
of the magnetic field, we find that (7.25) may be written in the form 
(7.28) - fy VF = —i(#_M,+#,M_4+ #,M,)F. 
Using the coordinate representations 
ace HP 6 . alle 

Mas e (+ gp +icota ze), 


Soft 
~* ap? 


(7.29) 
M, = 


we find by means of (6.2) and (6.4) that 

M, Yns= — (s—) Yn’, 
(7.30) M_ Wins = — (8+ 8) Yims’ 

Mg Yims = h-Yms- 
Now using (6.1), (7.28) and (7.30), we obtain 

fyiV,.F Siet8S) SY Sha ee 
(7.31) m=0 s=0 k=—s 
x [(s + 4) HE, Ys’ + (8 — k) 0 yet! — kh #yyh, J]. 

Using (6.11), we may write (7.31) as 
(7.32 fn VF =e > YD NS GRD ve 


where § is a linear operator which will be given explicitly in (7.36). 


1 hout loss of generality we convenie y p 
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D. The Differential Equations for the A;,* 


Combining (7.12), (7.21), and (7.32), we obtain for DF in the case of forces 
arising from an electromagnetic field, 


s 


(7.33) BAe) ue > NE bees + Do(Ant) + E(Amt) + §(4nd)] ve 


k=—=s 


where 


hy i 0 =k —k+1 
DolAme) (25+ 1) ea | hAg- 1 541 than oa 


~ Naha (SHA) (SAA) Amista + (+1) (SHR) (SRA) Awtt sa] + 


+ oe [PAmtit es —h Amis + 


(m-+s+—) (s—k) (s—k—1) Az) —(m-+1) (s—R) (s—h 1) Anta sa] + 


(7.34) 


tay [hm s+1 7 hAgeer 


+ (mt s+) (899) Anh 1 — (m+ 4) (st — 22) Ants, all, 


CA a Pa. Ls =) (s Ris ee AV AS 4 eh Ae is4] + 


(7.35) + 6 (m+s+—) (s+) (s+k—1) A a than tan| = 


and 
(7.36) $(4p%) = —i[(s—k) H. Ane +4 (s+ 2) #.A,t 4 RAAZ". 


By (6.8) and (7.33), we find that the Maxwell-Boltzmann equation for Max- 
wellian molecules is satisfied by the series representation (6.1) if 


(7.37) m=0 s=0 k=—s 


is an identity in e. To obtain the differential equations for the A,,*, we equate 
coefficients of like terms in (7.37). This gives 


k 
(7.38) an Dae OA) OO — Ane: 


The special case of (7.38) for which m=s=k=0O is the equation of con- 
tinuity. The cases containing the local time derivatives of AgimAoviald Ae: 
comprise the equations of motion. Proper combinations of these with the equa- 
tion for the case m = 1, s = k=O give the equations of thermal and of total energy. 
These equations are also consequences of the exact Maxwell-Boltzmann equation. 
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What is interesting here is that the exact form, rather than an approximation 
to these equations, follows from the present theory, based on the linearized collision 


integral. 
VIII. Recurrence Relations for the Eigenfunctions of J 


From the generating function (2.21) for the Sonine polynomials Sip (0), CCD 
are readily derived by elementary manipulations the following recurrence rela- 


tions: 
(8.1) a(x) = Si (x) — Sa (x), 
(8.2) x S™ (x) = — (m+) S™*1 (x) + (n+ 2m +1) Sh (x) — (m+ m) SE" (x), 
(8.3) x Sw (x) = (m + m) Sts (x) — (m+ 1) Seta (4), 
d m—1 
(8.4) ax. —1(%) = — oy (x) ) 
and 
(8.5) x2. Si (x) = m Sz (x) — (nw + m) Se (2). 
In each of these equations x may readily be replaced by hc?, dx by 2hc dc. 
In order to establish (7.7), (7.8), and (7.9), which give (c; icy) yh, and Cams - 
as linear combinations of the Weiss we shall use the following known recurrence | 


relations for the associated Legendre polynomials [13, pp. 107—108, eqs. (3.48), 
(3.49), and (3.50)]: 


(8.6) (2m +1) w P(t) = (mw — m +1) Pra (u) + (w+ m) Pia (u), 
(8.7) (2m 1) Ps Pe) ce Pe (eee ae 
(2m +1) 1 — py? Bi (u) = (m + m) (n + m —1) Pa" (u) — 

— (nm +1) (m — m + 2) PMG (y). 


(8.8) 


Expressing £;=c,-+7C, in spherical coordinates and using (6.2), we have 
(8.9) yk, =csin de? y* ,—T(s —k +1) si4 (Ac?) co*4 sin & PF (cosd) th? 


We substitute into (8.9) from (8. 7) for the quantity sin @ P* (cos 9), obtaining 
a term in P*,,(u) and a term in P* 1(u). In the former term we substitute for 


Se 


Ss'-3 (he) using (8.1) and in the latter we substitute for c?.S”,, (ic?) using (8.3). 
The result is : 


k s—k+ (RF) @ 5 Ss 
ae yi “15 z ——* fhe PAH (u) [S@,q (he?) — SM} (he2)] — 
c~} P&S} (a) [Gnas hays (hc?) = (m1) Sy (he*)}}. 


Applying to (8.10) the definition of yf,, in (6.2), we obtain equation (7.7). Taking 
the complex conjugate of (7.7) using (6.13), and replacing k by — k, we have (7.8). 


We may express c3;=€, in spherical coordinates and make use of (6.2) to 
write 


(8.11) ) Sy. cos Py), = Ts — hE) S34 (hc?) c*7 cos 0 PF (cos 8) ef*?. 
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We substitute from (8.6) into (8.11) for cos # P¥(cos 9), obtaining a term in 
P*,(w) and a term in P*,(u). In the former we substitute for S34 (2c?) using 
(8.1) and in the latter we substitute for c?S%"., (ic) using (8.3). The result is 


Here ab 5 m m— 
(ean) Y= “Gs eayh — {hls — B+) c*PA s(n) [S75 (et) — SP (he®)] + 


+ (s+ hk) co P*, (uw) [(m+s+ 4) S®, (he?) — (m+1) Soe ie ye. 
Applying (6.2) to (8.12), we obtain (7.9). 


The following known relations [/8, pp. 483 —484] will now be used to establish 
equations (7.17), (7.18), and (7.19), which give his (ena. 


0&; ms} ° 

es eho P# (4) R 

Peay, eee 
: = eil(k+l) p f[dk Ss R41 dR (s+1) fi 
Base ae — oR [Zot RIP 

and 

ON tk pk 
(8.14) °° see DRE MING es ee) : 

Scam te Des c RIP a(W) (Ser i llperaseh saa RIP}, 


where &, is defined in (7.5) and R is any differentiable function of c. Taking for 
FR the function 


(8.15) Risé CEna es s+y(hc?), 
we may write 

d — Ss —_ ,—he? ~s—1 a m -~ 22cm 2 | 
(8.16) aa Tag 5 Re Ca iC lc a Seay (hc?) — 2h c? Sst 4 (hc?) 


and 


d R | (s+1) R 


ms | C ms 
(8.17) 


aoe he sot le= tn Sig (2c?) — 2he? Sy (he?) + (28 + 1) Ssiah | : 


Substituting for hc? S?'.4 (hc?) from (8.2) and for gue 


re s+4(4¢7) from (8.5), we 


may write (8.16) and (8.17), respectively, as 


(8.18) ae Rus = RygHe 6 \2(m-+ A) St (hc?) — (2+ 28+ 3).S5)9 (he?) | 
and 


(8.19) Rpt SEY R,, = 2(m $1) oP OSH (he?) — SU, (ho). 


Applying (8.3), with » replaced by s+ 3 and x by hc?, and (8.15), we obtain 
from (8.18) 


(8.20) =~ Rhys = Rs =— 2h Dates 1: 
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Using (8.1) and (8.15), we obtain from (8.19), 

d (ssP 4), = 1)R a 
(8.21) qe Rmst cars Rys = 2(m +1) Rng, s—1 


From (8.13), (8.20), and (8.21), we obtain 


3 Ren at 7 pilktle 
tk pot 
E, [e Ey (4) Kn Al 25 +14 x 


x le h Ee+i Re (u) re. (m y 1) Y Pon Ean ges (0) - 


Multiplying (8.22) through by I’(s—k-+-1) and applying (6.2), we obtain (7.18). 
Taking the complex conjugate of (7.18), making use of (6.13) and (7.5), and re- 
placing k by —, we obtain (7.17). 

From (7.5), (8.14), (8.20), and (8.21), we obtain 


(8.22) @ 


es 
AE. [e a Bee (1) Ne Ail 


ery [—A(s—k+1) Ry iRise) + (m+ 1) (s+) Ripe 
Multiplying (8.23) through by J(s—k-+ 1) and applying (6.2), we obtain (7.19). _ 
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Naberungsverfahren hoberer Ordnung 


fiir Gleichungen in Banach-Rdumen 


Ly GOLLATZ 


Herrm ALWIN WALTHER zum 60. Geburtstage am 6. Mai 1958 gewidmet 


1. Einfiihrung 


Bei den Iterationsverfahren sind die ,,geometrische Konvergenz’ bei Ver- 
fahren 1. Ordnung und die ,,exponentielle Konvergenz“ beim Newtonschen Ver- 
fahren bekannt. Diese exponentielle Konvergenz gilt auch fiir andere Iterations- 
verfahren héherer Ordnung und im folgenden soll eine allgemeine Theorie fiir 
Konvergenz dieser Verfahren, Existenz von Lésungen, und fiir Fehlerabschatzun- ~ 
gen von Naherungswerten aufgestellt werden. Damit die Ergebnisse auf ver- 
schiedene Bereiche der praktischen Analysis anwendbar sind, werden die Methoden 
funktionalanalytisch durchgefiihrt und die zu lésende Gleichung 7 #—=# in einem 
Banachschen Raum mit # als Nullelement betrachtet; w ist gesucht, der im 
allgemeinen nichtlineare Operator T sei stetig und zwei- oder mehrmals differen- 
zierbar. Die Lésung uw wiinscht man sich als Grenzelement einer Folge /,, von 
Elementen, derart, daB T/, gegen # konvergiert; es soll ¢,=||Tf,|| gegen Null 
gehen. Von dem Iterationsverfahren werden zwei sehr naheliegende Annahmen 
getroffen: 

Es soll die Anderung 6,=/,,;—/, mit ¢, klein werden und es soll bei jedem 
Schritt 4, abnehmen nach einem Potenzgesetz, dessen Exponent eine fiir das 
Verfahren charakteristische Konstante, die Ordnung & ist [genauere Formulierung 
fiir beide Annahmen in (2.7) und (2.8). 

Nach einigen Vorbereitungen wird in Nr. 4 der allgemeine Existenzsatz mit 
Fehlerabschatzung bewiesen; Nr.5 und 6 geben Hinweise fiir die praktische 
Durchfithrung der Rechnung und die einfachen Beispiele in Nr. 7 sollen zeigen, 
da8 man scharfe Abschatzungen aufstellen kann und daB es Falle gibt, bei denen 
sich die Mitfthrung hdherer Ableitungen empfiehlt, z.B. bei Matrizen-Eigen- 
wertaufgaben, bei denen die zweite Ableitung des geeignet gewahlten Operators 
konstant ist und die dritte Ableitung identisch verschwindet. 


2. Bezeichnungen und Hilfsmittel 
Vorgelegt sei in einem Banachschen Raum B eine Gleichung der Form 


(2.1) Tu=%. 


Dabei sei T ein in einem vollstandigen konvexen Bereich F des Raumes B defi- 
nierter stetiger, linearer oder nichtlinearer Operator und # das Nullelement in B. 
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Es wird nach ,,Lésungen“‘ oder Elementen wu aus F gefragt, die die Gleichung (2.1) 
erfiillen. Nun sei T ein oder mehrere Male im Frechetschen Sinne differenzierbar : 
die Ableitungen werden durch 7’, T’’, ... bezeichnet. Fiir Gleichung (2.1) werde 
ein Iterationsverfahren benutzt, welches eine Folge von Elementen fo, f,, ... 
liefert. Es bezeichne 7,” die k-te Ableitung von T an der Stelle f,: 


(2.2) lie pa Dis fie ne (7) : 


Es wird zugelassen, daB bei dem Iterationsverfahren zur Berechnung von f,,;, 
die GréBen f,,, gewisse Ableitungen 7," und die Inverse J,’ ~* des linearen Opera- 
tors J, benutzt werden, und zwar werden Iterationsverfahren der Form unter- 
sucht 


(2.3) p= fatr— = Oi a Lil Ly 5 eas) 


mit gegebenem Ausdruck ®. (3.2), (3.5), (3.6) geben Beispiele solcher Ausdriicke. 

Im folgenden wird der Mittelwertsatz benutzt. Es seien g und g+k zwei 
Elemente von J’, und das ,,Geradenstiick“ /" der Elemente g+¢k fiir O0<t<1 
gehort auch zu I’. Ist T langs £’ mindestens (p+-1) mal differenzierbar und gibt 
es eine Zahl M,,, mit 


a parte lls Maun fir fer, 
so gilt? 

|Tet+ aT @et sl) k k++ TM @)k...k-TE+D)| 
(2.5) irl, 


= My, ||A|/?™ 


(T’(g) k bedeutet, daB der lineare Operator 7’(g) auf k angewendet wird. 
T’'(g)-k- m bedeutet, daB der bilineare Operator 7’’(g) auf die beiden Elemente k 
und m angewendet wird). 


Es werden die Abkiirzungen verwendet: 
(2.6) b= (Tr "lL, = [ITIL 
Es zeigt sich, daB bei einer Reihe bekannter Verfahren Abschatzungen folgender 
Gestalt gelten: 
(2.7) tn41 SG (by, t,) th 
(2.8) IL onll SA Ons tn) tn 
wobei G und H angebbare in ihren (nichtnegativen) Argumenten monoton nicht 


fallende Funktionen sind und k eine positive Konstante =1, die sogenannte 
, Ordnung“ des betreffenden Verfahrens ist. 


3. Einige bekannte Iterationsverfahren 
Im folgenden sei 


(3-1) grlIT? OSM, fir feF, 


alle vorkommenden Argumente f,,, 7, + 6,,/,+¢,,--. mogen in F liegen. 


1 Vgl. etwa L. V. Kantorovicu, Functional Analysis and Applied Mathematics, 
engl. Ubersetzung, herausgeb. v. G. E. ForsyTHE, 1948, Ss G2. 
Arch. Rational Mech. Anal., Vol. 2 a) 


68 lik (Copies Naw75 


I. Beim gewdhnlichen Newtonschen Verfahren lautet die Iterationsvorschrift: 


(3.2) gn=hu-bh=—-h*t, oder 1,4+%,6,=0. 


Ww 


Hier wird T als mindestens zweimal differenzierbar und die Existenz einer Zahl 
M, nach (3.1) vorausgesetzt. Nach (2.5), (2.6) folgt dann 


(3.3) | Sal] SUZ WM Zal SS bn bn 
Weiter ist nach (3.2), (2.5) 
(3.4) tyt4 7 |Zna <1 iL az bh é,|| < M, || 6, ||? = M, Di, ae 


II. Es gibt zahlreiche Verscharfungen des Newtonschen Verfahrens. Man 
kann z.B. ein Element ¢, wie beim gewdhnlichen Newtonschen Verfahren und 
dann die Korrektur 6, mit Hilfe der 2. Ableitung berechnen: 


(3.5) T+ IT en) 
(3.6) T+ Tn bn + 2 In bn bn = 0. 
Nach (2.5) gilt 


|Z L 1 On, Z Ee 6, 6n|| S M||6, 


ww 


, 
also mit (3.6) 
Fria 3 Ly On (On = Gn) oe Dx (Op Gea ee | 


oder 
(3.7) (IZ < M,(||6,,|| a aD Il, =, | ai Ms || 6,,||°. 
Nun gilt wie bei (3.3) 
Il Sn || Sb, bn 
Aus 
Tr (6, — oe) Hele Cn = 0 
folgt 
|6, —¢,|| S68 MN, eee 
und damit 
(3.8) Onl S16, — Call + [Sul] SG, + Me 23 t,) t,. 


SchlieBlich ergibt sich aus (3.7): 
(3.9) tn41 <= Mz b;, (2 + M, di:t,,) t, + MG|I 6, ||. 


III. Eine vereinfachte Variante erhalt man, wenn man in (3.6) die zweite 
Ableitung nicht bei jedem Schritt neu berechnet, sondern fest wahlt, also nach 


(3.10) T+ Tron = 0, 
(3.41) Dat tn One to enen 0 
rechnet. Wie bei (3.4), (3.7) gilt dann 


(Za—- f= Tr 8n|| S Mg|| 6, D lem ll Ont, 
tn41 S23 (|Zo || ont, + MG|| 6, |]? 


[Onl] SG, + 3 || Zo'|] 2 t,) ty 
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4. Der allgemeine Satz 


Satz. Iiir das Iterationsverfahren (2.3) zur Lésung der Gleichung (2.1) mégen 
zwei fiir nichtnegative b,,, t, monoton nicht fallende Funktionen G und H mut (2.7), 
(2.8) extstieren. Der erste Iterationsschritt von einem gewdhlten Element f, zu f, set 
ausfiihrbar und by=||Tp ~1|| bekannt. Es sei ein vollstandiger konvexer Bereich F 
und eine Zahl y>by mit folgenden Eigenschaften vorhanden: 


1. In F ist T definiert, $\|T’|| <M, und dort existieren Ableitungen so hoher 
Ordnung von T, wie sie in (2.3) auftreten. 


2. Mit 
(4.1) o=Hy,b), »=Gly,), B=vt*, g=2Mot, 
gilt 
(4.2) 0O= oy = 1p 
und 

y? 

(4:3) ert (y— 9) bo =: 

3. Es gehéren fy, f, und alle Elemente s der Kugel S 
(4.4) Ils— Alls ot 
zu EF, 


Dann gilt: Die Iteration (2.3) ist unbeschrankt ausfiihrbar, es existieren alle 
Inversen T,~1 und sind beschrinkt, b,<y; die Gleichung (2.1) besitzt in F eine 


—— 


Lésung u, die sogar in der Kugel S liegt, und es besteht die Fehlerabschatzung 


to Ibe 
(4.5) Ile — fall S- re ena.) 
mut 
(4.6) Bo Ree re 0B =e 


Beweis durch vollstandige Induktion. 


Fiir ein bestimmtes natiirliches ~ (und alle kleineren natitirlichen Werte von 1) 
gelte: 


a) es ist 
(4.7) Pea Brie 


b ths ele 


C\af eq CXIStierL, 


e 


f 
(4.8) b, Sy — (vy — 4) B, Sy- 


aes 
Re 
d) man bilde y,=20, M, ct,; es gilt p,S9<1, 
Wie, Wegt in a, 
) 


es ist 


Diese 6 Aussagen sind fiir »=0 nach den getroffenen Annahmen erfiillt. Es 
wird gezeigt, daB die 6 Aussagen auch fiir +1 gelten. 
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a) T,,, existiert nach c) und es gilt nach (2.7) 


a+ 


(4.9) Le GO ent, = Gly, to) Povi svi i B= Blo pa ere =) ByiiSty 


Das ist Aussage a) fiir m+ 1. 
b) Es ist 
|Z a (Ce “3 Tr+1) He b, ||Tn is, T,+1|| S26, M, || nl] = Yn 


(4.10) < 2b, M, Hy, to) (20 VEO lL, Oo, 


wobei die Bezeichnung g; eingefiihrt wurde. 


Fiir den linearen Operator W,—T, ~* (TZ, — Ty41) existiert wegen || W,||< 14 
nach einem bekannten Satz? der inverse Operator zu E —W,=Z,,; es ist Zp 
Lod wt 2, hat dain awe ui T,.1 einen inversen Operator Z,, 1, ~\undes gilt 

bale 


* ? 


Wom Pn 
also 


c) Mit 7,4" existiert auch f,,9. 
d) Nach Definition, (4.9) und (4.2) ist 


k—1 
os 2041 Me Ob ay = Sf a ove, ty 


— Pn Tae tee aE Pn De Pn y iio ee ps Pn: 


1— Py eee | an Vink ened ae 
Die @, bilden also eine monoton nicht wachsende Folge. 


e) Fir 7,4 gilt nach (2.8), (4.1), (4.7) 


n+1 n+1 nord " oly B 
Ivo Alls 216s Lots X ot B,=oto[B } BPR oe BERT BS | ee 1— pe: 


Nach (4.4) gehért also f/,,,. zur Kugel S und damit zu F. 
1) Weeen 0= 9, =@y =< 11st 


jae Pn 
= 7 hee 1— M 
(denn ae sc get fears 2) < 1 Pot Pn(Po= Pn) __ | Pn 
>on (1—¢,,) (1 — qo) crs (1— 9) (1— @) 05 ss 


Damit wird 


by, n 
Ont S =, (44 ae | = 6, (4 + yqB,) 
1— 


> Pn 
mit 
= 2M, oly ee Po 1 
tO 1—P bo” 


2Vgl. etwa L. A. LyusterNnIk-W, I. SOBOLEW, Ele 


mente der Funktionalanalysis, 
Berlin 11955; 53401, Satz 2: abe 
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Es folgt nach (4.8), (4.3) 


ry —(Y — 05) Bi tg Dy) =» —(y —b,) B, g* 
mit 


= 2 2 2 
ey By ee ie Dee peape 
7 vq yi 06 on y— bo (vy — Bo) bo 1— =p26 
oder 


TA =y— VY  %) Bay Ged 


Somit gelten die Aussagen a) bis f) fiir alle natiirlichen Werte von n. Fiir m=>0 
gilt weiter 


lFntnta — Jal 
(4.12) nm nm n n+1 n+m—1 
= 2 ||6|| Sot: B,, (4 suo age: sake ‘aban ea ’) = oe ° 


Die Schranke geht fiir 7— oo gegen Null; die (in der Kugel S liegenden) Ele- 
mente /, konvergieren daher wegen der Vollstandigkeit von S gegen ein in S 
liegendes Element w. LaBt man in (4.12) m— co gehen, so ergibt sich die behaup- 
tete Fehlerabschatzung (4.5). Da B,—>0 fiir noo strebt, geht ¢, nach (4.7) 
fiir »—> co gegen Null; es ist daher uw wegen der Stetigkeit des Operators T eine 
Lésung von (2.1). 


Zusatz. Es gibt Falle (z.B. Beispiel I in Nr. 7), in denen man wei8, daB 7’ 
in ganz F eine Inverse besitzt, die in F gleichmaBig beschrankt ist; dann kann 
man diese Schranke als y wahlen. Die Voraussetzung (4.3) und der Teil f des 
Beweises kénnen dann fortfallen. 


5. Rekursive Rechnung 

Es kann eintreten, daB das Ausgangselement /) zwar zu einer gegen eine 
Lésung w konvergenten Folge /, fiihrt, daB aber der Ausgangswert ¢) noch zu 
groB fiir eine brauchbare Abschatzung nach (4.5) ist, daB z.B. der Wert von B 
zu nahe an 1 (oder gar oberhalb 1) liegt oder die Ungleichung (4.3) nicht erfiillt 
ist. In solchen Fallen [insbesondere wenn die nach (5.1) fiir 4; erhaltene Schranke 
wesentlich kleiner ausfallt als ¢)] kann es zweckmaBig sein (vgl. Beispiel II in 
Nr. 7), einige Betrage ||6,|| direkt abzuschatzen und erst dann die allgemeine 
Formel (4.5) zu benutzen. 


Aus (4.11), (2.7), (2.8) gewinnt man, falls gy <1 ist: 


b 
< : 
(5.1) i= 1—2b)M,|| dol] ’ 


4=G(b,t) 8; ||6,||S 2,4) 4. 


Man kann dann b,—0¢, t,=%, 6,=069 setzen und im Satz von Nr. 4 statt von 
by, ty, Oo von bg, to, 09 ausgehen und erhalt dann statt o, v, 8, B, natiirlich neu 
zu berechnende GréBen o*, »*, B*, BX; es wird dann z.B. 


at if Be 


(5.2) I~ — All Sell + 1— pe? 


Notfalls kann man diesen Schritt mehrmals wiederholen und einige weitere || 6,]|| 
abspalten. 
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6. Zusammenstellung 

Allgemeiner Rechnungsgang. Es liege eine Gleichung der Form (2.1) vor und 
es wird ein Iterationsverfahren der Form (2.3) benutzt, etwa eines der Ver- 
fahren (3.2) oder (3.5), (3.6) oder (3.10), (3.14). Zur Anwendung des Satzes aus 
Nr. 4 sind folgende Schritte durchzufihren: 

4. Wahl einer méglichst guten Ausgangsnaherung f,, Aufstellung von f, durch 
Ausfithrung eines Iterationsschrittes und von 69 =/f,— fo: 

Die Anwendung der Theorie ist nur sinnvoll, wenn 69 Sehr klein“ ausfallt, 
vel. hierzu Punkt 5. 

2. Aufstellung der fiir das gewahlte Verfahren geltenden Beziehungen (2.7), 
(2.8). Es ist: 


Gewohnliches Newtonsches Verbessertes Newtonsches | Vereinfachte Variante 
Verfahren (3.2) Verfahren (3.5) (3.6) vom vorigen, (3.10) (3.11) 
G (ns ty) My bi M2b8(24.M,b24,)-+MyH® | 4||74'|| 024M, H2 
(On tn) bn by tM, b3 ty Gye eos 
k 2 3 | 2 


Man braucht hier die 2. Ableitung T’’ [und bei (3.5), (3.6) auch die 3. Ab- 
leitung|; man muB einen konvexen Bereich F wahlen, der f/, und /, enthalt und 
vollstandig 1m Sinne der verwendeten Metrik des Banach-Raumes ist, und die 
Maximalbetrage von T” und T’” in F abschatzen, da sich aus diesen M, und VM, 
nach (3.1) berechnen. 


3. Aufstellung des inversen Operators 7j~1! und einer Schranke fiir seinen 
Betrag i 
bo 2 ||Zo II. 


Ferner braucht man eine Schranke ¢) = ||J||. 


4. Wenn die Voraussetzung des Zusatzes in Nr. 4 erfiillt ist, kennt man Y, 
berechnet o,, 8, y) nach (4.1), und wenn (4.2) gilt, liefert (4.4) eine Fehler- 
schranke fiir f,. Wenn der Zusatz aber nicht anwendbar ist, muB man y wahlen ; 
es muB y> by sein. Dann findet man o, y, B, gy nach (4.1). 

5. Priifung der Ungleichungen (4.2) und (4.3). Sind beide ,,reichlich“ erfiillt, 
so benutzt man (4.4) als Fehlerschranke fiir f,. Sind die Ungleichungen nur 
, schwach“ oder nicht erfiillt, so versuche man die rekursive Rechnung nach Nr. 5. 
Versagt auch diese, so ist wohl die Wahl eines besseren Ausgangswertes /) zu 
empfehlen. 

7. Beispiele 
I. Als ganz einfaches Beispiel sei 2 zu berechnen. 
Vorgelegt sei die Gleichung 


Du =g(u) =w—2=0. 


Das verbesserte Newtonsche Verfahren (3.5), (3.6) lautet hier: 


(74) tein =P,) mit Yu) = — ae lew) + Fe" (SS), 
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1. Es werde fy =1,4 gewahlt, dann wird 
Ty = 8 (fo) = — 0,04; To = g'(fo) = 2fo= 2,8;  g'"(u) = 2. 


Nach (7.1) wird dann /, = 1,414212828.... 


2. Hier ist T”=2, T’’=0, also M,=1, M,=0. Als Bereich F werde etwa 
das Intervall 1,4<4u<1,5 gewahlt. 


350 5) 1 5 
4Ye Th oy bia he fA 397 tes ty = 0,04. 


4. Der Zusatz von Nr. 4 ist hier anwendbar, es ist 


1 1 
he Ave a 
ew) | = 2, fiir uwCF. 
Weiter wird nach (41) »=G(y, 4) =0,0326, o=H(y, 4) =0,360, p=r = 
0,0000522, gy=0,0103 und schlieBlich nach (4.4) 


|~ — f,| <0,000000752. 
Der tatsachliche Fehler von f, ist fast so groB wie diese Schranke und betragt: 
u — f, = 0,000000 734.... 
II. Vorgelegt sei eine Matrizen-Eigenwertaufgabe 


Ax =AB% 


(2) CL) 
=i 4 = 2 Gl 


Bei einem so einfachen Beispiel mit zweireihigen Matrizen kann man natiirlich 
die Eigenwerte bequem direkt ausrechnen, es wird A=2(9+)/21), aber man 
kann hier die allgemeine Methode gut illustrieren. Man verwendet als Banach- 


mit den Matrizen 


% 
Raum die Kombination u = ( ; aus Eigenvektor und Eigenwert. Zur Normierung 


des Eigenvektors nimmt man ein Funktional s(x), z.B. ee Komponente von x, 
die bei dem gesuchten Eigenvektor nicht verschwindet, hier z.B. die zweite 
Komponente %,= s(x). Dann liegt die Gleichung vor 


- = x\ Ax—iABx 4 
var() (ici -s 


Die erste Ableitung ist ein linearer Operator T7’(u), der auf ein Element ‘| 


angewendet lautet: k Ak—wBx—ABk 
T’(u = j 
ie ( | ae, 


Die zweite Ableitung ist ein konstanter bilinearer Operator T’’(u), der auf Ele- 


k 
mente und () angewendet ergibt 
lu y 


real) (n°) 
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Die dritte Ableitung ist ein Nulloperator, der auf drei beliebige Elemente 
() (’) (") angewendet stets # liefert. 


a — 0,56 
4. Als Ausgangsnaherung werde /)= 1 , also A>=0,74 benutzt. Ein- 
setzen ergibt 0,74 
ees 27h Ae 
Th=| —0,0088], A—0,74B= a a 
0) 
Ray 2,74 Ray + 4,52 Rey — 2,56 4 
Te | Rey ) = | 0,48 Ray + 0,26 Re, — 2 : 
ie Rs) 


Zu Ty ist die Inverse zu bilden: 
55) 264, 103392 
AVAS Ty ea 0 OUPMIaS 
12 —=68,5 —0,43 


Es wird 
; 20 0,001 746725 
e=()=-ho m= gee (  of=| 0 
Mo - NeeaaS — 0,003 755.458 


Das gewohnliche Newtonsche Verfahren wiirde als nachste Naherung ergeben 


== 0;5582535 
fotlo= 1 ? 
0,736 2445 


wahrend das verbesserte Newtonsche Verfahren nach (3.5), (3.6) 


TREY EE — 0,000 004 297 
Oo a Co ar To : 0 = lg aE 0) 
— 0,000007161 


mit nur geringfiigig vermehrtem Rechenaufwand eine viel bessere Naherung 
liefert : 


— 0,558257572 
h=ft+ One 1 
+ 0,736 237384 


k 
2. und 3. Als Norm eines Blementes( werde einfach Max (| &q)| -| Ray, ¢| |) 
lu 


fiir c=1 verwendet und als Betrag von T’ baw. T’~1 jeweils das Maximum der 


Summen aus den Betraégen der in einer Zeile stehenden Zahlen (Maximum fiir 
die verschiedenen Zeilen). So wird 


f= ||Tol| = 0,0144; by = || Ty -1|| = 23+ O4 + 63,92 


114,5 = 1,58. 
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3 ||T’|| wird das Maximum der Zeilen-Betragssummen bei B: 
Veer |g teva Se Mes 


4. Der Zusatz von Nr. 4 ist hier nicht leicht anwendbar. Versuche mit 
y=1,8 und 2 zeigen, da8 dann zwar die Ungleichung (4.2), aber nicht (4.3) 
erfiillt ist. Da aber ¢, und 6, bereits ,,klein“ ausfallen, liegt es nahe, die rekursive 
Rechnung nach Nr. 5 anzuwenden. Mit || 69|| = 0,003 763 berechnet man nach (5.1) 


Phe AtOA hoe Lied 626-40", 1.4; || = 2,68-1057% = || 63||- 


Diese GréBen liefern dann o* = 2,000004, v*= 288,0006, B*= 7,65 - 10-2; = 
3,91-10~®; die Ungleichungen (4.2), (4.3) sind fiir die GréBen mit Stern ,,reichlich‘‘ 
erfiillt (die Existenz einer Lésung u bei f, ist damit gesichert) und (5.2) ergibt 


Il“ — All| S 2,7- 107-7. 
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On a Partial Differential Equation 
with Two Singular Lines 


E. T. Copson & A. ERDELYI 


§ 1. Introduction 
In an earlier paper! one of us investigated the partial differential equation 


ie eU , 2a aU @U , 26 aU 
(1-4) one | 4 Ox ay? | y Oy ~ 


This partial differential equation has two singular lines, the axes of coordinates, 
and its investigation might be expected to throw some light on the properties 
of partial differential equations with intersecting singular lines. 

In I, those solutions of (4.1) were investigated which are continuously dif- 
ferentiable'* in the closed quadrant x20, yO, and twice continuously dif- 
ferentiable in the open quadrant x>0, y>0. It turns out that such solutions 
must satisfy the supplementary conditions 


ou 
Ox 


==) “for ) x= 0,0 yes03 S =0 OPER 0 ay — 03 


(1.2) 


If the axes were not singular lines, one would expect to be able to prescribe a 
further boundary condition along one of the two axes, giving, for instance, one 
of the two functions 


(13) U(x,0) =F(x) for x20; U@,y)=G(y) for yo. 


Such is not the case here. One of the most noteworthy results of I states that 
in the event 6>«a, the additional boundary condition should be prescribed on 


the axis of x. This is so because /*(x) and G(y) are connected by the very simple 
functional equation 
1 


(1.4) G(y) = nena x f On — a Fyn dt 
0 


from which it follows that G(y) must be more regular than merely continuously 
differentiable. 


* Corson, E. T.: On a singular boundary value problem for an equation of hyper- 
bolic type. Archive for Rat. Mech. and Analysis 1, 349—356 (1958). This paper will 
be referred to as I, so that I § 6 and I (2.1) refer respectively to section 6 and equation 
(2.1) of that paper. 

ee The condition of continuous differentiability on the boundary could be relaxed 
See in this connection J. B. Diaz & G. S. S. Luprorp, On the Euler-Poisson-Darboux 
equation, integral operators, and the method of descent. Proc. of the Conference on 
Differential Equations held at the University of Maryland, March Wey ise I GOSS 
pp. 73—89, in particular p. 81. eh . 
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In I the partial differential equation (1.1) was investigated by RIEMANN’S 
method, and the explicit form of the solution was obtained from the known? 
Riemann-Green function of (1.1). The relation (1.4) was then deduced by means 
of some transformations involving special functions. 

In the present paper we propose two other methods for investigating the 
partial differential equation (1.1): integration of fractional order, and the Mellin 
transformation. Each of these methods will enable us to recover (1.4) and to 
obtain some other results, complementing those obtained in I. 

It has been shown in I (3.5) that (4.4) may be written in the form 


EGOS res B 7P-a7p—5 F(Ve)] 
de z 
ees ain x48 : 


where J°~* denotes the Riemann-Liouville operator of integration of fractional 
order 6 —a, and Ij_, ,_, denotes KoBer’s operator®. This form, in conjunction 
with a known result* on generalized axially symmetric potentials suggests the 
application of integration of fractional order to (4.1). It turns out that an 
appropriate integration of order «’—a with respect to x? changes (1.1) into a 
similar equation with «’ in place of «, and integration of order 6’—f with respect 
to y?, into a similar equation with f’ in place of 8. Thus, integration of fractional 
order may be used to correlate two partial differential equations of the form (1.1) 
differing only in the values of the constants « and f. 

We shall use in particular two such correspondences. If B>a, we may 


correlate (1.1) with 
(1.6) eU ! 20 LU eU ! 2a CU) 

On" OG Oy? y oy 
The solutions of the latter equation are symmetric functions of x and y, and this 
remark will at once lead to (1.4). Secondly, (1.1) may be correlated with the 
wave equation (« =f =0), leading to a “general solution” of (1.1) involving an 
“arbitrary” function of a single variable. 

Our second method involves transforming (1.1) to polar coordinates 7, # and 
applying the Mellin transformation with respect to 7. The resulting ordinary 
differential equation can be transformed to the hypergeometric equation, thus 
leading to an explicit solution of our boundary value problem. In this method 
(1.4) appears as a consequence of Gauss’ formula for the sum of the hypergeometric 
Benes at 21. 


§ 2. The classes ©, 
For the sake of brevity, let us denote by ©, the class of functions U(x, y) 
satisfying the following conditions: 
(2.1) U is continuously differentiable in the closed quadrant x20, y=0; 
(2.2) U is twice continuously differentiable in the open quadrant *>0,.y>0; 


2 Copson, E. T.: On the Riemann-Green function. Archive for Rat. Mech. and 
Analysis 1, 325—348 (1958). In particular § 5. 

3 Koper, H.: On fractional integrals and derivatives. Quart. J. of Math. (Oxford) 
11, 193—211 (1940). ; 

4 Erpéxyi, A.: Singularities of generalized axially symmetric potentials. Comm. 
Pure and Appl. Math. 9, 403—414 (1956). 
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(2.3) U satisfies the partial differential equation (1.1) in the open quadrant 
x>0, y>0; and 
(2.4) U satisfies (1.2). 


It may be noted that the first condition (1.2) is automatically satisfied if 
a> 0, the second condition (1.2) if B>0. (2.4) has been included in the definition 
of ©, because we should like to consider this class of functions when «20 and 
f=0; and in the cases « =0 or f =0 the additional hypothesis (2.4) is appropriate. 

The boundary value problem solved in I consists in finding a function U 
in ©,, which satisfies U(x, 0) =F (x), and it was shown in I that such a function, 
at least when }<a<f, exists provided F(x) is continuously differentiable for 
x= 0 and satisfies f'’(0) =0. 

It will now be shown that the application of the operator J; _; ,,_,, of integration 
of fractional order «’—« with respect to x? transforms a function of Cg into a 
function of C,,,. 

We assume that OSa<a’, and that U(x, y) is in ©; we set 
U(x, ¥) 

Ox? 


OU (x,y 
Up ae, 


and consider the function 


U'(x, 9) = f Ba — B94 UE, 9) dé. 


Uy9(%,Y) = 


> 


Here differentiation under the integral sign is permissible (even under the 
less restrictive assumption —}<a<a’). Clearly U’ satisfies conditions (2.1), 
(2.2), (2.4). Moreover, the following identity may be verified by an elementary 
computation: 

iG fae £2)¢ a4 {Us (E x, y) = = Ce (E x, y) OU (Ex, y) Kak 2a’ OU (Ex, a 

S 


Ox? % Ox 


AD Ed , 
aera [ead 62 Seaaea y)]. 
Now® 
a 
1 d 2a 2) a’ —a 
pies yo Caer en 
0 
since U,,(0, y) =0 by (1.2), and we thus see that 
il 


ie 20) COU d - ae F 
Ox? ' 4 Ox = ih dé Lore Giace ye 2% Uy (E x, y) | dé + 
0 
1 
ae het = £2) 2-1 | Ug (E i y) ot = U,9(é a y) de 
0 


ele =P) Lae oy pe 


ean C ie ees 
Gye Neva ovis 
* This holds for «> — 4 since Uy (Ex, y) =O(&%) as €>0. 
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thus, U’ satisfies the partial differential equation appropriate to €,, g- Lhere are 
similar results for integration with respect to y? and double integration with 
respect to both x? and y?. Summarizing, we obtain: 
If OSa<a’, OS B<f’, and U(x, y) is in Cg, then 
il 


(2.5) Jeet = 6), Sa Ulery dé 

1s 1n Syg, 

(2.6) iene =) vo UN) ay 

is in Cy, and 

(2.7) Sf BP? (1 — 2)" (1 th 4 (Ew, ny) Edy 
1s in Vy pr. 


The above formulae may be written as fractional integrals. Indeed, with 
KOoBER’s operator of integration of fractional order, 


28) Tat@ =a fer tee ae, 


0 


the expressions (2.5), (2.6), (2.7) are constant multiples of 
jE eG NC Th 46-8 Wis Vek ey he: 3,6’—-6 U(%, ¥) 


respectively. Our result may then be written symbolically 


(2.9) ee Oe Uns OSa<a’, 
(2.10) Te —p Gap < Cape 05 = 6; 
{2.11) Vhs La’ aie au ay p Sap C Care (OS me <n Os Hap 


The symbol © in (2.9) indicates that while we have established that every 
expression of the form (2.5), with U in ©,,, belongs to G,,, it has not been 
proved, and we believe it is not true®, that every function in C,,, can be represented 
in the form (2.5). Similar remarks apply to (2.10) and (2.11). 
It follows from (2.8) by the change ¢ =z? of the variable of integration that 
r( 


(2.12) Ef) = peggy tO Ae 


Applying this formula to the mapping (2.9) of ©,, into ©, , we find for 


Oa ea ie 

, ne I'(a +4) 
(2:43) U"(0, ¥) = Ta’ +4) U(0,y), 
(2.14) Uk G0 R= A Pero x, 0): 


6 For the reason that even if U’ satisfies (2.1) and (2.2), and is of the form (2.5), 
it does not follow that U(x, y) satisfies the regularity conditions (2.1) and (2.2). 
Moreover, it is not true in general that a function satisfying (2.1) and (2.2) can be 
represented in the form (2.5) for all non-negative «<«’. 
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§ 3. The class ©,,, and its connection with other classes 

If «=f, our partial differential equation is (1.6), which is symmetric in % 
and y; and it is natural to conjecture that solutions of (1.6) satisfying both 
conditions (1.2) will be symmetric functions of » and y. Before we verify this 
conjecture, we shall represent explicitly that solution of (1.6) which satisfies the 
conditions 
(Godly Unit(s), xa =O. for 7205) =e. PUL 20 dot. i= O29 220 

oy OX 

This we shall do by transforming (1.6) to the Euler-Poisson-Darboux equation, 
and using the known solution of the latter equation. 

Assuming 0X y<~ in (1.6), we transform this partial differential equation 
to variables 
(3.2) E=+y, yn=2xy, OSnSE, 
and obtain 

eU  e8U , 2« aU 


6-3) GE Gif ek On 


the Euler-Poisson-Darboux equation with the singular “‘initial’” data U =F(/é), 
oU/én =0 for 20, 7=0. The known solution of this problem for «>0 1s? 


(2) P(x) 
=i 
or 
hte 
6.4) UG) sar, pe) Blond 
=i 
where 


wt = x%* + y24- 2%. 


This can also be written in the form 


x+y 
D(a 5 : =) 2 
(3.5) Ul%,9) = Fate (2xy)* ** [ R™ Flo) odo, 
x—yYy 


where 
R* = [w? — (x — y)"] [(«@ +)? @?). 


From this formula we can deduce that a solution of (1.6) and (3.1) which 
vanishes everywhere on the characteristic x = +, vanishes identically. For, we have 
from (3.5)8 

Bee Vane 
Ona) = T@Pa) Qua fo * (4x? — w)*-1 F(@) do. 
0 


* WEINSTEIN, A.: Sur le probléme de Cauch : i : 
; y y pour 1’é t , 
tion des ondes. C. R. Acad. Sci. Paris 234, basa foe SPER ga Ce 


* The formula for U(x, x) given in I § 4 reduces to this when piso 
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Now, this is, apart from some obvious changes of variables, ABEL’s integral 
equation. If U(x, x) vanishes identically, then® so does F(x) =U/(x, 0), and, by 
(3.4), U(x, y) for OS yx. A similar argument with the roles of x and y inter- 
changed shows that U vanishes also for O<*<y. This proves the statement 
when «>0. Ifa =0, (1.6) is the wave equation and 2 U(x, y) =F(x +) +F(«*—y) 
for OS ySx, 2U(x, x) =F(2x) +F(0), and clearly, U(x, x) =0 for all x implies 
F(x) =0 forall x. 

We can now prove that for «20, the functions belonging to C,,, ave symmetric 
functions of x and y. For let U(x, y) be such a function. Then V(x, y) = 
U(x, y) — U(y, x) also belongs to ©,,, V clearly vanishes on the characteristic 
x =, and hence vanishes identically. 

Since U is a symmetric function of x and y, we have the following representa- 
tion. Let U be a function of ©,,, with «= 0, and let U(x, 0) =U(0, x) =F (x). Then 
(3.6) 2U(x,y) =F(a+y)+F([x—y|) if «=0, 
and U 1s given by (3.4), or by (3.5) with the lower limit of integration replaced by 
|x«—y]|, af w>0. 

We may now use the mapping (2.6) in order to establish a connection between 


C,,, and C,, if B>a. Let U be in ©,,, and let V be that function in ©, for which 
Via, O)i—= U (4; :0) 1%). Then 


1 
(3.7) U(x, 9) = Fen aie @ imc — 9 )P-** V (x,y) dn. 
0 


Moreover, V, being in €,,,, is a symmetric function of x and y, V(0, vy) =V(y, 0) 
=F (y), and hence we have (1.4) for G(y) = U(0, y). 

Representing V(x, y) in terms of F(x) by means of (3.4), we also obtain from 
(3.7) a double integral representation of U in terms of I. 


§ 4. A “general solution” of (1.1) 


We first find the general form of a function W in ©),. Such a function satisfies 


the partial differential equation ew ew 


On Oy? 


and hence is of the form 
W(x,y) =f(~ +y) +8(*% —9). 
We know from § 3 that W is a symmetric function of x and y, so that g(x) must 
be an even function of x. Clearly, # and g must be continuously differentiable. 
Setting W(x, 0) =W(0, x) =@(x) x20, we have from (1.2), 
t(™) +e(*) =@(x), F(x%)—8'(%) =O *%20. 
By adjusting an additive constant, we may take 2/(x) =2g(x) =p(x) for x20. 
Since g is an even function, g(x) =q(|x|) for all real x. We thus have 
2W(x,y) =o(x +¥) + (|x — 9) 
for ~=0, y=0. By the definition of €)y, v(x) must be continuously differentiable 
for x=0 and gy’(0) must vanish. We shall now show that in the present case 


9 Compare a similar argument in A. G. Macxktz, The generalized radially sym- 
metric wave equation. Proc. Roy. Soc. A 236, 265—277 (1956), in particular p. 269. 
Arch. Rational Mech. Anal., Vol. 2 6 
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must possess two continuous derivatives. Indeed 
fi ow (x, ow (x, v) 
(e+) =o zs 

is differentiable in the open quadrant x>0, y>O and hence @ is twice con- 
tinuously differentiable except possibly at the origin; and 

elle—9)) ee 
is also differentiable in the open quadrant, and hence wy must also possess a second 
derivative at the origin. 

Let p(x) be twice continuously differentiable for x20, and let y'(0) =0. Then 
y(x+y) +(|x—y|) belongs to Gog, and conversely, every function im Go can be 
so represented. 

Thus, in the case of ©)», the hypotheses (2.1) and (2.2) result in the function 
in question being twice continuously differentiable in the closed quadrant x20, 
y=0, in particular, having twice continuously differentiable boundary values. 
This is the more noteworthy as for B2a«> 4% it has been proved in I §§ 4 and 6 
that the boundary value U(x, 0) =F (x) may be any continuously differentiable 
function whose derivative vanishes at x =0. 

Having a representation of functions of ©)), the mapping (2.7) may be used ~ 
to find a representation of functions of ©, with a>0, B>0. The representation ~ 
in question is 


1 


a 
Ue.) PS Ome mn ele ny) Pte @ am alae 
and may be re-written as 
(4.1) Oey) EY Y= ae pea ig) aaah 


where A is that part of the rectangle -1<6<1, —1<7<1 lying in the half- 
plane x +yy=O0. An alternative form is 


1 il 
(4.2) U(x”) =4S f 1-2) (1 9)" (Ex tn y|) dé dy. 
=i 
Let p(x) be twice continuously differentiable for x=0, let p'(0)=0, and let 
“>0,B>0. Then (4.1) or (4.2) represent a function in C,, in particular, a solution 
of (4.1). Moreover, 


(4.3) F(x) = U(x, 0) =O (1 — 1 (En) a8, 
0 
1 

(4.4) Gly) = U0, 9) = SSE f(t — 08 pny) an. 


(4.5) Bee 
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and it follows from the known identity!° 
Leg patos, pl 40 1 BS &, 


that F and G satisfy (1.5). Furthermore, (4.5) also shows that the function » 
involved in the representation (4.1) may be obtained from F(x) by differentiation 
of fractional order. 

Not every function in ©, possesses a representation of the form contemplated 
here. Since p(x) is twice continuously differentiable and «>0, F(x) will at least 
be twice continuously differentiable; but we know (I §4) that at least for 
B>a>z, this will not be true for all U in @,,. Of course, (4.1) represents a 
function in ©,, under more general conditions than those contemplated here. 

The question as to the most general conditions under which (4.1) will represent 
a function of ©,, is an interesting one, and we hope to return to this point in a 
subsequent publication. 


§ 5. Application of the Mellin transformation 
We shall now transform our partial differential equation to polar coordinates, 
y and #, and apply the Mellin transformation with respect to 7. In doing so we 
Shall assume that all infinite integrals converge, differentiations under the 
integral sign are permissible, efc., and will not discuss precise conditions under 
which these assumptions are tenable. 
Het Us Set 


(5.1) GIN ee St Vit, SOS t= cosh = 


n (4.1). After multiplication by v*, we obtain the partial differential equation 


eU 2 aU 22U 
gertenymy Gore a 
pe?) aU adglous 
(2B —20 + 2t—1) 7 + 2(20— 2at-+-2pt +1 — 844-88) ——=0. 
We call 
(5.3) u(s, t) ee Udr 


the Mellin transform of U; we assume that Re s and U are such that this integral 
converges and the following operations are justifiable. By integration by parts 


CO lo) 
[Soar =[r Ue —s [FU dr =— su 
if « 
0 0 
[o.e} lee) CO 
‘ Sie ; 
[eh yedr=— (+1) le = 7 —=S(S = 14) [re U dr =s(s +1) 
Y A “ 
0 0 0 
loo) (oe) 
perl ae, ae | seo ey a 
ie po er mak eae 
0 0 
co 
Hie 1 eU y 
y y= — St. 
Or at 


10 KoBER, H.: loc. cat., Theorem 4. 
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Substituting all this in (5.2) we obtain for w the ordinary differential equation 


BE 
oe 


p+4 Se a 
t 


CO ersae partners |= Bie— 1) E—¥) 


This ordinary differential equation has four singular points: t=O and ¢=14, 
corresponding to the singular lines x =0 and y=0 of the partial differential 
equation (1.1), t=4, corresponding to the characteristic ~=y through the 
singular point x =y=0, and f=co. The last of these singular points has no 
connection with the partial differential equation, and it turns out to be a removable 
singular point. Setting 

WS te Diss 


we find that v satisfies the differential equation 


le ie eee St fat a4 (95h Te 20) so ey 


(5.5) i +( i ee he (ee: ' 822(t—1) (¢—3) 


The last differential equation has only three singular points, the point at 
infinity being an ordinary point of the equation. Moreover, each of the three 
singular points is a regular singular point, and the differential equation is of 
the hypergeometric type. Determination of the exponents belonging to the 
singular points leads to the Riemann symbol” 


0 1 > 
ng a 0 Os 
s = 1 4 : jee cer ea 
By a linear transformation of the Riemann P-symbol, 
: ee 1 
(5.6) w=t-ty=ze-ts Pp) O = 0 eet 
Grea > + Lime +B-—s 


and either by the symmetry properties of (5.4) when ¢ and 1—¢ are interchanged, 
or by a further linear transformation of the P-symbol we also have 


0) oo 4 
ig oe t 
(5.7) u=(1—f- Pp) 9 A 0 ree 
1 Al 
Fo te eat ee ae 


In (5.4), s is any complex number with a sufficiently large positive real part, 
and we are seeking solutions of this differential equation which are continuously 
differentiable in the closed interval O<t< 1, and twice continuously differentiable 


4 Corson, E. T.: Theory of functions of a complex variable. Oxford 1935, § 10.2 
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in the open interval O0<¢<1. Moreover, with the notation (1.2) we have 


a 


(5.8) ust) =f(6), u(s, 0) =g(s), 
where ; 
(5.9) i) =P PO) ar, 96) = fre w)ar. 


Now, ¢ takes the value $ when #=4z7, and t=# is a singular point of the 
differential equation. The only solution of (5.4) which is regular at ¢=+4 is a 
multiple of 


eR (, <5 LS tap; aes =i 
and this solution is certainly not continuously differentiable at t=0 if a>—H4. 
Thus we shall not be able to use a single analytic solution of (5.4) throughout 
our interval but rather two distinct analytic solutions in the two intervals 
OStS< 4 and $=tS1, “matching’”’ these solutions at :=#. 

There is a solution of (5.4) which is continuously differentiable at t=0. If 
a> —¥, this solution is unique up to a factor, which is independent of ¢; this 


factor is determined by the second equation (5.8), and the solution is 


(5.10) u/s, ) =¢(s) (=) HFS Nestea ag a S3 a OSt<s. 


2 


Similarly, if B>—4 


65.11) u(sd=f(rhR(S, S*-a p+t;4S), f<est 


t 2 
is the only solution of (5.4) which is continuously differentiable at ¢=1 and 
satisfies the first condition (5.8). 


If «+ >0 and we choose s in the strip 0<Res<a-+f, both functions 
(5.10) and (5.11) will be continuous at ¢ =}. They will have the same value there if 


HS) F(Z, S*—a; B+: 1) =e F(S, “S48: a+s, 4). 


De ay Reo 


Using Gauss’ formula! for the value of the hypergeometric series at z=1, we 
can re-write this condition as 


2 


eee ROE es tues oe Ei 
r(p+->*) r(a+8- 5) r(a+) rat 6-5) 


r(a+*>*)r(b+5 
(5.12) ae -< ee) 
rat 5) r(B+ 2 


r(p+5)ret+B—s) rat >)r+B—s) 


or 


f(s), 


and by Koper’s theorem on Mellin transforms of integrals of fractional order?®, 
this is equivalent to (1.5). Thus, if (1.4) is satisfied and «+f>0, (5.10) and 


12 Copson, E. T.: Complex variable, § 10.32. 
13 Koper, H.: loc. cit., section 4. 
Arch. Rational Mech. Anal., Vol. 2 6a 


, oy ae q he - + ‘ ie pais 
an Corsow & ik Eeoeree we Ne iain ei ‘Two Sing gul ya 2 


ae . e 


ea)" 


, - *s th. 
%: te. a continuous function in vies which is analytic, and satietied 


Pt ae 
6. 4), except possibly at ¢= 
- To investigate the dtierant ability of u at t=4, we use the ee: for the 


analytical continuation of the hypergeometric series, together with linear bitin | 


formations of that series, and obtain ee, | 4 


rat a)ete-9 Pe ee hear at) 
r(a+ oa (a+ i see 5 ; 


Dag ee ‘ _ 


u(s, t)=g(s 


ee 


Paes se 
_ 


The first series appearing in these two expressions of w in the neighbourhood ~ 
of t = are identical, and we see that w will be continuous at t=4 provided (5.12) — 
is satisfied and 0< Re s<a +, continuously differentiable if O0< Re s<a«+f—41, 
and twice continuously differentiable provided 0< Re s<a% + —2: in the latter 
case u also satisfies the differential equation at t=4. 

The formulae (5.10) and (5.11), in conjunction with the convolution theorem 

for Mellin transforms could be used to represent U(x, y) explicitly, in terms of 
F(x) if OS yS*x, and in terms of G(y) if O<x<y. Since explicit formulae have 
already been given in I, we shall not pursue the matter further here. 


14 Corson, Er. T.: Complex variable, § 10.33. 
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